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APPROXIM AT ION FOR SMALL Y -1 
II* A- Bethe 

£.1 GENERAL PROCEDURE 

The solution given in Chapter 2 is only valid for an exact point source 
explosion, for constant $ , for constant undisturbed density of the medium 
and for very high shock pressures. It is very desirable to find a method 
which permits- the treatment of somewhat more general shook wave problems and 
thereby comes closer to describing a real s hoc k wave. The clue to such a 
method is found in the very peculiar nature of the point source eolation of 
Taylor and von Neumann. It is characteristic for that solution that the den- 
sity is extremely low in the inner regions and is high only in the irmnediate 
neighborhood of the shock front. Similarly, the pressure is almost exactly 
constant inside a radiis of about .8 of the radius of the shook wave. 

It is particularly the first of these facts that is relevant for construct- 
ing a mora general method. The physical situation is that the material behind 
the shock moves outward with a high velocity Therefore the material streams 
away from the center of the shock wave and creates a high vacuum near the 
center. The absence of any appreciable amount of material, together with the 
moderate 3ite of the accelerations, immediately leads to the conclusion that 
the pressure must be very nearly constant in the region of low density. It 
is intereeting to note that the pressure in that region is by no means sero, 

f 

but is almost l/2 of the pressure at the shock front. 
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pronounced for values of 
the specific heat ratio close to 1* It is well known that the density at 
the shock increases by a factor 


evacuation of the region near. 




. X-*-- 

_/°o T-l 


This becomes infinite as X approaches unity. Therefore, for / near 1 the 

assumption that all material is concentrated near the shock front becomes 

«*• 

more and more valid* The density near the oenter can be shown to behave as 

r 3 /( y~V' 


The idea of the method proposed hsre, is to males rspeated use of the 
fact that the material is concentrated near the shook front. As a consequence 
of this faet, the velocity, of nearly all the material will be the same as 
the velocity of the material directly behind the front. Moreover, if X is 
near 1, the material velocity behind the front is very nearly equal to the 
shock velocity itself; the two quantities differ only by a factor 2 / ( +\) • 

The acceleration of almost all the material* is then equal to the acceleration 
of the shook wave; knowing the acceleration one can calculate the pressure 
distribution in terms of the materiel coordinate, i*e*, the amount of air 
inside a given radius. This calculation again is facilitated by the fact 
that nearly all the material is at the shock front and therefore has the same 
position in space (Eulerian coordinate^* 

» The procedure followed is than simply this* Tfe start from the assumption 
that all) material is concentrated at the shook front* Ws obtain the pressure 
distribution* From the relation between pressure and density along an adi- 
abatic, we can obtain the density of each material element if we know its 
pressure at the present time as well as when it was first hit by the shock* 

By integration of the density we’jo&tT ■Sherfifiaft* a'more accurate value for the 


I?* 


j j < 
? 4 


.1 .* / 
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position of each mass elemsn1?|* P^Pf'kaT.^d be repanted if required; 

it would then lend to a power series in powers of if - 1 . 

The method leads directly to a relation between the shook acceleration, 
the shook pressure and the internal pressure near the center of the shock 
wave* In order to obtain a differential eouation for the position of the 
shock as a function of time, we have to use two additional facts* One is 
the Hugoniot relation between shock pressure and shock velocity. The other 
is energy conservation in some forms in some applications such as that to 
the point source solution itself, we may use the conservation of the total 
energy which requires that the shock pressure decreases inversely as the cube 
of the shock radius (similarity lawJ * On the other hand, if there is a cen- 
tral isothermal sphere as described in the last chapter, no similarity law 
holds, but we may consider the adiabatio expansion of the isothermal sphere 
and thus determine the decrease of the central pressure as a function of the 
radius of the isothermal sphere* If we wish to apply the method to the ease 
of variable If without isothermal sphere we may again use the conservation 
of total energy but in this case the pressure will not be simply proportional 
to l/Y 3 . 


As has already beer, indicated, the applications cf the method are very 
numerous. The case cf not very high shock pressures can also he included; 
in t.hiv case the density behind the shock wave dees not have the limiting 
value of Equation (1) but depends itself on the shock pressure. This does 
not prevent the application of our method .is long as the density increase 
at the shock is still very- large so that most of the material is still near 
the shock front. 


The only limitations of the method are its moderate accuracy and the 

possible comp? lent ions of the numerical work. The accuracy seems satis fa c- 

• » • • • • • * • * 

v* •»••***• 

tory up to a about l.i. For the.*bni r.ll.s oiSrc a direct, comparison with the 
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exact solution is possible. ; .<*. •„ 

5.2 GENKHAL EQUA.TI QMS 

We shall denote the initial position of an arbitrary mass element by 
r, and the position at time t by R. The position of the shock wave 
will be denoted by Y. The density at time t is denoted by , the init- 
ial density by The pressure is p Or ,t) and the pressure behind the 

shock is p s ( Y ). 

The continuity eouation takes the simple for® 


J> R* d R * * 2 dr 

o 


( 2 ) 


From this we have 


'b R 

b r 


jF 


(3) 


The equation of motion becomes simply 


dfR 

dt* 


1 

J* 


K* 

im yr x *** 

^r 9 


(4) 


The pressure for any given material element is connected with its den- 
sity by the adiabatic law (conservation of energy). The particular adiabat 
to be taken is determined by the condition of the material element after it 
has been hit by the shock. If we assume constant 6 the adiabatic relation 
gives 

X - / 

(r ,t) « t» 5 ( r ) ^Q(r ,t)^ 


(5) 


We shall use this relation'mostly to determine the density from the given 
pressure distribution. Using Squat ion (lj^for^the density behind the shock 
, and the continuity Equation wfe find - 

3 Ml ••• •# 


APPROVED FOR PUBLIC RELEASE 



APPROVED FOR PUBLIC RELEASE 


b R 

b r 


= JSt 

R® 



■ » » * « # 

*> M 

A 

if it 

» 

* H» 

* «*. • 




a 

no 

Hr* 


V *(r,t) 7 


V* 


(c) 


The three conservation laws, (2), (4) and (d), must be supplemented by 
the ffugoniot equations at the shock front which are knOwn to be, thews elves, 
consequences of the dame conservation laws. These relations givfe for the 
density at the shock front the result already quoted in Equation (1), for 


the relation between shock pressure and shock velocity I t , 

7 

M 12 2 

" = 


: rv 


if -a 


'2 


. * 

Z 1 ' " ov 


(7) 


and for the relation between the material velocity behind the shock, R, and 

- f'f) is V - Q a / ' ' V y>r" 

/ Vtf 


the shock velocity, Y > 

' R * 2 Y/(Y ♦ 1) 


. / M- 


( 8 ) 


The problem will now be to solve these eight equations for particular 
cases with the assumption that if is close to 1, Then Equation (4) reduces 
to 


• • 

1 d P ^ Y (9) 

r * £ r ° Y a 


Oh the right hand side of this equation we have used the fact discussed in 

the last section that practically all the jpaterial is very near the shock 

front. Therefore the position R can be identified with the position of 

• • • • 
the shock Y , and the acceleration R with the shock acceleration Y. 

Since the right hand side of Equation (9) is independent of r. It integrates 

immediately to give ’ \ i 

* * 

t>(r,t) a 1»3 (y) ♦ J> .JL- <y 5 - p 5 > (10) 

3 Y* 

If we use the Hugoniot relation (7) and put * 1 in that relation we find 
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H r ,t) 
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? - 6 




( Y 


-r«) 


( 11 ) 


This equation gives the pressure distribution at any time in terms of the pos- 
ition, velocity and acceleration of the shock. 

Of particular interest is the relation between the shock pressure and the 
pressure at the center of the shock wave. This relation is obtained by put- 
ting ■ r. a 0 in I qua t ion (1$). Then we get 


(c , t ) j + Y Y/3 (12) 

The pressure near the center is in general smaller than the pressure at the 

• • 

shock because Y is in general negative. 

It can be 3een that the derivation given here is even more general than 
was stated. In particular, it applies also to a medium which has initially 
non-uniform density. It is only necessary to replace ^ r 5 by the mass en- 
closed in the sphere r (except for the factor 4TT/3). 

From the pressure distribution (11) we can obtain the density or the 
position R using Equation (6). Th(£ remaining problem is now to calculate 
this density distribution explicitly, and to determine the motion of the 
shock wave in particular cases. 

5.3 THE POINT SOURCE 

H 

The simplest application of the general theory developed in the last 
secticn is to a point source explosion. In this case, the theory of 
von Neumann and G. I. Taylor is available for comparison. 

Equation (12) gives a relation between various quantities referring 
to the shock and the pressure at the center of the shock wave. To make any 
further progress we have to use the conservation of total energy in the 
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Table £.3 


# r ‘ 

* 

1 \ 1*1 

1 

1*2 

1.4 1 1.67 j 

2 

> v - 

9 

*600 | *462 

*424 

*366 *306 1 

.261 

( If - 1 i tfyy A 

- 

2*094 j 2*162 

| 2.148 

2.144 2.06 

1.96 j 

With tbs relation 

of internal and 

shock pressure known, we can 

now cal- 


culate the total potential energy content. We know that the potential energy 
per unit volume it p/(^ “l) • We further know from Equation (ll) that the 
pressure is constant and equal to p(o) over the entire region which is near- 
ly free of matter* Moreover, we know that all the matter is concentrated in 
a very thin shell near the shock front* Therefore, with the exception of a 
very sna 11 1 fraction of the volume occupied by the shock wave, the pressure 
ia eoual to the interior pressure* The total energy is then 


i * £T ** Jf<pi * U dti 

s T-r s jrr 

« I^^Ly 5 * 2 = 2JL JLa 

3 T-T 3 


(17) 


( of. Equation (13) ) 

In the last line of Table 6*3, above, we give the exact numerical factor 
in last sxpression in (17), according to calculations of Hirschfelder* It 
is seen that this factor is very close to Zw/3, for all values of ^ up to 
1*4. This is due to a compensation of various errors* The internal pressure 
is actually l*ss than l/Z of the shook pressure, but this is compensated by 
the fact that the pressure near the shock front is higher than the internal 
pressure* Indeed, the ratio of the volume average of the pressure to the 
shock pressure is much closer to l/Z than the corresponding ratic for the 
internal pressure (cf .Equations 31a, 31b). A further error which has been 
Tade in Bquetion (17) is that the factor z/( ^ +l) has been neglected in 
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V - 7 

shook wavs* Since there is no characteristic length, time, or pressure in- 
volved in the problem, the blest wave from a point source explosion must 
obey a similarity law as has been pointed out by Taylor and von Neumann* 

In other words, the pressure distribution will always have the same formj 
only the peak pressure and the scale of the spatial distribution will change 
as the shock wave moves out. Now the energy is mainly potential energy (1) 

W ; 

This assumption is not necessary for the validity of the following 
equations* 

/ 

if 'i is close to 1* the potential energy per unit volume is p/ ( ^ -l) 

and therefore the total potential energy will be proportional to p g Y®/( -l) 

Therefore p 8 and Y 2 (of* Equation (7)) will be inversely proportional 
to r • This gives immediately the equation 

Y 2 - AY“ 3 (13) 

where A la a constant related to the total energy. Integration gives 

Y 9 (s/2) 2 / 6 A 1 ^ 6 (14) 

and differentiation gives 

Y Y - -(3/2) Y 2 (16) 

Inserting this in Equation (12/ we find immediately 

*(0,t) 1 ;2 1 l=s 

~ ~ '2 ^ T 

Therefore in the limit of V close to 1, the internal pressure is just l/2 
of the shock pressure* This can be compared with the numerical result of 
von Neumann's theory which gives the following values for the ratio of in- 
ternal pressure to shock pressure* 
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Equation (7). On the other hand, the kinetic energy has been neglected* 
This kinetic energy is very nearly equal to 


E 


kin 


2 it 


A 


_3 *2 
T Y 


(X -1) E 


(18) 


because all the material moires almost with the shock velocity Y* It is 
seen that this kinetic energy is small compared with the potential energy 
by a factor H -1# this justifies our neglect df the kinetic energy along 
with a large number of other quantities of the relative, order 2^-1* It is, 
of course, only an accident that there is almost exact compensation of all 
these neglected terms up to values of X as high as 5/3* 

Ws can now use our result to obtain the density distribution of the 
matter behind the shock front* We need only apply Equation* (l) and (ll) 
to (16> and find 

' ^ ^ 

. *4. (jiuLEi) .> 

7TT V r,Y) 1 T*1 \ 


A 

s 


t 

x(l*x) 


■) 


with 


r3/ y 3 


( 10 ) 


(19a) 


Setting also 


Equation (6) beootnes 


JL~ 

dx 


3/3 

R A , 


# ll ( 2 

♦1 v x(l+x) / 


TZ 


i/if 


(19b) 


( 20 ) 


to integrate this eqi*aticn, it is convenient to distinguish two oases; 
(l) If x is not too small, more precisely for 


( 20 «) 
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y 


we may cone ider the exponent in (2o) as equal to 1 sine* 

Y is assumed email* Then, negleoting quantities of order 
2 

( if -l) , integration of (2 g) gives 


- 1 - (tf -1> log i-JLS * 1 - (ST -l) log -- * -f ■ (2l) 

Y° 2 x 2 P 3 

(2) If 

x « 1 (21a) 


we may pegleot x compared to 1 • Then, neglecting quantities 
of relative order it -1, we get 


dy • ( Y -l) x dx 

y . , A 


( 22 ) 


where A is a oonstaht* The regions defined by (20a) and (21a) 
overlap very considerably. Comparing (2l) and (22) we find 

that '** 



(22a) 


negleoting a a me 11 term of order Y-l« This value of A will 
make (22) sensible for small valuds of x. Inserting (10fc), 
(19b) , we get 

R - Y 1 ^ r 1 “ (23) 

m 

or 

r/y . (r/Y) (a ' 1,/ ' (23.) 

From the position of any point we can deduce the velceity by a simple 

V 

differentiation with respect to time. In this process, the material coordinate 

i' a 

r should be kept constant* Equation (23) gives for the dfiiterial velooity 
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(neglecting terms of order Y -l) 







( 24 ) 


Over most of the volume the materiel velooity is nearly linear in R whioh 

* 

is borne out by the numerical integrations of the exact solution (see Charter 
2). Over most of the mess the material velooity is nearly equal to the veloc- 
ity of- the shock wave. 


6.4 COMPARISON OF THE POIRT SOURCE RESULTS WITH THE EXACT SOLUTION . 

The results c&tained in the last section can be oompared with the exact 
solution described in Chapter 2. The results of that chapter can very easily 
be applied to the special oase when Y is very nearly 1. 

In going to this limit one should keep the exponent of 9 correct because 
this quantity goes from 0 to 1, and if it is close to 0 a factor 9 * -1 
will matter. In all other factors the base of the power becomes (0 ♦ l)/2 
in the limit Y «1, whioh goes over the renge from 1/2 to 1 and therefore 

t 

never becomes very small. Consequently 5 -1 may be neglected in the exponent 
of those other factors axoept if higher accuracy is desired. 

Regleoting small quantities in this manner, Equation (2. 37) redoes s to 

r /Y , 2 . $ 1 6 l ^ Z (26) 

2 and © being the notations used in Chapter 2. A ^ 

(26) may be rewritten 

9 m 2 3 . <r/Y> 3 (26a) 

Eouation (2.38) becomes then 

l* if -1 

b/y - f » § z '/*i . ( r / Y )Y^T ( Z6 ) 

This result for the Eulerian position i» identical with that obtained from 

i 
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our approximate theory in Iguatlon (25a) • A more accurate •valuation, keeping 
terns of relative order -1 throughout, gives 

2-1 

a» 

(26ft) 


(rrr) 


W7T 


For the pressure we find fron iquatien (2*41) in the limit Y »1 the re- 


suit 


# ♦ 1 
** r~ 


* T*» 


[‘ •«»] 


( 27 ) 


Thig again Is i&ehtioal wifch the result of our approximate theory given la 

' •*** ' ’ 


Bquation (ll) end Equation (16) « Again, a more accurate evaluation, neglect- 

2 

ing only terns of relative order ( 2-1/ , gives t 




•|-i [l*«{ U —1 ♦f( r i- T -jj[ 


(28) 


with the abbreviation 


£ m H -1 


(28a) 


Of particular interest it the relation between the total energy and 
such quantities as the shook pressure* We shell therefore calculate this 
relation including: terms of relative order Y - 1* First of ell, we shell 
oaloulate the potential (heat) energy * 


S pot. * *V *‘ dl */(*-!) 


ly 


• ¥ hr- [ d(F5) TT 


(29) 
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Using (26a) end (28*), we have 


a(p ! ) : d (e fM 


(* ) 

[ d 1 - " i~\' ] 6 T-) 


(29a) 


neglecting^® terms of, relative order ^ 2 * Inserting (29a) and i 5 from 


7 

It ie permissible to set the factor 0 which should appear in the 
, •> second term in the square bracket, equal to one; the error in (29) is 
only of order of^‘ • 


(28) into (29) i we get 

d (F°) 


I X 2 £ * 45 


■t [ d < 9<r > 


:] 


e+n 


["‘(htP+tfr -*)] 


(30) 


r 

This integral can be evaluated very easily* ffls note that 0° changes from 
0 to 1 at very small values of 0 so that in first approximation for this 
part of the integral, the integrand should be taken at 0 ■ 0* (This 

corresponds to the physical fact that most of the material is near the shook 
front, F 3 » 6^ becomes close to 1 already for relatively email value* 
of 0 or of the material coordinate Z - 0^ 3 ) . Evaluation of (SO) g|v** 


1 ♦ * <A± + 1 - I) +£-£ 

2 2 


(3Qa> 


or 


B pot. * -y- 


T y* f [*+£< | -A> 2 )] 


(3l) 


This result, except for the last factor, is identical with the result 
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of our approximate theory, Iquation (17). The laet factor ia seen to differ 
only very elightly from 1, the factor of $ being only -0.2. It is of 
some interest to define the average pressure (volume average); this is ac- 
cording to (3l) » 


: *,<> * £ *’• [ C| -&> s )] : } *. f 1 - •«»*' 

This T»y be compared with the central pressure (of. (2$) 

1 +• <T'( 1.- bn s 1, ^ 1 - .886 


(31a> 


♦(0) 


1 K 
y • 


(31b> 


The average pressure is, of oourae, higher than the central pressure; it dif- 
fers from it only in the order as is to be expected; and it is muoh closer 
to one-half the shook pressure than the central pressure is. 

How let us calculate the kinetic energy. According to (2.4g), the ratio 
of kinetic to potential energy in any mass element is 8 , therefore 





y 3 * a 

Jf 1 d(F 8 ) 

9 

"a -t 

2 IT 

~T~ 

y-i 

4 

2ir 

1 

* ! v 



(32) 


The simplifications in this integral, l*e» negleot of the last square bracket 
ih (3o) and replacement of d (0^) by &&G/&, are possible because of the 
fact'3r 0 in th« integrand. This also makes the integral of order S • 

The result (32) agrees with that of ths approximate theory, (18). 

Adding (3l) and (32), we find for the total energy 


1 -p H 




*e T* 

Y=r 


[l + S(§ -U 2)] 

fit , D 1 


( 33 ) 
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This gives the shock pressure as a function of the radius including terms 
of relative order 2f -1 which will be useful for the calculation of the waste 
energy* We may also replace by the shook velocity Y according to (?) j 


jA o *)L 


i' [»+/ 


(33.) 


Here again the correction faotor in the square bracket differs only slightly 
from 1, in agreement with the nuaerioal results reported in Table 6*5* 

A further quantity of interest is ^r/ &r • dF/da for which Iquation 
(2*i7) gives the result 


p R . = * j f. *- i 

77T w ' 


7 TT 


Ffom this expression or direotly from Equation (2*39) we can find the density 


which turns out to be 


4i' 4^ ® *** 


We can also express this density in terms of the Eulerian position in which 


case we get from (26) 

_ 3/(**l) 

f* s ^■+1/ , R , \ 1 ( 1 + Z^ 

“5* yrr-(W r (I+I) 


(35a) 


This equation shows that the density becomes extresfcly low for all points 
away from the shock front even if they are only moderately close to the 
center of the explosion. This is in agreement with our basic assumption 
that most of the material is concentrated near the shock front* 

Finally oombin-'ng (2*38) and (2*40^ We find in the limit If' - 1 
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However, from the similarity solution wo know that 

_ */* 

T » a t 


where « is • constant, and therefore 


jt 

¥ 


T 

Y 


(36b) 


This result is again identical with ths result of our approximate theory 
given in Equation (24.) . 

WO see, therefore, that our approximate solution is identical with ths 
limit of ths sxact solution of the point sourcs for # ■ 1 if tsrms of thft 
relative order -1 *ro oonaistontly noglootod* 


S.5 THE CASK OF THE ISOTHERMAL SPHERE 
T/*s shell now consider the somewhat more complicated problem of the pro- 

*' . w 

duction of a shook wavs by e qjpher# which is initially hasted to a high uni- 
fora teaperature and then exof^s pressure on its surroundings* The relevance, 
of this problem of the ieotherqial sphere has been disoussed in Chapters 1 and 
4 and is connected with the great influence of energy transport by radiation* 
The problem new no longer permits the ‘application of similarity argument** 
For this reason we can no longer use the conservation of total anergy to ad- 
vantage* Instead of this we can now assume adiabatic expansion of the iso- 
thermal aphere* This is completely equivalent to an application of the energy 
conservation law because the adiabatio lap itself is based on the assumption 
that there is no energy transport out of the isothermal sphere* 

Let ua assume that the materiel coordinate of the surface of the iso- 
thermal sphere is r 0 * The initial position of this surface is then equal 
to . r 0 * At a later tine When the isothermal aphere has expanded to R a its 
average density has decreased by a factor (r^R^ 3 • if we assume that the 

' - * . 
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density and pressure in the isothermal sphere ere uniform* the pressure will 


be equal to 


K • * 'i**7*o ) 


3 r 


(37) 


where P is the initial pressure in the isothermal sphere whioh la related 
to the total energy by the equation 

5 






( 37 a ) 


Ws shall now prooeed in two steps* First of all we shall ooniider the 
case when the radius of the shook ware is not extremely great compared with 
the initial radius of the isothermal sphere, more precisely 

T/r 0 «ft 

We shall show that with this assumption the solution with an isothermal sphere 

¥ 

approaches the point source solution as f^r 0 increases* Afterwards, for large 
values of l/r 0 where (37b) is not valid, the approximations used in the first 
part of our calculation will break downj but we then use the results of 
the last section to obtain Bg/Y for the surface of the isothermal sphere, 

and this will enable ua to solve the- problem for the ease of large values 

of Y/r 0 . 

(l) Case I: Y/r 0 Moderate 

If Y/r 0 is not very great we oan replace the exponent 3 S' 

in Equation (37) by 3* At the same time we oan use our general assumption 

that praotioally all material is close to the shook frost and that, therefor*, 
R 0 ia very nearly equal to Y* With this approximation we find from Squation 
( 12 > 




P r. 


-3 
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Thia aquation can ba integrated without difficulty by catting 

y ’ # p It) » 

and 

Then 

V . i if. 

2 d T 

Then Equation (9ft) become a 

^ 4 i 1 - AT 

Thia can be integrated and give a 
<p Y* s 2AT ! +B 



(39a) 


(39b) 


(to) 


(4l) 


where B It t new oonetant* 

Hew the initial condition ie, for Y 


r o * 


a - 12 t 

r “ 

o 


-?r = 1 


or* with (39) and (39a) : 


(42> 


<f> ( r o) 


(«•) 


ao that 


w. 


B « - A r. 


(42b) 


and finally 


s^ c Y a fk> A ( iy - £$7- ) 
sor-l) 1 / r 0 3 \ 

^ “ r? 


2 TT 


Y& 
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Equation (43) show’s that for large values of Y/r 0 , the shook pressure 
£s approaches the value given in Equation (17), i.e. the value corresponding 
to the point source solution* This result appears rather important because 
it shcws?(l) that the point source solution is stable and is approached even 
under conditions whioh initially deviate strongly from those assumed in the 
point source solution, and ( 2 ) that the simple and -well-known point source 
solution oan be used at late times for our problem including the isothermal 
sphere* 


(2) Cane II » y/r 0 Large 

We know from the disouselon of Case I that our solution ap- 
proaches the point souroe solution es soon es y/r 0 ^ 1* We oan than usa 
Equation (23) for the position R^ of the surface of the isothermal ap Hare, 
and obtain, neglecting terms of second order in i^-li 



This is the sane expression whioh n» used in Case I and which we then justi- 
fied simply by neglecting y in the exponent of (26). Therefore the further 
development is identical with that leading to Equation (43). 

We have thus shown that Equation (43) is valid both for small and for 

* 

large expansions of the isothermal sphere* It is possible to derive the den- 
sity distribution, the position end the velooity, as we did in the previous 
section for a point souroe case* However, the analytioel expressions are 
fairly involvad and thare does not seem to be any particular application for 
them* The ratio of the shook pressure to the central pressure in the iso- 
thermal sphere is according to (43) and (36) t 



2 - r 0 v r> 


(«.) 


4 


f 


APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 
V - 20 

It changes from 1 in the very early stages to f 2 in the late stag*#. 

At any reasonably late stage the energy in tha blast Is related to the veloo- 
ity of the shock front in the sane way as for the point source solution. 

j 

6.6 VARIABLE GAMMA 

The theory developed here can be used to solve the problem of a ehock 
wave in a medium with variable Y • The assumption is, of oourae, that -1 
still remains small throughout. We also assume that the shock pressure is 
high enough so that the HUgoniot relations hold in their limiting form. 

We shall make the further simplifying assumption that ^ is a function 
of the entropy only, so that it remains constant for any given mass element 
r as soon as that element has been traversed by the shook. This assumption 
is fairly well fulfilled by air, with the value of )f decreasing from 1.4 
to about 1.2 with increasing entropy, and later on increasing again to 1.67. 
The more general problem in which % is a function both of the entropy and 

the density can also be solved by the seme method, but the algebra becomes 

. \ 

so involved that it seems hardly worth-while to use the present method in- 
stead of direct numerical integration. 

•v ^ 

The pressure distribution Bquation (ll) will still be valid. However, 

• # • 

the relation between Y, Y, and Y will no longer be given by Bquatlon 
(16). We introduce the pressure at the front and the pressure at the center 
of the shock wave separately by writing 

: * -h 

Ho) -f* -£-p ■ <«> . 

where and (. 3 are slowly variable with the shock radius Y. In the case 
of constant Y we have (cf. 17) 
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^ -1 


It is our aim to calculate oC end j& for a given variation of 7f with 
the material coordinate r. 

Using Squat ion (llJ we have 


ITS T 14- < 4r > 

where J) is an abbreviation for Y • Inserting the expression (45/ and re- 
number ing (?) we get the following relation between and ft • 


f) S °<+-r 


d e< 

T T og T 


It will be shown in the following that °< changes very slowly with log Y# 
in fact dc*/d log Y is of the order Y -1 relative to ©< itself. Therefore, 
in our theory in which # -1 is considered as small we have 


Therefore, even with variable ¥ the ratio of the shock pressure to the in- 
ternal pressure is equal to 2. 

In order to find o< for a given function Y (r) we use the fact that 

i 

both the geometrical and the materiel coordinate must be equal to Y at the 
shock front* We shall calculate the geometrical coordinate R as a function 
of r with the help of the density distribution. The required condition 


is then 


R 3 (r 


i 

- Y) - J~ d(r 5 ) 


As in Equation (6), the density of the material element r at the time 

J 

when the shock wave is at Y, is given by 
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l/v 

' TT r ~ ; -i -5srkr ‘( 1+ V) ( t ) Vif 


1/V 


( 51 ) 


In this equation we have made use of the pressure distribution (11) and also 
of (45) and (49). Furthermore,** have put in the exponent equal to 1 

in all these terms where this makes an error of the order V -1. Inserting 
(6l) we find R as a function of r as follows 

r, 3 r - 1 ] / r 3 \ 

/ i(r) 4-) °< lr) ^7") 

° (62) 


•1 


E s - 


No appreciable error is made by neglecting the last factor in this expression 

because it is different from 1 only over a region in R of the order )f -1 

■y 

In fact, only by neglecting this last faetor do we get the correot result for 
constant / , owing to other negleoted terns in our theory* In any casejthis 
last factor ie no different for constant and variable % end ean > therefore, 
not be relevant for the theory of variable % • 

With this simplification we obtain 


^ _ 1 
- ^rrr ✓ 


/ - 

and in particular for r * T we find, after dividing by T 3 

~ (I1 : / rfe, Y ' S 0 H* (r) ■ 


This is the desired equation determining o< (y)» It is seen to be a linear 
integral equation. 
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Let us first consider the esse of constant Y . 
(3‘i) reduces to 


JL 

<*(Y) - f*<(r) 


3(ar-i) 


In this case, B^uation 


(65) 


which is solved by c< > constant, with the value of the constant arbitrary. 
From previous considerations, particularly Equation (17), we know that in 
this case ■ Y ~1* 

Let us now assume that Y is constant and equal to , for all values 
of r up to rj i thenoC -1 in that' region. Further, for 

Y y ri the integral Equation (54) reduces to 

M-i) 


(y) = o< x ( ri /r) 


r dr 3 -3 4 3/jf(r) 

■{ P55TF) T 


. ' (Y(r) -l)®C(r) 


(56) 


In order to solve this equation we proceed in two steps, similar to the cal- 
culations in Section 6.6. In the first stsp we shell consider Y/r^ as 
moderate so that the Y "1 power of the quantity can be considered equal to 
1. In this cafe we shall be able to obtain a general and rather simple dif- 
ferential eouation for o( which can be solved by ouadratures. As a second 
step we shall then admit large value*, of Y/r^; in this Case we shall obtain 
a solution only in the special case of having a ptep-wise variation of / « 
(l) Case I: Y/V^ Moderate 

In this case we may replace the exponent 3/V in Equation 
(66) by 3« Also we shall expand the first term on the right hand side of 
that equation in a tfaylor series. Then Equation (66) reduces to 

£ 

«=>< (Y) -3(/ 1 -l)o/ 1 log (Y/r x ) (G7) 

Y 

3 J~ d log r (Y(r )-l)o<(r) 
r l 
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To solve this integral equation we need only differentiate with respeet to 
Y or better with respect to 

Is 3 log (Y/r^) 


( 68 ) 


Then we obtain 


do( 

nr ' 


o<' 1 (^ 1 .l)4o<(x)(Tf(x)-l) 


(69) 


This equation can be integrated with the result 

°<(x) = e • 


f Ml* )-l) dx* ^ 


1 


dX* C 


/ « 

tyCx^ )-0 dx< 


(60) 


The result (60) as well as the differential Equation (69) show that 
does not have a discontinuity at a point at which Y has one, but only 
do(/dX has a discontinuity. This nay be set even more In evidence by 
solving (69) for small values of X j i.e* for points just beyond the place 
at which Y bejins to change. In this case we find * 

*IX) - [it / dx' MX') J^)] ( 61 ) 

A special case which is of some interest because it is the simples? 
possible model of a subctanoe with changing / , is obtained by assuming 
that Y has the v e lue $2 ? or values of r y r^. In this case (69) 
can be solved explicitly with the result 

(tz-V* i 


<=>< s Vr 1 

T 


m 




(62) 


Fcr small values of X this reduces to 

r 


°< • (^-D 


4 . . . 


1+ <«> * ] (63 > 

Another result which follows from Equation (59) ie that dof/dX is of the 
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order =*< ( ^ -l) • This result has been used above in obtaining the relation 

jS *o<. 

(2) Case lit Y/ r i Large 

We shall consider this problem only in the particularly simple 
case when 2^ hAS the constant value for all values of r ^ r^. In '•* 

this case Equation (56) reduces to 

^ (I) » o^ i ( ri /y) 

Y 


3(^-1) 


/ d(f) ° <(r) 


(64) 


Similarly, as in Case I, we solve this equation by differentiating with respect 
to X* Then we obtain 

s p \ s (*l’* 1 ) 

= ‘ (m." 1 ' °*1 ( ■+ (ig-l)«<(x) 


^ <Ur*) *<»> 


(65) 


The integral in this equation can be expressed by means of Equation (64). 
Neglecting terms of higher order in X -1 we get then 

rip 


(66) 


This equation is actus 1 ly even simpler than the differential Equation (69) 
which we obtained in the approximate theory of Case I* Using the boundary 
condition o(, « / 1 -1 for Y * r^ Equation (66) integrates immediately 
to 


o( r V g -1 - (*-*) (r a /Y) S(i l 1} 


(67) 


For small values of X this reduces to 


o<s (V X > ( 1+3 Vr x ) 
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whleh is identical with the result (63) • 

Equation (67) shows th e t the value of goes gradually from the orig- 
inal value )^-l to the value # 2*1 which corresponds to the new value of 
Y • It is seen, that this asymptotic value is reached only for ertrerrmly 
large values of l/rj. As long as y/ti is moderate, the shock pressure 
is still influenced largely by the previous value of Y instead of by the 
present one* For air in particular we may take * 1.2 and tfg . 1.4. The 
shock pressure in a substance with Y constant * 1.4 should be twice as 
great as that for Y constant and equal to 1.2, for the same radius Y of 
the shook wave and the same energy E* Actually, when the shock pressure 
falls low enough so that Y increases to 1.4 the coefficient c< does not 
immediately increase by a factor 2, but increases very slowly. Physically 
the reason for this is that the interior part of the shock volume still has 
the low value of Y and therefore has a high internal energy for a given 
pressure. Only when the hot gases which possess the low Y fill a small 
part of tie volume included in the shook wave will the Y in the outer re- 
gions of the shook determine the shock pressure. 

6.7 THE WASTE ENERGY 


G« I* Taylor has introduced the concept of the waste energy, i.e. the 
energy which remains in the hot gases traversed by the shock wave after an 
adiabatic expansion to a pressure of 1 atmosphere. The knowledge of this 
waste energy is useful because it permits one to calculate the energy which 
remains available to the shook wave at small overpressures* 

The waste energy can be calculated very simply for the point source 
solution. Let ub consider a material element which is traversed by the 
shock at a shock pressure t> g « When this element has been expanded adia- 
batically to atmospheric pressure p c , its density will be 
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( tt ) 


Var 


/> 


( 69 > 


Its temperature will he fj>/^0' l and its energy content per unit mass: 


£ : 


(70) 


>° 0 (*Vl) 

In calculating the energy content, we here ueed the specific heat 8 t conetent 

S 

pressure j the reason for this is that our final state is obtained from normal 
air by heating it at constant pressure p o to the temperature l> o /jz . Striotly 
speaking, in r~der to get the waste energy, we should subtract from (7o) the 
• expression (X -l) , but we shall confioe our discussion to the case 

when • 

The total energy wasted In the shock wave is then 


Vf - 


4TTyO 


* 2 


d Y 


(71) 


We can now use the relation between shook pressure and radius, (33). We are 
using this fairly exact relation because it will turn out that we have to know 
the waste energy including terms of relative order V -1. We solve (33) for 
Y* inserting the result and (7 q) into (7l) we obtain 


W = C 


e ip \ Y 
\° I Kt>, 


f 


d T>_ 


(72) 


with 


C • 2 (*-l) [ 1- (tf-1) ( | .in 2)] X A&l) (72a) 

(rf-1) [ 1 - {X -1) ( 1 -in 2)1 

A 

neglecting terms of relative order (X ~V • Equation (72) gives the energy 
wasted up to the time when the shock pressure has fallen to the value Kp 0 « 


Equation (72) can be integrated immediately and gives 

y / *i \X -] 

w = C Ej 


■m 


( 73 ) 


= J 1 + (&-l)£n ij E K ” 
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Subtracting this expression from the total energy E we obtain the energy 
which is still available* 

E .ff. -- ‘[‘-wi't] (74) 

It is clear that for small values of X -1 and moderate K, this expression 
is proportional "bo K -1. T! is shows the necessity of knowing the relation 
between p g and 32 up to terms of the order V -lj i.e. of using (33) rather 
than (17). 

It is somewhat problematic what to use for K. Clearly the relation (33) 
will break down at too low values of p g i namely, when the limiting form of the 
Hugoniot relations ceases to be valid* This requires 



Setting 

K = 2 n/ftr'-l) (?6 ) 

then 1 

end using the fact that if -"l it small* Equation (74) reduoes to 

Vf. - E (!M) 1o *?=i <7,) 

Of course the available energy will be further reduoed as the shock pressure 
it reduced closer to atmospheric pressure. In fact, Penney has shown that 
the dissipation of energy continues indefinitely as the shook wave expands 
(see Chapter 8 , Section 8 . 7 ). It is therefore not possible to give any accurate 
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definition of the percentage of energy wasted but it would be necessary to 

I 

specify the shook pressure for which this statement is made* In order to 
obtain the percentage wasted for shook pressures of the order of one atmos- 
phere, it would be necessary to follow the shook wave through the region of 
intermediate pressures which can adequately be done only by numerical msthods 
(see Chapter 7) • 

However, it is clear from our Equation (77) that the energy available 
for the shock wave is smaller the smaller )( -1. The fact that the air has 

a small value cf if at high temperatures leads to increased dissipation of 

_ + 

energy at the shock and, therefore, to a relatively smaller blast wave at large 
distances* This is the main reason why the blast wave from a nuclear explosion 
is less strong at a given distance than from a TNT explosion liberating the 
sane total energy. In the latter ease the temperatures reached are only mod- 
erate, and tha anargy wasted is, therefore, less than for the nuclear explosion. 

The greater wastage of energy for smaller )f is^elso related to the re- 
sult obtained in Section fi*6; namely, that the shock pressure after a sudden 
change of if does not oorreapond to tha new value of if but is rather dose 
to the value for tha original value of )f * This slow variation of tha shock 
pressure is in turn important ^ if one wants to calculate the waste energy for 
a gas with variable . In v fact, if it were not for this gradual change, 
gases might occur in whioh the waste energy would be greater then the total 
energy available which would be obvious nonsense* The change of Q< derived 
in Section 5*6 is just sufficiently gradual to keep the waste energy always 
below tha total available energy* 
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‘ CHAPTER 6 

% 

EFFECT OF VARIABLE DENSITY ON THE PROPAGATION OF THE BUST WAVE 

— — — — — ■ — — - 

K* Fuchs 

6.1 INTRODUCTION 

We h«ve seen in preceding chapters that there is a considerable range 
in -which fairly reliable predictions about the propagation of the blast wave 
may be made* The range extends from somewhat below one million degree shock 
temperature to about 6,000 degrees. Above 1,000,000 degrees the isothermal 
sphere extends up to the shook front. If r Q is the radiue of the shock 
front at the time when the isothermal sphere separates from the shock, than 
the shock pressure at a given radius Y > r D is proportional to (see Equa- 
tion E.43> * 

3 

- (1 - r ° /2Y 3 ! (1> 

i 

and therefore the effect of the isothermal sphere is negligible even. for 
moderate values of Y/r 0 « ' 

Si- 

Below 5,000 degrees the formation of an opaque layer at a slight die** 
tance behind the shook front may lead to absorption of radiation whioh might 
otherwise escape, and therefore radiation transport of energy has a consider- 
able effect on the propagation of the shock* Furthermore, be low 5,000 degrees 
the luminosity comes from some distance behind the shock front and this makes 
photographic observation of the shock front difficult. 

Between about 200,000 degrees and 20,000 degrees shock temperature the 
ionisation of the L-«leetrons proceeds and within this range the variation 
of is not very pronounced. This , therefore, appears to be the most useful 

VI -1 
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range. 

The corresponding shock pressures ere 


T 1,000,000° 200,000° 20,000° 

p s 360,000 60,000 1,600 


5,000° 

200 Atmospheres 


Observing thet the shook radius is roughly proportional to the inverse oube 

* - 

root of the shock pressure, this gives a fairly comfortable range of useful 
observation* 

However, there are some factors which limit the range. One factor is 
the height at which the explosion takes place. When the shock reaches the 
ground a reflected shook goes back and only those pertc of the shock sphere 
which have not been reached by the reflected shock, can be compared directly 
with the theory. This was particularly serious at Trinity, but would also 

‘present some limitation in future tests, since it is impracticable to raise 

•» 

the gadget to a great height without interfering with other experiments. 

At Trinity the gadget was set off at a height of 30 meters. 

The other limitation arises from the fact that initially the propaga- 
tion of the shook is affected by the material of which the gadget is composed 
Although the dimensions of the gadget are rather small, the effect persists 
over a considerable distance, since it is th« mass in the gadget Compared to 
the awss of air engulfed by the shock which matters. 

In particular, with a view toward an evaluation of the energy release 
in the Trinity test, ws shall attempt in this Chapter to get over the second 
limitation. Wfe consider a blast wave originated by a point source, travel- 
ling through material of variable density. We shall make some simplifying 
assumptions : 

(l) The density is supposed to depend on the radius only, and the 

variation of density is assumed to be continuous. We disregard. 
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therefore, Any •symmetry in the construction of the gadget 
and the exact details of the transmission of the shock from 

9 

the gadget into air are neglected. Neither of these two 
factors can have appreciable effect at a sufficiently large 

i 

distance. 

(2) We assume that - 1 is small compared to 1 and use the 
essential ideas on which the approximation of small Y -1, 
discussed in Chapter E, is based. Insofar as this approxi- 
mation can be compared with exact calculations of a point 
source it was found to be very good, even if Y is s/3» 
However, it does not necessarily follow from this fact that 
the approximation is equally good in more general cases. At 
present we have no means of estimating the error. 

(3) We assume that Y is constant. This is not a bad assumption 
sinoe we are mainly interested in the region from about 
200,000 degrees to 20,000 degrees, where Y does not vary 
too much. 

6.2 METHOD OF ESTIMATING ENERGY RELEASE BY OBSERVATION OF THE 
SHOCK RADIUS 

Be’fore proceeding to the analysis of the problem with variable density, 
let us consider the application of the method of estimating the energy re- 
lease in the simplest case when the similarity solution for constant density 
holds. The derivation given below is due to Bethe* 

If -1 is smell, the kinetic energy may be negleoted and the total 
energy is given by 

E * | — E— dv 

J ar-i 
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as pointed cut in Chapter 5, the pressure is constant and equal to l/2 the 
shock pressure p 8 over the greater part of the volume* Hence (l) can be 
replaeed by 


E 


1 P« 4ir 3 

T T^i T Y 


(3) 


where Y ie the shook radius. Using the Hugoniot condition for strong 
shocks 


Ps 


Yil * * 


(4) 


we obtain for ^ • 1 


2 ir J> n U 2 Y 3 

tt V=r~ * 


( 6 ) 


Here U is the shook velocity and ~P Q the normal density of air. 
Since E is constant, we can integrate (6) and find 


• f{ 


\T * /2 


If the finite value of -1 is taken; into account, the factor 2ir/3 
should be replaced by a constant B (i ) which differs but little from 
2tr/3. B(j^ ) ie the quantity tabulated in the last line of Table 8.3, 


( 6 ) 


Chapter 8* 

S/2 

If, therefore, we plot the observed values of Y • against the time 
t, they should lie on a straight line and from the slope of the line we can 
immediately derive the total energy. 


6.3 INTEGRATION OF THE EQUATIONS OF VOTIOH. 


We shall now turn to the case of variable density. If -1 is small 
the kinetic energy is small cor, pared to the internal energy. The internal 
energy per unit volume is 
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p/<tf-l) (?) 


Now p ia practically constant over the larger part of the volume of the 
shocked sphere; it varies appreciably only in the region where most of the 
mass ia accumulated, i.e* near the shook front* If we integrate (?) over 
the volume, we can disregard the small shell of variable pressure near the 
shock front and identify p with the- pressure p(e,t) at the center of 
the sphere. Then the total energy E is given by 

E . (8) 

Using Euler variables R, and Lagrange variables r, the equation of motion 

ia 


( r) r 2 R + R S T s 0 


(9) 


where ^ is the original density of the mass element which was at the radius 
r, before the shock reached it. Integration of (9) yields 


p(r,t) 


I (IX 


J* (r) r dr 


( 10 ) 


where r/r 2 is to be considered as funotion of r for fixed time t. p g 
is the shock pressure and Y the shock radius. Now, nearly all mass is 
near the shock front* Hence r/r 2 is practically identical with y/y 2 . 

This will no longer be true if r becomes very small, but then the contri- 

"P 

« 

bution to the integral (lo) will be 6mall. Hence we may write approximately 

Y 

p(r,t) « P ® + Y 2 ^ 1-2 dr (ll) 

r 

In particular at the centre 


p(o,t) 


Ps 


♦ ' 


12 
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where 


m(y) 


. r J» <rJ t 


dr 


is the total mass per unit solid angle within the shock radius< 
The shock pressure is given "by the Hugonict condition 


7TT * ^ 


and U is the shook Yelocity 


U 


Since we assumed Y near 1* we may omit the factor 2/(2f *l) 
Considering U as function of f, we have 


U 


dU 

dY 


Substituting (14) and (16) into (12) one finds 


p(o,t) - / U Z 


U . d£ 
Y 2 dY 


i i. iA 2 ) 

2V.Y 2 dY 


Combining this equation with (8), we find 


dY 


(uV 5 ) 


— (Y -l) 1 E 

2 n Y 


and therefore 


U" 


!r <*- lJ b' j 


«(r) 


dY 


The lower limit, in the integral is derived as follows! 

If Y is small, J* may be considered constant ar# equal to its 
the center* Then (19) reduces to 


IT 


— (Y-l) JE-3 

2tt j> Y° 


(13) 

(14) 

(16) 

in (14). 

(16) 

(17) 

(ie) 

(19) 

value at 

(20) 
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which is the Squat ion (6) derived for a point sourea and constant density, 
as it should be* If a oonstant of integration were added to the integra- 
tion in (19) we would not obtain the oorrect behavior for small shock radii 

2 

In add it ion, the kinetic energy 4K MU would be oorae infinite for Y -*■ 0 
contrary to the assumptions made* 

6*4 EFFECT OF VARIABLE DENSITY NEAR THE CENTER ON THE AIR SHOCK 


We assume next that the density has an arbitrary distribution f (Y) 
up ta a certain radius Y* and constant density f 0 beyond. Vie wish to 
know the propagation of the shock, after it has reaohed the region of con- 
stant density* 

If Y > r», we have 


SV i Y 2 dY mi f 9 (Y 3 ♦ Y 0 5 ) 


(2l) 




where Y 0 Is the radius of a sphere occupied by the excess material at the 
center, if it were spread out at the density f Q i i*e» 




o Y 0 


^ £/>(y) " ./ej y2 « 


( 22 ) 


Similarly 

'(Y* 


h 


dY 


7 A 


• Jr l ( 


f(r> - j>oj Y * dY ( 23) 


Partial integration yields 


I 


KC lj. r . 1 p y 3 

Y ^ m 9 Jo ^ * S 


s fo (l» 


Y - in Y \ (24) 


} 


where is the logarithmic average radius of the excess material 


Y’ 


Jfa Y - |jb Y(/(T> r fo) Y 2 &tj J ( fix* - fo) Y 2 dY 

Substitution of (2l) and (24) into (19) yields 
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sEcaf-u * 

2*>o y3 


f (y) 2 



whe re 

f(T)* . 

7 1 ♦ (V Y) J 2 

If f • 1, the formula (26) reduces to the equation (5) for constant 

density* f , therefore , is the correction factor by which the shock velocity 
is changed* 

It is of interest that the factor f depends strongly on the radius 
Y 0 * The latter depends only on the total mass of excess material at the 
center, but not on its distribution* This is rattier fortunate, because the 

i 

distribution of matter in the gadget is changed during the implosion and 
it would not be simple to calculate the correct distribution* The radius 
Y on the other hand depends on the distribution of matter* However, since 
it is only a logarithmic average, it is net very sensitive to small errors 
in the assumed distribution of matter* 

For the Trinity test the gadget was located on a tower and we are most 
interested in the exper.sion of the shock before it hit the ground* The ef- 
fect of the tower on the average shock radius may be bracketed between two 
limits. On the one hand we may neglect it* On the other hand, we may as- 
sume that the matter in the tower is spread with spherical symmetry around 
the gadget* This leads to the following problems 

The density distribution is arbitrary up to a radius Y' as before. 
Beyond that radius the density is given by 



y i 2 Ay for Y > Y' (23) 

The formulae can be evaluated for this case in a manner similar to that em- 

p loved above j one fines for the correction factor f 
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1 + 


3 Y, 


jin (Y/Y) + 



Y* 


Y 0 


( 1+in 




? 


( 29 ) 


la order to be able to neglect the effect of the t^ssr, we require that in 
the region in which measurements are used, the formula (29) should be for 
practical purposes identical with the formula ( 2?) • 

If Y is sufficient*^; large compared to Y * the excess mass near the 
center beoomes negligible compared to the total mess within the shock radius 
and the solution must approach asymptotically the similarity solution. This 
is indeed the case since f tends to 1 for large Y. 

For the similarity solution we have the exact formula valid for any 
value of Jf • It can be written in the form 


U 2 


e ( y'-i) 

B ( y } Y 3 


(39) 


where B ) differs only little from the value 2ff/3. B as function of 
y is tabulated in Chapter 5l It is identical with the quantity given !n 
the last line of Table E.3. 

It seems reasonable to fix up the formula (26) in*euoh a way that it 
gives the correct asymptotic behavior; i«e» we put 



5 ( JMJ. -2 

3 (JT > J°o Y 3 


( 31 ) 


In particular, if the formula is used in the region where f approaches 1, 
the effect of a finite on f will be a small order effect and then 

it- is perfectly justified to include the finite ( i -l) in the main term but 
to neglect it in f. 
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6.5 APPLICATION TO THE TRINITY TEST 

The correction factors f have been evaluated for the conditions in 
the Trinity test. Four functions f have been calculated as follows: ' 

In the first instance only the material in the gadget itself was in- 
cluded. For the second curve the X- unit was included* This unit does 
not have spherical symmetry. However, the difference between the two carves 
was found to be small. 

Next, the platform on w/iich the gadget was located was taken into account, 
and finally the tower* Again these two curves did not differ appreciably 
from each other. 

• A 

However there was quite an Appreciable difference between the curves 
which did or did not include the platform. Since the platform has no spher- 
ical symmetry, the difference between the two ourves represents an unavoid- 
able experimental error* Clearly for future teats it would be desirable to 
eliminate this source error* 

Values of the shock velocity on an arbitrary scale calculated with the 
help of these four correction factors are 3hown In Figure 1* For comparison 
a curve which neglects the variable density is also shown. The numerical. 

data for these curves are: 

1 

Gadget only: Y 0 * 8.56 m, Y * 0*36 m 

Gadget and X-unit: Y Q - 3.75 m, Y * 0*37 m 

Gadget, X-unit and platform: Y 0 * 11.9 m, Y ■ 1.0 m 

Tower: Y' * 2 m, Y^ * 5*3 m* 

It will be noticed that the correction factor f is below 1 for snail 
shock radii and above 1 for large radii. The explanation is as follows: 
When the shock passes from the gadget into air, the shock velocity increases 
rapidly- The gadget material, which wa3 shocked by a slow shock, must now 
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Figure 1 

Shock Velocity as Function of Radius 
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be speeded up and this requires that the pressure gradient behind the shock 

ba reversed* Thus the average pressure, instead of being one half the shock 

% 

pressure , (see eouations 5-16) , will now be higher than the shock pressure. 
Eventuollyjas the «ir shock slows down, the gadget material must again be 
decelerated. Since there is now more material, the pressure gradient re- 
quired to do so must be larger end the average pressure drops below one half 
the shook pressure. Hence the equation (8) may be written in the form 

4" Y 3 2n Y 3 „ /,2 

E * — —z — P average - — -rr~r P*/ f w 2 ' 

3 a'-i 3 

where f < 1 at small distance and f > 1 at large distance. With the help 
of (14), it. is found immediatley that f as defined here is identical with 
the correction factor f» 

In the early stages the heavy material at the center ha s, the ref ore } the 
tendency of holding the shock back, until it is accelerated sufficiently to - 
follow the shock. In this stage the assumptions on which the theory is based 
are not very well satisfied* Combining (3l) and (27) on* finds for Y Y 0 


3F. ( Jf -l) 


BOP A V 


in (4-)f for Y« Y 
Y 


and the shock velocity varies little with the shook radius. However, the 


essential reason for the accumulation of matter near the shock front lies 


in the rapid decrease of shock velocity with radius, which is accompanied 
by a corresponding decrease in entropy. Hence, there is normally a region 
of high entropy and low density in the inside of the sphere. This is no 
longer correct if the shook velocity remains constant over an appreciable 

distance. For this reason the theory is not reliable until the second phase 

m 


is reached, and the shook slows down again. In this phase the heavy material 
oontlnues to press outward and , therefore, rsises the shock velocity above the 
value expected from the similarity theory* The magnitude this effect 
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can be fudged from Figure 1 by the fact that the shock velocity ourve at 
one point touches the similarity curve for 2.8 times increased energy 
release. It is, therefore, easily possible that the energy release be over- 
estimated by a factor of this order, if the effect of the gadget on the 
propagation of the shock is neglected. 

From the shock velocity we obtain the shock radius as function of time 
% means of one integration. Using the factor (27) we find 

l/2 


J* 


« • < ^ TJ^ 


'I 


-n . 


Y 

T 


(34) 


F (t, e 0 ^ 


I 



3 

3'fn z/ tc) 


ds 


(36) 


If z is large compared to 
the similarity solution 



1, we obtain asymptotically the formula for 

• ! .. * if <36) 


The function F has been evaluated by numerical integration for two 
values of i 0 , obtained if either the effect of the gadget only is incl'jded 
or if the platform is also taken into account. They are shown in Figure 2 

c;/o 

where Y'"' is plotted against t for an energy release of 21,000 tons 
of TUT. l/( $ -l) was assumed to be 4.3. The value of 21,000 tons was 
chosen in order that the experimental points for shock radii between 20 
and 30 meters should lie between the two theoretical curves. In this re- 
gion the experimental points lie on the dotted straight line shown in 
Figure 2. 

The exper imental data were obtained by Mack's group from photographs 
of the ball of fire by means of Fast*x cameras (see Chapter 18). 

The shock pressure is obtained from (14) as soon as the shook velocity 


APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 


VI - 14 


Figure 2 

*5/2 

as Functions of t. 



APPROVED FOR PUBLIC' RELEASE 









APPROVED FOR PUBLIC RELEASE 


VI - 15 


is Tcr.own* It varies between 15,000 - 26,000 atmospheres at Y • 20 inters 

to 4 , 000 . - 5,000 atmospheres at 30 meters and is } therefore., well within the 
range postulated in the introduction. 
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CHAPTER 7 

THE IBM 3QUJTI0M CF THE BUST WAVE PROBIiM 

W. Fuchs 

7.1 INTRCDUCTICK 

The discussion in the preceding chapters shoes that the problem of 
the propagation of the. blast wave from a nuclear explosion is quite complicated. 
Even if we disregard any transport of radiation, appreciable complications arise in 
view of the large variations in temperature and entropy. We mav roughly divide 
the range of shock temperatures into two regions. Firstly,the region fromabout 
one million degrees to about three thousand degrees absolute. Here dissoclattlni^ 
of molecules and ionization of the atoms takes place. Consequently ^ | is 

fairly small but varies with temperature. a elow 3000° ^ is less variable and 
approaches eventually the value 1.4 in nonnal conditions. 

The temperature varies of course also along an adiabatic. However# the * 
total variation between the shock pressure and one atmosphere is not excessive - 
about a fhctor 2 at shock temperature of 3000° and slightly more than a factor 
10 for a shock temperature of 1,000,000°. The effect of decreasing temperature 
along the adiabatic on the degree of ionization or dissociation is partly 
balanced by the decrease in density. For this reason the variation of ^ along 
an adiabatic is not as pronounced as one might expect from the change in 
temperature. For this reason qualitative statements made about the conditions at 
the shock front hold to a large degree also for the subsequent expansion behind 
the shock. 

A temperature of 3,000° K is reached in the shock when the shock pressure 
is about 80 atmospheres. (For an energy release of 10,000 tons of TNT the 
shockmdius is then 80 meters). We are however more interested in the pressure 
region from about one atmosphere down (corresponding to shock radii of 500 

meters and more). It was felt that the exact energy distribution at this early 
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ftauM'not have an appreciable effect on the pressure distribution in 

stages as long as the energy distribution is at least roughly correct, 
ioyt! trills reason it was decided to start an IBM calculation at this point and to 
£ saleolate the initial conditions for this instant by means of the approximations 
|'i|^|iopad in the preceding chapters. The approximations made include* 

I see section 2) 

' (1) in approximate treatment of the isothermal sphere. In 

Mi- ' ■ ' 

aotual fact the isothermal sphere at this late stage has no significant influence 

on the propagation of the shock. However, in the first instance we were 

initxrfested in the isothermal sphere as such# in the second instance the 
I.tj . 

| i«o thermal sphere greatly facilitates numerical computations, since it reduces 

.[total range of temperatures and entropies which exist at any given moment, 

eliminates the singularity ^at the centre. 

‘ I , (2) Y was assumed to be constant and an average value of 1.25 was 

'J r;j4 p *• • • . ’ 

aipfumd. This is probably the least satisfactory of th* assumptions made# but 


ft.##* 


I vM 


i ‘ ;i 



of the fact that we did not require a very accurate estimate of the initial 
, and that a calculation with variable $$■ Is hardly feasible without 


1 ' ’ IBM machines or a great amount of computation, this assumption seemed justified. 

(3) tf’-l was assumed small. This assumption is not essential, but 
error introduce! thereby Is small and it has the advantage that the 
isothermal sphere can be included as an Integral part of the calculation. 

for the IBM run it is of great advantage if the variation of the pressure 
|ft«n|g an adiabatic is a simple function of the density, though the variation from 
M Xdiabatlc to another may be given in numerical form. The reason is that in 
<Ef|iiiif case" the adiabatic of each mass point is given by one or two constants and 

»i* :: *?', iiif'V; 1 : 

;s|p require only a table of these constants as functions of one variable. 

||itfObleiBS in which the equation of state is given completely in numerical form 
! have not yet been tried on the IBM machines). 
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It would be difficult to find a staple presentation of the equation 

l 

of state covering this whole region of adiabatios between the shock pressure 
curve and normal pressure. However, with the limitation to shook pressure 
below 80 atmospheres the demand of a simple representation becomes feasible. 

In particular thi3 is true for the adiabatics which start at shock pressures 
below 80 atmospheres, since then no ionisation or dissociation occurs. Ve 
require of course also the adiabatics of the inner mass points, which were 
shocked by stronger shocks. However we require only the tall end of these 
adiabatics. Fur the more, the highest entropies are eliminated by the equalisation 
of entropy inside the isothermal sphere. The equation of state which has been 
used is discussed in Chapter 3. 

The requirement of a simple equation of state is the principal reason 
for starting the IBM x*un at such a comparatively late stage. 

Although radiation transport has been taken into accountMnsofar as 

it is responsible for the formation of the isothermal sphere, no allowance has been 

made for the radiation transport from the isothermal sphere into the region in which 

NOg is formed. As explained in Chapter 4, Section 4, this transport of energy 

becomes important when the shock radius has reached about 100 meters and it 

should affect the propagation of the shock shortly thereafter. ^ie opacity data 

required for the purpose of calculating this transport are not sufficiently well 

known. 3h neglecting the radiation transport altogether we are pessimistic, since 

it is of advantage to have the energy close to the shockfront. Then the 3hocj£ 

pressure decreases less rapidly than it would otherwise. The increased shock 

pressure would naturally lead to a greater degree of dissipation of energy by the 

shock so that at larger distances tho shock pressure might drop again more rapidly 

and at sufficiently large distances the effect of the radiation transport on the 

shock-pressure would be reversed. At present we are not in a position to make any 
definite statement about this possibility. ! 
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The IBM run was intended for an energy release of 10,00b tons of TNT. 

Coring to son*? unfortunate circumstances related in Section 3, no definite energy 
can be attributed to the run throughout its whole history. At sufficiently large 
distances the energy should be assumed to be 13,000 tons. Other energies can, 
of course,be obtained by the usual scaling laws. 

7.2 THE INITIAL CONDITIONS CP THE IBM RUN 
The initial conditions of the IBM run were prepared by Hirschfelder and 
Ibgee. The principal data are sunnarized below without going ijito the details 
of the calculation. All data are for an energy release of 10,000 tons of TNT. 
7,2-1 The Isothermal Sphere 

The Iagrange radius r© of the isothermal sphere is given in terms of the 
shock radius Y by Rjuation 4 of Chapter 4, flection 3. It can be written in the 

form „ . 

r -[ 11.85 Y* m6 -601. 5l X / 3 100 cm 

° L J (1) 

Here both r Q and Y are given in centimeters. It was convenient to chopse a 
simple value of Y/r Q and the value 4 was chosen, which correspond® to a shock 
pressure very near to 80 atmospheres. Then 

r o = 1997 cm, Y > 7987 cm, Y/r 0 = 4 ^ 

The actual radius R Q of the isothermal sphere is obtained from the conservation 
of mass. Since we assume constant density and constant pressure in the 
isothermal sphere one has 

r o 3 (3) 

r as function of Y is shown in Figure 1 of Chapter 4. It will be seen that r 

o 

varies very slowly after the shock radius has reached about 10 meters and the 
effect of the isothermal sphere on the shock is negligible a short while 
thereafter. B e yond a shock radius of 80 meters r 0 varies very slowly indeed. 
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Temperature 

Pressure 

Density 
Entropy AS/k 


According to Figure 2 of Chapter 4, it increases from 20 meters to 2? meters as 
the pressure in the sphere decreases to 1 atmosphere ( the analysis is no loiter 
valid for such small pressures, but it indicates the order of magnitude). 

For this reason it is a good assumption to assume that r 0 is constant. 

The actual radius R 0 then varies in accordance with Equation 3 only because the 
air in the isothermal sphere expands. initial value of R 0 is 60 meters (see 

Figure 1 Chapter 4). 

The initial condition of the isothermal sphere is given by the following 
quantities t 

» 49,000°* 

« 37,0 atmospheres ','5 

s 0,0392 x normal density 

= 85 rl ' 

Internal energy and enthalpy E/R = 1.487 x 10^, H/ft a 1.782 x 10^ 

• * 

R 0 a 60,23 meters, r 0 a 19,97 meters. 

These data were obtained from a calculation indicated below, 

7.2-2 Initial Pressure and Density Distribution 

It his been shown in Chapter 5, Section 5 that the isothermal sphere can be 
treated on the assumption of small |f— I. The small }f — f approximation has been - 
checked far a point source solution (see Chapter 5, Section# 3 and 4) and it was 
found satisfactory in the region in (Which we are interested. Since the solution 
which we require is in any case very close to the point source solutions except in 
the neighborhood of the isothermal sphere, the error of the saallCjC-f) approximation 
is soil. 

Wo snail not go into the detail# of the calculation. The analysis is ratter 
involved, but the lines along which it proceeds are sufficiently indicated in 
Section 6 of Chapter 5. The resulting equations had previously been evaluated 
for two values of K in order to see how sensitive they were to a change in Y . 
The values chosen were Y - 1.2 and Y ~ 1. 3, r 
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The values for ^ = 1.25 were then obtained by interpolation. In this way 
we find the initial pressure and density distribution as well as the velocities 
and the Euler coordinates of all mass points. 

Sons adjustments had to be made on the initial conditions. A minor 
adjustment arose from the fact that the small IS — I treatment of the isothermal 
sphere does not agree exactly with the treatment given in Chapter * Section 5. 

The error, which may be due to either method, may be seen from the values of 

R 0 /T* the small — 1 treatment gives R 0 A = 0.850 compared to 0,761 by means of the 

other method. The latter value was assumed to be more reliable. 

Furthermore , at the start of the IBM calculation the value of y at the 
shockfront is larger than 1.25# instead of a compression ratio in the 3hock 
A /A a 9,as would be expected for y = 1.25, the value obtained from the correct 
Hugoniot curve for a pressure of 77.25 atmospheres is J~l/ = 7.24. The 
density contour was, therefore, adjusted to give the correct compression ratio at the 
shock, and the correct radius of the isothermal sphere. This required also an 
adjustment in the Euler coordinates R, since it is essential that the initial 
conditions satisfy the equation of continuity* 


A 

/> 


ft 

~jr 


Aft 

dr 


(4) 


where r is the Lagrange coordinate. 

The initial velocities were then calculated directly from the equation 


u=Y (/ 

(UiV 


. 3(y-o// 

• + c * /Y) 


(5) 


which follows from the small # — I approximation. Here the correct shock 
velocity V for the given Shock pressure was used. 
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7.3 THB TOTAL B NERGY 

mm Mina 

The methods used to establish suitable Initial conditions are in parts 
somewhat arbitrary, For this reason it is not surprising that the total energy 
corresponding to these conditions turned out to deviate appreciably from the 
assumed value of 10,000 tons of TNT, Unfortunately the energy was recalculated from 
the initial conditions only after the IBM-run had been completed. It was then found 
that the total energy was 13,500 tons TNT, j 

Since the initial shock pressure of 77,25 atmospheres at the initial shock 
radius of 79.9 teeters corresponds to an energy release of 10,000 tons, we have no 
remedy for the discrepancy. All that can be said, is that the shock pressure versus 
distance curve corresponds to 10,000 tons up to 80 meters shock radius and to 13,500 
tons for large radii, For intermediate radii it should slowly change between these 
values. * 

Actually the discrepancy is slightly less. A check of the total energy 
at a shock radius of 2000 meters gave only 13, 100 tons. The "loss* of 400 tons is 
entirely due to errors of the IBM-run and is of the order of magnitude to be 
expected from this source. 

For most purposes the total energy in the IBH-run should be assumed to be 
about 13,000 tons, except at small shock radii, where 10,000 tons is more appropriate 

t V 

7.4 THE IBK-RUN £ 

The hydrodynamical equations are 



In addition we have the equation of the adtabatica, which were put in the form 
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where g and h are numerically known functions of the entropy. 

The boundary conditions are that the relocity dR/d t vanishes at the 

center and that at the shock^adius the Hugoniot conditions are satisfied. The 

\ 

latter determine also the entropy of any mass point as it passes through the shock, 
front. The entropy is assumed to remain constant, so that essentially g and h are 
given functions of the Iagrange variable r 6 

The methods employed in solving the system of partial differential equations 
by means of the IBM machines are explained in Volume 2 of this series. The 
method used first was suitable as long as the shock pressure differed appreciably 
from one atmosphere, but it became erratic as the overpressure became small. The 
method was therefore changed, so as to calculate changes in density and pressure 

’ 'm 

rather than their absolute values. Ails change of procedure quickly suppressed 
the erratic behavior of the pressure. 

The run was continued until the shock radius had reached a tsalue of 6^270 
meters. At that instant the overpressure in the shock waa 0,0251 atmospheres. The 
positive pulse was 290 meters long and the negative pulse 760 meters. Since the 
further propagation of the shock is influenced at these low overpressures only 
by the positive pulse, the approximations on which the semi-acoustic theory of the 
next chapterare based, are well satisfied. They are (1) that the overpressure be 
small compared ko 1 atmosphere, and (2) that the length of the pressure pulse be 

« 

small compared to the shock radius. Even the application of the semi-acoustic 
theory to the negative phase is not bad. Hence, the IBM- run was discontinued and 
the semi-acoustic theory was used for the purpose of continuation. 

v 

7,5 RESULTS 

The shock pressure as a function of the distance of the shock f resit from the 
center of the explosion is shown in Figure 1. 

In this graph all data have been collected from the various chapters. From 
a shock radius of 10 to 80 maters the similarity solution has been used. The 
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Figure 1 

Shook pressure versus distance for 
nuclear explosion. Tins is given 
In seconds, t * 0 is 0,012 after 
start of explosion. 
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dotted lines show upper and lower limits for t-he effect of the gadget material on 
the propagation of the shock, as calculated for the Trinity gadget (sealed down 
to 10,000 tons of TOT), In this region the curve corresponds to an energy release 
of 10,000 tens. 

For shock radii from 80 to 6300 meters, the shock pressures are obtained from 

the IBM run. Here the total energy is between 10,000 and 13,000 tons, the upper 

. »\ 

value being correct at sufficiently large distances. The time of arival of the 
shock is indicated at various points, t * 0 is the start of the IBM run which was 
0,012 second after the explosion. 

Beyond a radius of 6s3Q 0 meters up to 67,000 meters, the semi-acoustic 
theory of the next chapter has been used. 

A number of curves showing the pressure at a fixed distance as function of 
time are shown in Figures 2 to 9, Graphs of the pressure versus distance at a 
fixed time are shown in Figures 10 to 18, 

The duration of the positive phase of the pulse as function of the shock 
pressure is given in Figure 19 0 

For comparison with measurements at Trinity we have also made a graph of 
the arrival time of the shock at various distances. It is most convenient to plot 
the shock radius divided by the arrival time as function of the shock pressure 
divided by the normal pressure, since such a graph is independent of the energy 
release. This graph is included in Chapter 19 * It depends, of course, 
on the norm! velocity of sound in the given circumstances., c 0 * For the IBM- run 
c 0 = 347 centimeters! per second. 

Finally, the re are shown in Figure fcO the positive impulse 1+ and the 
fraction of the total energy which is left in the blast as functions of the shock 
pressure divided by the normal pressure# Tfoe latter is independent of the 
energy release. The positive impulse has beep soaled to an energy release of 
40,000 tons in free air (or 20,000 tons on the ground) for the purpose of 
comparison with observations at Trinity, 
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Figure 2 

* 

Pressure versus Tine for Radius 170.7 sMters 
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Wgur* 3 

Pressure versus Time for Rsdius 223.4 eaters 
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Figure 4 

Pressure versus Tine for Radius 332 inters 
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Figure 5 

Pressure versus Time for Radius 466*5 jesters 
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Figure b 

Pressure versus Ttws for Radius 682 msters 
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Figure ? 

Pressure versus Time for Radius 1224 mstsrs 
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Figure 8 

Pressure versus Tijoe for Radius 2031 meters 
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Figure 9 

Pressure versus Time for Radius 3586 meters 
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Figur* 10 
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Figure 13 
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Figure 14 
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Figure 15 
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Figure 16 
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Figure 17 
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Figur* 1 B 
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Figure 19 

Deviation of the positive phase of a pulse at a 

fixed distance 
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Figure 20 

The positive impulse and the energy in the bleat. 
(IBM) run scaled to 40,000 tens. 
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7,6 COMPARISON WITH TNT EXPLOSION. EFFICIENCY OF NUCLEAR BOMB 
One purpose of the IBM-run was t o find the ef fi ciency of a nuclear bomb 
compared to an explosion from an equivalent charge of TNT. In the nuclear explosion 
a greater amount of energy is used for the purpose of heating the air near the center 
of the explosion to high temperatures. A large fraction of this energy is useless 
for the propagation of the shock. 

Since no comparable IBU-run exists for a TNT explosion, we compared the 

v 

results for the nuclear explosion with experimental data. For this purpose, the 
experimental curve prepared by Hirschf elder. Lit tier and Sheard was used. It is 
based on experimental data for charges fired on the ground, and for shock pressures 
in the range from 15 to 2 pounds per square inch. The charges varied from 67 
to 55Q pounds. The curve is given by the expression 


A P - 


38,6 


85 


5760 


= Y/w 




(9) 

( 10 ) ^ 


Here I is the shock radius in feet* w,the weight of the charge in pounds, and 
the overpressure in psi. 

For lower pressures Hirschf elder, Sheard and Littler used the asymptotic 


formula 


35,0 


Ap - x yj i og x - 0.928 


( 11 ) 


The constants were obtained by fitting to the low pressure end of the curve (9). 

The analytical farm of the equation follows from the semi-acoustic theory presented 
in the next chapter, provided the pressure pulse has reached itaaeyraptotiV linear 


shape. 

From these formulae we calculated the weight of the TNT charge required to 
give the same shock pressure at the same distance as the nuclear explosion. 
Dividing this charge by the weight of 13, COO tons assumed for the IBM- run, and 
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multiplying by a factor 2, to take care of the fact that the TNT data apply to an 
explosion on the ground so that practically all energy is in a half sphere, we obtain 
the efficiency of the nuclear bomb compared to an equivalent charge of TNT, 

The efficiency defined in this way depends on the shock pressure chosen for 
comnariserv, and it 1 a not surprising that it should vary with the shock pressure* 
no w&ver'j it is surprising that the efficiency should increase with decreasing shock 
pressure as shown in Figure 21 over the range in which experimental data for TNT 
exist. However, the individual experimental points scatter appreciably, especially 
for high shock pressures and 7 theref ore, thi3 variation of efficiency with shock 
pressure is not necessarily correct. The most reliable data are those at the lower 
end of the experimental range, which would indicate an efficiency cf about 0,6, 

The curve has been extended on either aide to higher and lower shock pressures. 
The extension to higher shock oressurea by means of formula (9) is quite arbitrary 
and has been performed only since experimental data for the Trinity test have 
usually been compared with this formula. The extension to lower pressures depends 
on the assumption that the asymptotic formula is valid in this range, which is not 
necessarily correct. Again, the main reason for performing the extension is 
the fact that experimental 'data have been analyzed oy means of this curve. 

A similar comparison has been made by using the experimental curve of 
A. H.Taub,^ for half-pound TNT charges in free air. These differ quite 

rn 

Report NDRC-A-4076 

appreciably from the curve of Hirschfelder, Lit tier and Sheard, as may be seen from 
Graph 21. This discrepancy shows the difficulty of assigning any definite value 
to the efficiency of the nuclear bomb for the purpose of producing blast and makes 
It desirable to have a comparable IBM- run for TNT. 
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Figure 21 

Efficiency pf nuclear bomb as compared to an equivalent 

charge of T.N.T. 
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It should be noted that two factors, which might be important, hare been 
disregarded in the IBM- calculations. One of these is the radiation transport 
Of energy frost the bail of fire into the N©2 layer, discussed in Chapter 4. 

This has the effect of equalizing the temperatures in the interior of the shock 
sphere and it will make the nuclear explosion more similar to a TUT explosion. 
Consequently, the efficiency of the naclear bomb is Improved. The other effect 
is the radiation which escapes to l^trge distances. The arguments presented in 
Chapter 4 show that this radiation occurs at a time and place whsre it cannot 
affect the shock pressure in the region of practical interest, i.a. down to 
about one psi. However, it could reduce the efficiency at very low shock pressures 


7,7 SCALING LAWS 


% - u 


The IBM -run was made assuming a normal density of air - 1.163 and a 
normal prsssure P 0 » 1 bar (10^ dyn/cm 2 ). The sound velocity ofnonaal air is 
then c 0 s 347 centimeters per second. The total energy was Q x 13,000 tons. 

For the change to a different energy release we have exact sealing laws, 
we keep all pressures and velocities fixed, but change all radii and times in 
the ratio of the cube root of tbp energy release. 

For the change to different normal densities or pressures no exact scaling 
laws exist. However, within the accuracy of the primary data which enter the 
IBM-run, we can use the following scaling laws, which relate the dashed quantities 
for arbitrary normal conditions to the IBM- quantities t 


P 1 

3 JEL 

S' : Xfi 

Po 

Po 


1 


t » a JSjL X t 

u 1 

a U 


c o’ 

c o 

Q» a Pi q 


a 

Po 


( 15 ) 
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i.e. we scale all pressures, velocities and densities in the ratio of their 

normal values. The linear dimensions are scaled by an arbitrary factor A , the 

times by the same factor A and the inverse ratio of the normal sound velocities, 

s 3 

and, finally, the energy release is scaled by A and the ratio of the normal 

pressures. The scaling laws for pressure, velocity and density are, of course, not 

independent of each other, sin«^ # C^ C o V ( p °/j Po | 

These scaling laws are based on the assumption that the equation of 3 tate 

of air can be written in the form 

— - funation of ( 

0 o Ho 

which is true for a 2f-law, but only approximately correct for the true 
equation of state of air . (E is the interval energy per unit mass). 

For the Trinity test the scaling changes are appreciable, because of the 
low value cf the normal pressure (See Chapter 19), 
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CHAPTER 8 

ASYMPTOTIC THEORY FOR SMALL BUST PRESSURE 

H. A. Bethe K. Fuchs 

8.1 INTRODUCTION 

It is very desirable to have a theory of shock waves which is valid 
for small over-pressure . One purpose of such a theory is to provide a 
natural transition to the well-known acoustic theory, and to set in evidence 
the limitations of the latter. Another purpose is to give practical results 
for the pressure up to arbitrarily large distances and to make it possible 
to stop the numerical calculation with IBM machines (see Chapter 7) at some 
finite low pressure. This numerical calculation would become increasingly 
inaccurate and cumbersome with increasing radius of the shock wave; the use 
of asymptotic formulae therefore improves both the accuracy and the ease of 
the calculation. 

The theory presented in this chapter represents the first terms of an 
expansion in powers of the ratio I/* where Y is the radius of the shock wave ISP 
and L its length, i.e. the distance from the first shock to a point at which 
the pressure has decreased to a small fraction of the peak pressure . We shall IE 1 
in general retain the first and second pcwer of L/Y and neglect the third. Mil 

One strong reason for stopping at just this .point is that the entropy in the 
shock is proportional to p^, where p is the over-pressure at the shock front 
which in turn is proportional to l/Y. Therefore, if we neglect terms of order 
(LA ) 3 we can consider the entire process as adiabatic which involves a great 
simplification. At on© point in our development, however ,( Section 7) we shall 
calculate tfri© actual energy transformed into heat at the shock front and at 


tiii - i 
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this point wo shall,, of . course, carry terns of order (l/Yp. 

Ordinary acoustic theory does not represent the first term in an expansion 
of the type considered here. We shall find that the length of a spherical 
shock wave increases gradually as the shock wave moves away from its center, 
whereas an acoustic wave retains its wave-length all the time. Similarly, the 
decay of the pressure at the shock front is somewhat faster than that of the 
pressure in an acoustic wave. However, acoustic theory is very useful to 
provide some guiding ideas for our theory, and we shall therefore start by a 
recapitulation of the acoustic theory of an outgoing spherical wave. 

6.2 ACOUSTIC THEORY 


In the ordinary hydrodynamic equations (see below. Equation (12)) all terms 
should be neglected which are of second or higher order in pressure p or 
material velocity u. Then the hydrodynamic equations become 

s - c 


du _ _ „ 2Lflt 


Tt~ 




*fr- -•(& 


( 1 ) 

where <j" is Riemarm's quantity, (cf. Equation (10))which for small pressures 
is given by 


r r t/fi % 


Combining the two equations (1), we get the wave equation 



fr* v 32, <C ■ 

Hi r 2 




( 2 ) 


which has the solution 
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In which e Is the sound velocity which itself depends on the density* 
In (lo) we have used the definition 

c2 s dp/d J* 


( 11 ) 


Inserting (lo) into (9) we obtain 


b u 

yt 


+ U 


bu 

Sr 


o 


fttr 
b r 




( 12 ) 


In the acoustic limit the terms containing u as factor can be neglected 
and (12) reduces to Equation (l)« 

We know from the acoustic theory that <r and u behave approximately 
as 1/r at large distances. We therefore introduce the abbreviations 


2 = 

U -- 

’With this notation Equations (12) 


r <r 

r u 

become 


b U 

af 


+ e 




0 


du 

aT ♦ c 57- 


+ 


u 


fciL + c 

1) r 


u 

r 



0 


.(13) 


(14) 


We know further from acoustic theory that u becomes asymptotically equal 
to T for large distances. We therefore further set 

U 2 • D < l6) 
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and we expect tnat D will become small compered wi th i for large r. 
Equations (14) now become 

> (c»4 u) If. +c (i£ + 4-) - T# = 0 

/ 

- (c ~ u) - c - o JBj 3 0 

In the first approximation we nay neglect the last three terms in the 
first equation (16). We should expect that d/dr of ary quantity is of the 
Order of that quantity divided by L. Therefore the term U/r, which is 
approximately equal to x/r , is small co^jared dz /dr, the ratio of the 
two terms being of the order L/r. The term D/ r is small because D 
is expected to. be small compared to j as will be proved below. The last 
term is in turn small compared to the second last term,. 

In the first approximation, therefore, the first equation (16) reduces to 

+ (c +«) = 0 r ; , (17) 

• . . , : ■ <• 

This equation means that Z propagates With the velocity c - u. This is 

1 % 7 :\ '***" 

the analogue of acoustic theory in which r propagates with velocity c 0 * 

The replacement of the ordinary sound velocity c 0 by the effective sound 

" 7 'i 1 

velocity c - u is physically obvious and ; is analogous to the well-known , 

Rieman method for treating plane problems. 

The eolutioa of (17) is 

z * zr ^ ( is ) 

with 

t * t - f -..dr 

J «> + u (19) 
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The integration in the expression (19) is to be extended over e path 
along which L is constant. The integral therefore depends on l and 
the equation does not give an explicit solution of the Equation (1?). The 
lower limit in the integral may be e function of £. If we choose the 
same lower limit r o for all £, then 't is the time when the value of N 
£ to which it belongs passes the point r Q . For the present we shall 
not specify the lower limit. The function £ (^tO must be determined 
from the shape of the shock wave at some initial time. This will be dis- 
cussed in Section 8*10. 

8.4 SECOND APPROXIMATION 

We shell now try to determine the function D from the second Bqua- 

N 

tion (16) and also to get a better approximation for I. 

In the seoond Equation (16) we can certainly neglect the last term. 
The first taro terns can bo transformed by introducing instead of t the 
variable / £ defined in (19). We have then 

( d“t“)r * ( 

( 6 D ^ „ fh p \ 1 

U r Jt ^ v ^ /r 0 (20) 

and the second Equation (16) becomes 

rh (|Vi- 

( 2 1 ) 

» 

* 

Neglecting all terms of smaller order of itegnitude } tlie right hand side re- 


\ 
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duces to 2 c i/r 
grated end gives 

D 


end the factor on the left to 2. 
r *r ^ 2 dT 


Then .(2l) che be inte- 


( 22 ) 


This result shows that D is actually of the order 2 l/r, end is there- 
fore small compered to 2 according to our assumptions. This justifies the 
first approximation given in the last section, as well as the treatment given 
in this section* It is, however, of interest to write down the right hand 
side of Equation (2i) to the next approximation. For this purpose we note 
that according to (22) 



T 



(23) 


Therefore (2l) becomes 

( JlIL) - C o 4 u v '^ 

in which only! terms of order uD/r are negleoted, but terms cf order cD fr 
and u l/r have been taken into account. 

The form of Equation (22) shows that D has the character of an after- 
flow term and therefore corresponds in all respects to the after-flow term 
in Section 2. 

The integral in (22) must be extended from the shock front to the point 
at which D is to be calculated. To prove this, it must be shown th*t D 
behind the shock front is of smaller order then (22). This can easily be 
done by using the Hugoniot relations and the definition of CT • We hove 
behind the shook front 


2 

u • 


(P-Po) (V 0 - V) 


( 25 ) 
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where the subscript 0 refers to the undisturbed material, and V is the 
specific volume* Further, we have generally 


<r 


-fl 


J> 

•> 


dj» 


4^ 




d V 


(26) 


The value of D behind the shock front is by definition 


D - r ( u - cr) 


(27) 


To evaluate (25) and (26) it is convenient to introduce the quantity 


x - V 0 - V 

TT* 


(28) 


and to express the pressure in terms of this quantity: 

P ■ P 0 - V x *4* l£_ £..2 r 2 - 

O o dTV d V 2 nr d y3 


(29) 


Using the definition of the sound velocity we have 


-u\-- 

*0 

and we write 

C„ 2 

v 0 2 

i 

(30) 

53l|, 

2<* c 0 2 



Tv® 

'6 

= T* 

' 77 i. J v ! 

v 0 



o 

- " / »■ \ A* 

Vo ; 

(31) 


where oC and are dimensionless coefficients. If the adiabatic law is 
of the usual fora. 


P * Po ( V/ V G ) 

the values of and /3 are 
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<* = c * + 1 ) A 

a s C x + 0 ( 7 + 2 ) 


(33) 


With these abbreviations Equation (29) becomes 

c *■ 

F ’ P 0 = ~ ( *+<* x 2 +(9 x 3 *..) 


(34) 


Inserting (34) and (28) into the Hugoniot relation (26) we obtain 

u * - «o x (14-£c*x + (£^- *£<* 2 ) X Z •+■••*) ( 35 ^ 

Further we find. 


~ T"v " "Y" ( l4c **4 (-»r- ** 


(36) 


or for the sound velocity 
c s e„ 


[l4 (c*- 1 ) Jc 4 ( z ft/2 - o?/2 -o<) x 2 ^ . .. 

•wa 


(3 7) 


Inserting (36) into the expression (26) we find 

(T : © 0 x £l 4 (o< /2) x 4 ( (3/2 -c* 2 /6 ) x 2 + • • ] 


(38) 


and subtracting we get 


u - <r* sc 0 <*, 2 x 3 /24 


( 39 ) 
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This relation holds generally behind the shook front. It is interesting 

that it depends on the quantity o( , i.e* on d^p/dV^ which is also the 

quantity determining the entropy change at the shock front (see below. 

Section 8 . 5 ). For our present purpose the important point is that u-<X“ 

3 3 

is proportional to x and therefore to CT » and therefore falls off as 
;i/y 3 where Y ie the shock radius. The quantity D .lust behind the 
shock front will therefqre fall of as l/Y^ which makes it of smeller 
order of magnitude than the expression (22). 

We note in passing that, according to (35) 
of wave propagation is in first approximation 

Ct-U S C fl ( Itat X) 


and (3?), the ▼eloc ity 


(4q) 


We shall now insert our expression for D baok into the first Equa- 
tion (16). We have to calculate the quantity 



Using (2o), (23) and (24) and the definition of 'XT » this becomes 

. 0 . 2_ - _ ® 

r r r 


(42) 


As was explained above, only quantities of the order (u/c) (D/r) are neglected 
in this expression* 

Inserting (42) into the first Equation (le) we see that the terms u/r 
cancel. Therefore the first Equation (16) becomes, negleeting the very 
small term D fc/r^ s v 



+ (c 4 u) 

it 




(43) 


Neglecting u 


in 


c * u, this can be integrated to give 

y ~~ V 
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In other words, along a line of oonatant 't , not £ itself is a constant 
but the somewhat more complicated quantity 

2 _ y 

(44a) 

£ itself therefore increases as the ware propagates outwards. 


8.5 THE VOT 101? OF THE SHOCK jKOKT 

'i " 

Let us assume that the shock front moves from a radius 7 to a radius 
7 + d7« The shock velocity can easily be shown to be in first approximation 
midway between the effective sound velocities ahead and behind the shock. 

The former sound velocity is c Q , the latter c ♦ u so that the shook veloc- 
ity is 

V = i (\ + c+ '0 

(45) 

v 

and the time needed for the shock to travel the distanoe d7 is 


ctt a 


xdL Y 


(46) 


The shook velocity is smaller than the velocity of sound waves behind 
the shook, c ♦ u« Therefore these sound waves will catch up with the shook 
wave, and if the shock is followed by a rarefaction as it is in the case of 
a blast wave, then the rarefaotion will gradually cut down the strength of 
the shock. In our notation, the value of X at the shock will gradually 
change $ if it has the value t when the shock is at 7 it will have the 
value T 4 cl “t when the shook has moved to 7 «- d7» 

In order to calculate the variation of along the shock front, we 
have to complete the definition of which was left somewhat arbitrary. 

It is convenient to define the arbitrary constant of integration in (19) 
as follows 
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^ = t - 



where Y is the shock radi ias at* t he t ime wH® n the 8 ns 1 1 * re a c he & t he 

shook front# Y therefore , is a function of i. and this is permissible 
since the constant of integration was allowed to depend on i*. If we use 
the definition above then^C + ^/c o is the time when the signal I. reaches 
the shock* or t is the difference in the actual tin* and the time which 
would be required for a signal traveling from the origin with normal sound 
velocity c w - 

The time difference dt between the arrival of the signals x end 
X dX at the shock front is therefore 


dt = <At 4- ^ 


(4?) 


Comparison with (46) yields 


d/t _ 2 , Cv-C-iA _ C» - C - H 

AY ' '"C ' * 


2 V 


(48) 


This equation describes the catching up of the rarefaction wave with the 
shock front# 

From (4 q) and (38) follows 

C. * U. s *♦* O< 0 - 


(49) 


Hence (48) becomes 


dt 

<£Y 


- oCcr 
xc* 


oL 2. 

xc 0 is f 


(50V 


where we have also used the definition of 2 « or Y» 
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8.6 RESULTS FOR VERY LARGE DISTANCES 


Up till now we have not made any special assumptions about the function 
Z(r), i.e. about the shape of the shook wave# We can get further results 
if we assume that initially (at small distances from the origin^ the shock 
wave was of short duration, and tlvat most of the duration-. »t ahnok radius 


Y is due t o . the <j x co 3 s„ of -the- -aback, Telocity, .over the unperturbed. _agugd 
ve locity This assumption is quite well Justified for blast waves at 

sufficiently large distance* If t 0 is the time at which a given v n lue 
of I exists at a given small radius r Q then we have for the arrival time 
of this value of I at the distance r (cf. (4o))^-^ 


»t # + / -.4*. - ^ 

*0 + 

f 

J Mu 

J 

n. 


'V 

/L 

Cl 


0 

% , t - 2 pc . 

" cT ° c 1 

log 

ft 

A 

0 


4v i 

C -+«s( 2 7 v 


I /X.T "ts 

c%> 


(5l) 


If we neglect small difference in. the values of t„ and r 0 for dif- 
ferent values of Z it follows that for fixed r, Z is a linear function 
of t and therefore also the pressure assumes a linear shape* 

We apply the equation (6l) in particular at the shook front r « Y, 

Z a Z g . Then t is the time of arrival of the signal 2 S ft t the shook 
radius Y* However, this time was also equal to\ +Y jc ( where is the 

value a corresponding to Z # . Neglecting again r 0 and t 0 as small 
* 

quantities, we find therefore 

1 J I s «6 y 

X -- - log JL 

* (62) 

0 c 

In the approximation in which r Q and t 0 are neglected 1 T,. I is by 

definition the difference in the arrival t^rae of the shock and of the signal 
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l m 0# since the latter travels with sound velooity o 0 ; i*e*,lT s l is 
identical with the duration of the positive pressure pulse* This is easily 
verified from the Equation (6l). 1 

In order to obtain *t s as function of 2 g ,only, we have to express 
the shook radius Y in terns of 2 g • For this purpose we substitute (62) 
lrtto (5o) • Then we get 

' 

AY tY log CY/%) ( 63> 

which can immediately be integrated to give 

where A is a constant* In other words, the value of 2 ft t the shook front 
is not constant but decrees** du* to the catching up of the rarefaction* 

This decrease is very slow; it goes only as the inverse square root of the 

’w. 

logarithm of the shook radius* 

The derivation given above may give rise to the impression that r e may 
be ohosen arbitrarily as long as it is small compared to all values Y for 
which the formula (54) is applied* This is not so* In fact, it is easily 
shown that we are led to contradictions if we allow r 0 to vary over any 
appreciable range* 

W* shall later that * detailed analysis will allow us to determine 
the correct value of r Q which should be inserted in the equation (54). 

For the present l*t it suffice to point out that a pressure" pulse of the 
asymptotio shape, requires two parameters to determine its initial conditions, 
namely th^ strength and duration of the pulse* We require, therefore, two 
constants in th* asymptotic law and thus both the constant A and the con- 
stant r 0 are determined by the properties of the pulse* 
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Going back to Equation (6l), we fled then that the duration of the 
pressure pulse will increase as the square root of log Y» The product 
of peak pressure and duration ® , or the so-ealled impulse of the 

wave will therefore behave exactly as i/y, namely 

fy & * const./ y (66) 

We have thus seen that in our approximation a shook wave will spread 
out gradually* in oontrast to the acoustic theory in whioh the wave retains 
its shape for all time. The spreading of a shook wave in 5 dimensions is 
very slow. It would be considerably faster in two and still faster in one 
dimension* On the other hand* for a epaoe of more then three dimensions 
this effect will not ooeur* 

Connected with the spreading of the wave there is a decrease of the 
front pressure whioh is faster than i/y» Again, lit three dimensions this 
effect is small* in one d linens ion f elei*entary acoustic theory would give 
constant front pressure whereas the actual behavior is as l//v' • The 
case of two dimensions is again intermediate* and for more than three 
dimens lone, acoustic theory beeomes the correot asymptotic limit. 

Another interesting phenomenon whioh is connected with the variation 
of the sound velocity c ♦ u is the fact that a seoond shook must be formed 
in the negative phase of a shock wave. The region* of negative pressure 
have a particularly small propagation velooity, c ♦ u,which is smaller than 
normal sound velocity e 0 » Therefore the very end of the shook wave tends 
to catch up with the regions of negative pressure and a seoond shock will 
result. This will be dlsousted in more detail in Section 8*9. 

Another interesting problem is the motion of the point 2 ■ 0, which 
marks the sod of the positive pressure phase ( p> ) and the beginning of 
the negative phase* In our approximation the propagation velooity of this 
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point is just the normal sound velocity o Q . Actually, however, p 1* 
somewhat greater than <r by the amount D/r* Therefore the point X ■ 0 
propagates with a velocity slightly greater than normal sound velocity* 
Since the very tail of the shock wave moves with velocity c Q , the end of 
the positive phase I - 0 moves with a velocity slightly greater than the 
end of the negative phase* Consequently the positive phase will tend to 
become somewhat shorter and the negative phase somewhat longer than the 
elementary theory indicates* This corresponds to observations* 

A very important remark ab6ut the shape of shook waves has been made 
by Penney* It la most easily deduced from the fact that the total mass of 
air* behind the shook must be equal to the original mass of air within tha 


radius Y. This means 

.V 


/ 


J> a/ <U = V /* 


(56) 


There are two regions in which the density ie appreciably different from 
the normal density * One is the central region in whieh the gases 
have been left at high temperature by the shook and therefore have low den- 
sity* If the shock wave is far out, these regions have returned to atmos- 
pheric pressure and therefore to a definite density* If X denotes a 
radius small compared to Y , but large compared to the region of the hot 


gases, we shall have 

x 


./ 


J*f i/'-ctA, = jo x /* - M 


e 

t 


( 67 ^ 


where M is a constant independent of X and of the time* Subtracting 
(57) from (66 ) we get 


S 




(58) 


The shook region itself may be assudmd to be small in extension oompared 
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to Y. Therefore we can replaee 5 in this region, r by Y» Furthermore the 
deviation of the density from normal, J* - # is proportional to CT * 

Therefore we obtain from (58) : 

Y* f ctcLju = V C 2a* = M 

shock -ihook (69) 

From this follows that 

/ 1 dtu ~ Zr ( 60 ) 

However, we know that the value of 2 at the front is nearly constant* 
Therefore the integral in (6o) must consist of two contributions which nearly 
cancel each other* The shook wave must consist of a phaee of positive or 
(over-pressure) and a phase of negative <r (under-pressure) such that the 
impulses of the two phases canoel each other in first approximation. This 
argument is also a proof of the existence of the negative jftiase; it was first 
given by Penney using the energy rather than the amount of material* ' 

In particular, in the limit in which the pressure depends linearly on. 
the time, the shape of the shock wave becomes symmetrical, with equal shocks 
at the beginning and at the end as illustrated in Figure 1* 

8.7 THE 3NBHGY 

The energy flux through a given surface consists of two parts, namely 

■fo 

(1) the work dona by the material on one side,, on the material on the other 
side and (2) the energy transported with the material itself* Since we 
have proved at the end of the last section that the pressure pulse in a 
shock wave has ti negative phase balancing the positive phase, the final 
displacement of any point is sero in first approximation 

J~ UL<tt w oO/Y) 
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Figure 1 

Symmetrical, Equal Shook Wave* 
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This want that the net transport of energy with tha material is sero of 
higher order than tha work term mentioned as (l) above* 

Therefore the energy flux In the shock wave through a sphere of radius 


Y ii in sufficient approximation: 

A 


W = HTT Y 


/o 


udt 


(62) 


If we set u - <r and p * J>C QT and <y = W and if we further assume 
that the pressure distribution is as indicated in Figure 1, the energy flux 
becomes 


W-¥TT/>C dt w HTTfcI^e j * X dx 


(63) 


•here & is the duration of the positive phaae which, according to Equation 
(52) , has. the value ** 


6 = 1^6 

e © 


(64) 


so that 


SC . f 

Using the relation (54) for £ this becomes 

A 3 

w tlJoLp - 

3 c Vlcg(Y/ft. ) 


( 66 ) 


( 66 ) 


Equation (66) shows that the total energy in the shock wave decreases 
slowly as the shock w«ve propagates* This fact was first pointed out by 
Penney and means that it is impossible to define in any general way the 
energy wasted in the shock wave, but that this waste will depend on the 
distance to which the shock wavs has gone* It is interesting to determine 

t • 
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the decay of the energy in the ehook wave also in a different way* 

This way is based on the fact that energy is irreversibly converted 
into heat at the shock front at all times* To determine this wasted energy 
or increase in entropy, we u»e the Hugoniot equation for the energy 


e -e„ = t < rr.ycv. -v) 


(67) 


where the quantities without subscript are behind, those with subscript o 
in front of the shook wave* The energy can be expanded in a Taylor series 
in the change of volume and of entropy, as follows 

E - E„ * Cv-'O + i ^ CV-V„) X 




Cm-vJ +■ 




( 68 ) 


If we remember that 






(69) 


and introduce the abbreviation (28) , (66) becomes 
, p £ a P ^ ~ V*** 


+ ir |^v„ s A.. + t.(s-s) + .. 


( 70 ) 
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Further, we he we 


!'v*v/x 


x x 


Iaeerting (7o) and (7l) into (67), we find 

T. (S-S.)'H !^V, 3 x S 


(7l) 


(72) 


This is ths snsrgv wwstwd per gram of material swept over by the shock* 
if we now use (3l) end (38) we obtain 


- 4ttjjY 1 T ( S-S.)= - 


- - ZTT J 3 at - V 

- T O- r 

o 


(73) 


v In our ehook wave, *t very large distances, there are two shocks of 
equal strength and therefore the energy loss (73) occurs twice. Therefore 
we find for the decree## of^ energy in the shock ware itself 
^ V/ _ ‘fTT l2r 

ct v 3 c D y ( 74) 

Comparing this with the expression for the energy itself, (66) , we find 

iW w 

77 * ~£rtty7Z T (»> 

o 

which integrates immediately to 


Y/aJ 


(76) 


This result is exactly the same as deduced above by explicit evaluation of 
W, Equation (66). VTe .therefore, have found the decay of the energy by two 


entirely independent methods. Both methods are only applicable if the shape 
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of the pressure pulse remains unchanged* % 

8.8 THE PROPAGATION OF THE SHOCK AT INTERMEDIATE DISTANCES 


The asymptotic form of the pressure pulse behind the shook is usually 
reached only at very large distances when the shock pressure has become ex- 
tremely small* There is an intermediate region, in which the shock pressure 
is quit© small, but the pressure distribution behind the shock is not yet 
linear. Y/e shall now consider this intermediate region, i.e.» we make all 
the assumptions which led to the first approximation considered in the pre- 
ceding section, but we allow, within certain limitations., an arbitrary shape 
of the pulse. 

Since we are interested to apply the theory to the IBM calculations, 
we shall slightly change the procedure used in the preceding sections. In- 


stead of introducing as initial condition the shape of the pulse a t a fixed 
point as function of time, we introduce the pulse at a fixed time t 0 as 
function of the radius y • Hence we write the solution of (17) in the 
form 


2 CO 

,-fc 


(77). 




( t+u) <&' 


(78) 


Then Z as function of y is given by the initial pressure distribution 
at time t Q « Using (49) one finds 


y - 


K 


o 


Ct - 1 ) - ( * rr 


(79; 


The remaining integral is a small term and therefore we may replace dt by 
&T/c . Also <T = ZJn- end £ is constant for the purpose of integration. 
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Thua 


f- M Itjy 


(80) 


Alternatively we might have replaced r in the integral by y . c Cft ) 
and found 


* * - c 0 - ^2 A 

w * y 


(8l) 


Sinoe 2/c^ mast be small compared to y the two equations differ only by 
terms of the order (l/y )^. 

We introduce the abbreviation 




(82) 


which is a oonstant. For convenience we define k in such a way th»t 

/ , . tt-M 

<9 ( = 0 

Then equation (80) takes the form 

In particular, if is the value of 2. at the shock front R 

R w C o (t; +* -y* + ^5 f?4kj 


(82a) 


(83) 


(84) 


Differentiation yields 

d R 

dt 




( 86 ) 


The shock velocity is also eoual to the average value of c ♦ u in front 
and behind the shock, i»e» 


iE .a- [«„+‘+“']* t .[i + ^t] 


(88) 
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Here we have made use of the Equation (49). Comparing the two equations 
for cl R I dt and observing that iZ/t, R is a small quantity, one finds 

( 8 ?) 


Z, id Zj f A. 


2 R 


5 4 

c 


+Jn R j = 0 

s 


or 


2 R d dj 


(88) 


This is an inhomogeneous linear equation in Mqy R and can therefore be solved 
with the result 

5 

di 


JU R * + 4;* 


r 


1 dl di A 


(89) 

I 


Here 2 0 is the value of 2 # at time t • t 0 , when R « R 0 . 
Partial integration yields 


(90) 


This form of the equation is more convenient if g is given as function 
of 2 in numerical form, since it avoids numerical differentiation. 

Either of the two equations gives the shock radius R as function 
of 2 ■ R<r • Now from (lo) j 


s f<, C o r 


(9l) 


and , s ince 


’ * t 


(92) 


it follows that 


** % *ir 6 


( 93 ) 
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Thu* we obtain the shook pressure Ap s a* function of the shook radius* 
Since q, (2) and 2 are finite the right hand side of (90) is finite 
except for 2 # » 0 * Bence 2 g tends to sero as R tends to infinity* 

V 

Let us writ* (89) in the form 






H ** + ff * R * <\ f 2 f ^ -&L \ 

I (89a) 


il 


- a z ft* 


Then the last term tends to O if 2 § tends to saro and we have asymp- 
totically ^ 


where 


A *2. 

a 1 


> = SMC*]**) 




A.R.* = -(jt) 




(94) 


(96) 


This is essentially th* equation deposed previously for eery large distances* 
However, we have obtained here the Constant a in terms of th* initial con- 
ditions and furthermore we have now a definite radius R* in terms of which 
th# radius R should be measured* The contradictions which arose previously 
are therefore avoided* 

, a 

It will be observed that th# asymptotio law (94) is obtained only if 
2 is sufficiently small; if we plot 2 against y % it will have some shape 
of the form shown in Figure 2, approximating roughly th* initial pressure 
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Figure 2 

Initial Pressure Pulse 
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pulse* As time goes on various 2 values catch up with the shock front 

and the range of values which initially extended to 2 contracts and ex- 

o 

tends only to 2 8 • The asymptotic law holds only if 2 g has become small 
compared to 2 Q or more precisely if 2 S is so small that 2 may be con- 
sidered a linear function of y between 2 * 0 and 2 « 2 e • The 
smaller 2 the longer it takes for this 2 to catoh up with the shock 
front, and therefore it may take an appreciable time* before the asymptotic 
law is established, unless the initial pulse is already close to a linear 

u 

function. 

The pressure distribution behind the shock in space is most easily 
- obtained from the Equation (8l) 

r . y+c 0 (t-t 0 ) v r + Qo^o) (9fi) 

c o y 

Asymptotically 2 tends to sero. (We have proved this only for the positive 
phase, but previous considerations have shown that asymptotically the pos- 
itive and negative phase become symmetric* We shall consider the negative 
phase -in sene detail below). At the same time y tends to a fixed value 
y Q • Hence 2 and, therefore, also the pressure, becomes linear in r» 

The length of the positive phase is obtained from (96) by subtracting r 
for 2.0 from R for 2 m 2 S ; i.e. 

L + * ys-yc+ ^f * - ^ (97) 

Asymptotically we have y* — }Yo and 

y 0 ^ 6 e ^-to) : K *^R , (98) 

Inserting this expression into the JCv\. of (97) and using (94) one finds 


JU (R/y.) (99) 

,uiw \ jlmi,— r 

sl "A (rA*) 
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If w© wish to obtain the pressure as function of time at a fixed point, it 
is more convenient to use the Equation (80) instead of (81) 

c o (t-t 0 ) s r - y - r/y 

°o (80) 

Following the same arguments as above one finds that the pulse is asymptot- 
ically linear* The duration of the positive pulse if given by 

t - — f" y. - Vf. + — <^v R/y.l 

+ - 0 0 [ y » °o '} (lOO) 


or asymptotically 


c* A JWfo/yp) 


r-f L 

c o 4 


( 101 ) 


as it should be* Alternatively we may use Equation (83) to determine the 
length of the pulse; with (82a) we find 

l + * -£[* i>»] 


(99a) 


and for the duration 


T. - 


<>o 2 


[g (2 g )♦ Z 8 i r iij 


(100a) 


The shape of the pulse is obtained from (83). If A t is the time 
interval between the arrival of the shock and the arrival of the signal 


£ at the point r, one finds 


«♦ 


At - -21, 

" o ‘ 


£g (<£, ) - g ( 2 ) 4 (5, - 2. ) In rj 


( 102 ) 


8*9 THE NEGATIVE PHASE* DEVELOPMENT OF THE BACK SHOCK * 

'Ye assume that the initial conditions are as indicated in Figure 2, 
sc that the negative phase gradually returns to normal pressure* The value 
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Z m 0 travels with normal sound velocity, whereas all other Z values in 
the negative phase travel more slowly* In faot, throughout the whole region 
in which dljdy is negative, the Z - values have a tendency of catching 
up with each otlier* At some stage a Z- value in this region will overtake 
another and then a shook starts* 

As long as r increases with increasing y the 2 - values are in the 
correct initial order* If r as function of y has an extremum, some 2" 
values have already overtaken others* A shock starts at the boundary between 
the two cases, when r as function of y has a point where both the first 
and second derivatives vanish* 

From (83) follows for fixed t 


dr 

w 


1 - 


C* I 


Co 


Lie 

r 




(103) 


The condition that dr/dy ■ 0 gives 

i* r * 0 

dy ay 

Now from (82) 



(104) 


(108) 


n nd (since is small), dg/dy and d^f/dy cannot vanish simultan- 

eously* Hsnoe they must both be finite in (104) and we may write 


0 


(106) 


If we differentiate (103) once more and put dr/dy - d 2 r/dy* - 0 , we find 


d_ dg • q or d 




u, = 0 

dX ‘ 
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If we neglect small terms of the order Z./c 0 y 
dition c«n easily be reduced to 

d 1 2/dy 1 = 0 

and similarly (10ft) yields with (106) 


end 


^/y $ this con- 


(107a) 




( 106a) 


The Eouation (10?) determines the value II • Ij , at which the back shook 
starts. Substitution into (106) gives the radius r » rj and finally (83) 
gives the time t « . 

The back shock starts somewhere inside the negative phase and for some 

time at least the pressure in the rear and in front of the baok shook differs 

from normal. If 2L. are the corresponding £ values, we have 

T ' A- 

for the shock velocity in analogy to (86) 


it 

d 


;• H'*-' 


z J 


(108) 


where Y is the position of the back shock* In addition we have two equa- 
tions similar to (84) 

<*-* )U 


E? ry- c .C*-^l - 


(109) 


Following the procedure used for the first shook, we find equations similar 
to (88) 

^ ~ 4. Y + + Jfcy\ Y * 0 

x y d 

jLY v _i9 +. im Y = o (uoJ 
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Combining the two equations we find 

^ s0 

Irry V+- ig ( 111 ) 


Eliminating Y by means of (109) one finds 


d 3(2 J - - Z e.)<* 3 / a 

3 C 2 *) - 3 GEf )+(£{: - di<j /dZ f 


dl 


(112) 


This equation requires in general numerioal integration. 

Since both Z-, and l f are negative the shock travels with a veloc- 
ity below the velocity of sound. Benoe the v B lue 1*0 will eventually 
catch up with the back shook* and then the baok shock leaves the material 
behind at normal pressure. Then the equations for the propagation of the 
back shock take a form similar to those for the first shook 


J + 41 * is o, £ - o 

£y A2*- 


with the solution 


w iy 


f 2 ?ia di-|| ju\ 

J «■ ^ 


(113) 


( 11 *) 


Here Ig, Y z are the values of l f and Y at the time when Z r reaohes O 
Just as for the first shock it follows that If approaches sero for 


large Y and one has asymptotically 

2 - r - ~ — • 

‘ TM v/7* 




f 


^ - a 


* * / - ai 


5 ±S A 2. + z 

i r- 


, \ 


J 


(j|) R * 

X- o 


(lie) 

(116) 

(liv) 


APPROVED FOR PUBLIC RELEASE 



APPROVED FOR PUBLIC RELEASE 
VIII - 53 


We have seen previously that asymptotically the strength of the beck shock 
must be the same as that of the first shock. Pence A should be equal to 


B. 


I.e. 





(118) 


Using thfc expression (82) for g we find by Mens of partial integration 


f* 


% n . y + 

A £~ 





(119) 


The last term in the integral is smell and we require therefore 

JL 


J 2 d y * re. ^ ^ X/ya. 


(120) 


Here is the value of y corresponding to £3 * n( * made use of 

the fact that y * R 0 for £ ■ £ 0 « If we choose our initial conditions 
at the time when £ iust vanishes behind the beck shock, they y^ ■ Y 2 
and the integral of £ taken over the whole pulse vanishes. The same is 
true if we ohooee the initial condition# at any later time. 

For later application we shall finally write down the equations which 


result if the initial pulse has linear shape. I«e«, we assume 




y+ k- R « 


-fot t - y & 0 ~ 2 ^ y ^ R * 


(121) 


where R e ie the shock radius and L 0 the length of the positive pulse 
at time t » t Q . £ Q is the value of £ at the shock front, which is related 
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tc the shook pressure by the Equation ( 93) 


> - R . C * 

* " 0( ^ 


( 122 ) 


Differentiating the funotion g(l) defined by (82) one finds with the help 
of (l?l) 

13 £4 !» LjLv 

ii « t . ( (m , 


Neglecting the second term in the bracket, which is email and observing that 
y is for practical purposes constant and equal to B Q , we fiad that tL^/d 5- 


is constant 


_4j 9 Sju Jr* - Jn.8. 
ai <* , 1 


(124) 


Here wo have used the definition (96) for R*« 

For the constant A 2 defined by (9 b) we find 

^ ^ 0 ("S'* * ^ * "Sf - ** = 


%r 




(126) 


The pulse is characterised by the two constants R* and A 2 i is essentially 
the product of the duration 0* L«/c & *nd the shook pressure , multiplied 
by the radius R. More precisely 


a 1 . S. 

°<r h 


The radius R is similarly given by 


(126a) 


R e 


-5f c - 6 Kf*J 


(124a) 
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Here the subscripts o hare been dropped* since these quantities are con- 
stents and apply therefore et any time. 

8.10 TWO PRESSURE PULSES CATCHING UP WITH EACH OTHER 

If the explosion of the gadget takes plaoe at some considerable height 
and the pressure pulse is measured by airborne instruments, the problem a- 
rises whether the reflected shook from the ground oatohes up with the first 
shock, and what happens if it does. 

We assume for simplicity thqt both pressure pulses have reached the 
asymptotic form. As long ee the front of the second pulse has not readied 
the rep r of the first pulse the two pulses behave independently and the cal- 
culations of the preceding sections apply, the mid-point of either pulse 
travels with sound velooity c Q and they keep, therefore, at a oonstant 
distent (apart from geometrical factor* arising from the faot that the re- 
flected Chock has « center different from that of the first s hock i these 
are negleoted in the following considerations). However* the length of each 
pulse increases indefinitely and therefore the refleoted shock will some 
time catch up with the rear of the first pulse. We then have three shocks 
as shown in Figure 3. 

The equation for the first shock, of course, remains unaltered* i»e» 

(126) 

where the suffix 1 refers to the first shook, and the suffix o to some 
specified time, e.g. to the time when the second shock just oatohes up with 
'the first pulse. R* is given by (124). Observing that the length of the 
positive pulse at time t 0 ie equal to l/2 (R 0 - Y 0 ) , we find 

If 

4 O 
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f ' figure 3 

Three Shooks where the reflected shock will 
sometime oetch up with the reer of. the first pulse. 
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Similarly wa have for the back shook of the eeeond pulse 


5 , J r 2 Lz'./z * 


AZ*.JnZ,-i £f %!• 


( 128 ) 


(129) 


For the central shook we have to use the equations devised for the case 
when both the sir in front and in rear of the shock deviate from normal* 
These occurred previously in the treatment of the back shook and are given 
by (110) 


IV 


dY 

d£. 


4 /^vY •» o 


^3 ~ dV + d c Y s 0 


in the first of these eouatlone 

of the first pulse 

v i»e. it is eoual to 

Similarly^in the second eouation * 

Z x -l, 

dY + ^ 4 i =0 

2Y 



dY « 0 

2-Y 

dl ^3 


( 13 o) 


is to be taken from the g-f unction 


These equations may be combined to give 

X V ~ «Li) | 1 J. * » Q 

-is 4 ~xz7ffC j^r 


(13l) 


(132) 


Integration yields 


y i x(i 

■*-3“ *~x to 


) Jk3inE3LLrIBD 
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If we substitute this express ion into one of the Equations ( 13> and 

t 

integrate, we obtain £2 end £3 separately. The result is 


S. , -A 


jlTcirp') 


+ 6 


Z,-A 


■y/z') -iM mn . + p \ 

i 


i^Y/2) 


where 


'lc, 


a. a, - U 7 ) 

infft* / /-*} ! 

r 


(3 _ VU : ^/M + ^30^ /^ .) \ 


N 


TT" .. .M M .. MN - . I 

-UiCY /«*) + /^(Y e /z*) J 


(134) 


(13 fi) 


if ^ (y/s* )>>!,&*. (h*/z* >, the equations o«n be simplified 

-7 ,. ±f? -> ) /jkCv/ftWr./z < T + B 

1 *■' JtuV7w'i 


y . 1 c? - 2 ) <rs£Ziii5tS5^ + b 

fc- , — -» V 30 JIO J 


( 136> 


3 £ 


#nCY/Z*) 


The shock strength of the central shock is related to £3 - £g* According 
to (133) it drops much faster than either the rear shock or the fpont shock. 
However, the pressure level at which this shock occurs tends ^according to 
(136), to a finite £ value. This implies that the oentral shock must either 
catch &he first shock (if 0 ) or the rear shook must catch the central 
shock (if B < 0 ) . 

7 Let us consider first the exceptional ease B » 0. Sinoe 


/ 
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B * O implies that 

(isr) 

Neglecting the small difference between Y 0 and R 0 , these quantities are 
identical with the quantity A defined by (126), which in turn is related 
to the " Impulse" of the pulse © multiplied with the shock radius (see 

Equation (126a)). We conclude that B » 0, if the impulse of the two pulses 

e * 

is the same at the same radius. It will be noticed that R and Z may 
nevertheless differ. Thus,ene pulse may have long duration and low shook 
pressure* the other short duration and comparatively large shock pressure. 

If B « 0, then both £<> and tend to 0 more rapidly than either 

£| or J- 4 * Thus the positive phase of the second pulse and the negative phase 
of the first pulse are eventually eliminated and there remains the positive 
phase of the first pulse, which combines with the negative phaae of the 
second pulse. All this, of oourse, takes some time, until Jtffv (Y/R* ) > > 

Jv% (15,/R >• If R* and Z* differ f rom each other the pressure grad- 
ients in the positive and negative phase will also differ, only the product 
of shook pressure and duration will be the same. If the two pulses were 
completely identical, the remaining pulse woi^ld be indistinguishable from 
either, so that the reflected shook has simply disappeared. 

i 

If B is positive, we require that the impulse in the second pulse 
is greater than the impulse in the first and in this case the central shoek 
catches up with the first shock. Eventually the first pulse has disappeared 
and only the second pulse remains. Similarly, if B is negative, the sec- 
ond pulse will eventually disappear and only the first pulse remain * 

Asymptotically the energy in the blast wave is therefore not equal to 
the sum of the energy in the two pulses but equal to the energy in either 
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the first or the seeond pulse, whichever is the greater* The remainder of . 
the energy has been dissipated at the central shock* This increased energy 
dissipation is duo to the fact that the shock front of the second pulse 
must go through the negative phase of the first pulse* E.g. the energy dis- 
sipation ,iust before the seoond pulse catches up with the first is propor- 
tional to 4(dp a ) , assuming for simplicity equal shook strengths in both 

pulses* The factor 4 arises from the 4 shocks* Just after the second 

* 

pulse caught up with the first, there are two shocks of strength dp 8 
one of strength 2Ap s and the energy dissipation is proportional to 

3 3 3 


and has therefore, increased by a factor 2*f>. 


If B > 0, the radius at which the central shock efttohe* up with the 
first shock is given by the condition ■ ^*1 * After some manipulation 
one finds as oonseouence that a £4 * 0 » the front shook, which is 
now £3 has the same value it would have had if there had been no first 

a 

pulse* Also one finds for the radius R tha condition 


r /r% * 


(158) 


6.11 THE CONTINUATION OF THE IBM-RUN 

The methods developed in the preceding sections have been used to con- 
tinue the IBM results of Chapter 7. The IBM machines give us the pressure 
£p as funotien of the radius nt a fixed time* From these data we obtain 
immediately the function 2(y) and we can deduce the function g (30 de- 
fined by (82). This funotien should , of course, be independent of the time* 

It is shown in Figure 4* The circles refer to the lest IBM cyele when the 
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shock radius is 314.7 (we use here the "IBM-unit" for the radius which is 
19.97 meters). The crosses refer to an early cycle, when the shock radius 
is only 200.2. 

For the numerical evaluation it was convenient to introduce , instead of 
the quantity , 


JUL , 


(139) 


and instead of g 


- *4!*) 

c„ 


(140) 


then the positive phase Can be represented in the form 


45 * 2.2 • + .075 a 2 

J 




NO 




Substitution into ( 69) yields^ \ 


MS'' 




L R • ( m»X + X [2^ + aoTs -] » 


Sf 

ft, ' R 


(141) 


(142) 


where Sp is the value of $ at the shock front. Here we have used the 
boundary condition that at R - 314.7, Ap 8 « .0261 Pq . k/V'T ’ 

As S becomes small, we have approximately y <-*< 

*■ -V 






ft NlUCR/R*) 


i %n R* r ^ i 13*03" 


(R in IBM units ) 


(143) 


Introducing the meter as unit, one finds 

— x >3 I- ' t '~ 

C260»^ 

nr 

to f 


( 144 ) 


”CV, . 


(R in meter) 
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The date from the earlier cycle are better represented by the formula 


g s 1.6 a-f 0.12 e 3 


(141») 




which leads to 


/ 

fa R » -f £ [lv6 -tv 0867 » f j 


(142a) 


and asymptotically, 

* 

K 


(R in meter) 


(144a) 


which differs by about 6 per oent from the previous formula* This is the 
error to be expected, since at the smaller radius the over-pressure is 4.4 
per cent. The data from the larger radius should be more reliable and heve 
been used for all calculations. 

For the negative phase, we get a good approximation by the formula 

( , , — I — ^ 

S • (lit 0.28 s) 4 (3 t ys + a U 6.78 s 

(s<0) 


( 


(145) 


The constants in this representation have been chosen in such a way 
that we get not only a good representation of the numerical data, but that 
also the slope of the function g is continuous at r » o and that the 
integral of the pressure taken over the pulse vanishes. In actual fact, 
the data from the last IBM-cyele are somewhat erratic in the neighborhood 
of s*o. This must be due to an error in the IBM- run, since the data 
from the earlier cycle are quite siwoth. For this reason no attention 
has been paid to this erratic behavior and the condition that do / d S> 
be continuous at s » o has been used instead, in order to get a reasonable 
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function g. 

77e have' seen before (see Kousticn (107)) that the beck shock will start 
at a value of s for which 

* *» l isX -° 



If we use the presentation (146) we find that this condition cannot be satis- 
fied for any negative $• However, this is due to the analytic representation* 
if we go back to the numerical data, we find that the curvature ohange6 sign. 
The exact value of s where the ohange of sign occurs is not easily deter- 
mined,; however the minimum slope can be determined fairly accurately and 
is 8.2* The Slope determines the radius at which the back shock starts 
(soe Equation (106)) 

£m r 8 *2 

s* V«a /minimum (147) 

r ® 1.4 % 10* IBM units = 2.8 x 10® meters 

.Ve shall not be interested in sueh large distances. 

For the duration of the positive pulse we find from (100«) 






JU K - 3 Cs 



(148) 


otf with (142) prtd (14l) 



170.4 

% 



(148a) 


The shape of the pulse at a fixed distance as function of time is similarly 
given by the formula (see Bnuation (102)) 

t e-irs. s> t*, R - $ (s)l 

t L ^ J (149) 
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where s varies from its value at the shock front to at. the pres- 

sure minimum and back to tero. In order to obtain the time in seconds, we 

should use for c_ the "IBU-value" ■ 17. 3 8. 

o o 

The peak pressure versus shock radius curve is shown in Figure 5, using 
"IBV-units". For comparison some points obtained, if we stop the IBM-run 
at an earlier instant, are also shown* 

The .duration versxis peak pressure is shown in Figure 19 of Chapter 7. 
The equations for the shape of the pressure pulse are used in the next 
chapter, to obtain data for the two combat bombs* 


( 
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Figure 5 

Peak Pressure versus Shook Radius 
— —— Seni-Aooouatio theory plotted at R 
y Semi -Aoc oust lc theory plotted at R 

— IBM Result* 
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CHAPTER 9 


THE EFFECT OF ALTITUDE 


K* Fuchs 


9.1 IKTRQDUCTIOH 

0n# of the easiest maasursasnts of the effeot of a nuclear bomb in con- 
bat conditions is a measurement of the pressure pulse by means of airborne 
gauges. The instruments used for this purpose are described in Chapter IS. 
Safety considerations for the planes from which these gauges are released 
make it impossible to get close to the bomb. Furthermore, since these gauges 
are released at great, height at about the sane time the bon* is released, 
but are attached to parachutes, they will still be at a considerable height 
when the blast reaches them. 

The information received from the air^prne gaugee was in faot the only 
quantitative information on the bombs droj^i over Hiroshima and Hagesaki, 

: r , ■ T 

until scientific te^ms could enter these'"*'' 


n 


towns after the oessation of hoe* 


height of about 30,000 feet, where 


tilities. The instruments recorded at a 
the atmospheric pressure is considerably reduced. They were at a distance 

.<■; ’ y 'r * ■' 

of 36,000 to 4o, 000 feet from the bomb. 

Such large distances can be oovered by the semt-acouatla theory dial * 
oped— in- $hfl prft c ed iag_jha£ter. However, we have to estimate also the effeot 
of the change of pressure and temperature with altitude in order to inter- ' 
pret the records* 

We shall make the assumption that tha energy is emitted uniformly in J. 
all directions. This assumption appears eminently reasonable, since the 
altitude effect will become .pronounced only after the shock wave has decayed 

IX -1 
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into a pressure pulse of length smell oompared to the distance from the bomb* 
If in addition we assume the acoustic treatment, we obtain the Rayleigh 
correction factor for the over-pressure 



where C j are density, sound velocity and over-pressure at a given height 
and J> c are the same quantities for a uniform atmosphere. The duration 

Qf 09 |0 

of che pulse a.id the energy in the pulse are unchanged. The latter follows 
immediately from the fsot that in the acoustic theory the energy in the pulse 
is constant* 

We have seen in Chapter 8 that the energy-dissipation in the pressure 
pulse is important. Clearly, if the_ft aerev dissipation is taken Into account , 
the pressu re pulse at a &l,ve.a_altit«die Will. AeasuadL-jaat.- on ly on the c onditions 
at t hat altitude , but elsfe of the pul se. For a linear 

pulse we Shall find that the correction factor is composed of two factors. 

The first is identioal with (l). The second factor represents a change in 
the form of the energy dissipation term, which in a uniform atmosphere is 
(log (F/r 0 )) . The log is replaced by the integral dR/R, weighted with 

a. 

a function which depends on the variation of the atmospheric conditions. 



The duration of the pulse, which in the acoustic theory is Unchanged, 
also undergoes a change, if dissipation of energy is taken into account* 
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We shall find for a linear pulse 



0 



so that 



O 


Ct» 
— ***»• 

c. 





( 4 ) 


(5) 


The duration of the pulse therefore increases with altitude* 

In Station 5.4 these formulae are generalised for an arbitrary shape of 
the pulse* 

Compared to the acoustic theory the pulse is therefore les3 strong., but 
longer* In principle It ie therefore easy to oheck the theory by considering 
the distortion of the pulse ahape* which ia approximately independent of the 
energy release* 

However, if we try to do so, we are confronted with the embarrassing sit- 
uation that the Hiroshima record and Nagasaki record oontradiot each other* 
This is the more mysterious, since they were taken at exactly the same alti- 
tude and approximately the same distance from the explosion* lienee, we should 
expect the shape of the pulse to be approximately the same in either case* 

In actual fact the Nagasaki record gave a higher peak pressure but exactly 


the same duration as the Hiroshima record* 

If this result is real, it cannot be explained as an altitude effect, 
unless the variation of temperature and preseure with altitude differed com- 
pletely in the two cases* Neither could it be explained on the basis of 
any effect which obeys the usual scaling laW3* 

If we trust the peak pressure measurements, the Nagasaki tomb was about 
5 times as powerful as the Hiroshima bomb, in moderate agreement with the 
factor 4 deduced by Penney from the blast damage* If the duration of the 



\ 
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pulse is used, the two bombs were equally powerful. 

The Hiroshima record gives a shape cf the pressure pulse in fairly close 
agreement with the theory* the Nagasaki pulse, however, is much steeper than 
expected. 

The values of the nuclear energy release derived from either theory, 
using the peak pressure or the duration of the pulse as criterion are given 
in the following table: 

Table 9.1 



Energy release 

cf combat bombs 

ir tons of TNT 



peak pressure 

obtained from 

1 

\ 

— 

duration 

Acoustic theory 
Theory with energy 

dissipation 

36.000 

51.000 . 

Nagasaki 

i 

■■ 

20,000 

11,700* 

I 

\ 

Ao oust ic Theory 
Theory with energy 

dissipation 

1 

| 7,700 

! 10,600 . 

Hiroshima 

19,500 

11,400 • 


9.2 ACOUST IC TH EORY 

We make the following assumptions: 

(1) The pulse is weak; mere precisely A 

(2) The length of the pulse is small compared to the radius to 
which it has penetrated* 

( 3 ) The length of the pulse is small compared to a distance over 
which the atmosphere changes appreciably. 

(4) The energy is emitted uniformly In all directions. 

An important consequence of the assumption (3) is that the distortion of 
the pulse is confined to a change of scale only. In particular the ratio 
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of the mean square of the pulse to the square of the shock pressure Is un- 
changed. This is not true if energy dissipation is taken into aoeeunt, un- 
less the pulse has the asymptotic linear form. 


(a py 


( A f0 


( a k.) 


( 6 ) 


Fere the suffix o refers to the seme quantities in a uniform atmosphere. 

In the acoustic theory ell signals travel with the weloeity of sound* 
If at time t - o, the end of the positive phase was at r • r 0 , then 


at time t it will have reached the radius r, where 

-K> ft* 

Ct« 


-/ 


Similar ly for the head of the pulse 

h,4 L <V- 

cU 

cC >") 


h.- 

/ 


V + L 
0 o 

where L q , L is the length of the pulse *t time t » o, and at time t. 


Subtracting the two equations we find 
-*t, + L o + L 

iiv - / ^ 

J co; 




(7) 


These integrals, however, are identical with the duration of the pulse* 
Fence we find that the duration cf the pulse is constant. 


o - e„ 


( 6 ) 


The energy per unit solid angle -in the pulse is given 

X 


ZO 


R j A |>u. cLt 


(9) 


to be integrated over the path of a volume element. For a weak pulse, we 
c an integrate keeping the position in spaoe fixed* Also j yQC'u. * Hence 


w ~ z R 5 C-AfO* sjj* c* 


(10) 
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The factor 2 arises from the negative phase* 

Sinoe in the acoustic theory the energy' per unit solid angle is constant* 
it follows with (6) and (8) that 


hso 


J> c 

o o 


( 11 ) 


This is the Rayleigh correction for the pressure. 

9.3 THEORY INCLUDING ENERGY DISSIPATION 

In order to include the effect of energy dissipation, we shall make 
one additional assumption to those mentioned in Section 2. That is 

(5/ The pulse has linear shape. 

This assumption is at least approximately true for the positive phase in 
the region in which we ere interested. Since the positive and negative 
phaeasdo not interfere with each other, the assumption la therefore approx- 
imately correct. 

One of the shock conditions is 

r * * ( 12 ) 

in view bf assumption (3), we can treat the history of any given volume 
element on the basis of a theory in a uniform atmosphere. Then the equation 
(12) is satisfied throughout the pulse 

e 

(13) 

The energy per unit solid angle is again given by (9 )* for a linear pulse 
we find 

W , R A J A u it ■(& R 11 j 
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The energy dissipation in the shook and back shook can be expressed in the 
form (compare Equation (74) Chapter 8) 


d R 3 f J c ♦ 


(It) 


Combining the two equations we find 

a ^ [r 2 ( A )>/ ejj> cj + R 1 (A f>/ //• 4 = 0 


The shock velocity is given by the Hugoniot condition 


4 B = C + 

dfe 


2^c 


(16) 


(17) 


So far we have been concerned only with the local conditions at the point 
where the pulse happens to be. We require one more equation, in which the 
variation of atmospheric conditions appears* 

For this purpose we observe that the average velocity of the pulse is 
given by the sound velodity* This follows from the feet that the velocity 
of the back shock lags as mubh below sound velocity as the velocity of the 
front shock is above sound velocity, and that tha pulse is symmetric. The 
latter was shown in Chapter 8 by means of the general argument that the 
total displacement of any volume element decreases as i/r^ and tharafore 
asymptotically 

ci t ~ J Ajz cLt'x 0 


Thus the midpoint of the pulse, which is elso the end of the positive phase 
travels with sound velocity. If & is the duration cf the pulse at radius 
R, we have therefore 



cl \ — Q 
t 



where t is the time when the shock front reaches the radius R. 
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We near have ell necessary equations end proceed to solve them* Dif- 
ferentiation of (18) yields 

cU a X - && 
cLR c 


(19) 


Obt erring that the second term ie small# we find 

1 .s c -+• c* d Q 

d. t R 

Comperieon with (17) yielde 


( 20 > 


d © 

~d R 


oc t 


(2l) 


This equation can be used to eliminate from (16), with the result 


(22) 


Integration of this equation gives 

H jot B ( ) « const - ^ to 


(23) 


A second integration yields 

.2 _ o* ! i£ 


0 


where r c is another constant of integration* 

Substitution into (2l) givee the shook pressure 

r~ir“ / T S' R 1 r\ 


A b .. ufii/r r _t 

' dR / I / C 1 

t n. 


dLR 

"k 


( 24 ) 


(25) 


In a uniform atmosphere of density ^ and sound velocity e 0# we find the 
asymptotic law of Chapter 8 

, . . ■> BC-40^* 

c ^ ' «tp^rtR/r.s 


w 

I 


(26) 


thus r a should be identified with the rediue R* in the asymptotic law 
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for a uniform atmosphere* The altitude correction factor ie the ratio of 
(26) and (26> 

A ^5 - / -J°C 06 

C&f&h 9 ( 27 ) 

where e/e. ie the altitude correction factor for the duration, which ie 
given by 

6 - 

© 

o 

9.4 ALTERNATIVE DERIVATION 

We shall now give an alternative derivation of the formula* given above. 
Thie derivation followe closely that given in Chapter 8 for a uniform atmos- 
phere* 

In this derivation we use the concept of 2. - signals. In e uniform 
atmosphere X reins ins constant along a characteristic. If atmospheric con- 
ditions change, this ie no longer true* We shall make the assumpticn that 
all £ - values change with the distance from the center of the explosion, 
by a common faotor |j\ , the same for all X • This would seen to be a reason- 
able assumption an long as the length ot the pulse is small compared to the 
distance over which the atmosphere ohanges appreciably.. However, the complete 
equivalence of the assumption about the behaviour of the 2L- signals with 
the approximation of a small pulse, is only .demonstrated by the fact that we 
shall obtain the same results as in Section 9*5. For this reason we preferred 
the previous derivation, quite apart from its simplicity. 

The derivation which we shall give below has on the other hand the ad- 
vantage that it applies to an arbitrary ehapa of the pulse and allows us to 
make s prediction about the shape of the pulse. The pulse is distorted by 
the change in the 5L - valytes by the factor V • The shock pressure, however. 
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is changed in addition by another factor, expressing the fact that the energy 
dissipation has changed and, therefore, the 2.- value at the shock front is 
no longer the sane as before. This second factor has no influence on the 
shape of the pulse, except insofar as it outs off the pulse at a different 
2 “ value 1 The derivation of Section 9*3 gave us only the product of theae 
two factors. Hoe we shall obtain each factor separately. 

Furthermore, in the derivation of Section 9*3, we had an arbitrary con- 
stant of integration r Q , which represented the radius at which the atmec- 
pheric conditions could still be considered uniform. We should have expected 
that we should take r 0 » 0* However, in order to get correct results 
we found that r c should be identified with the radius R * in the asymptotic 
law for a uniform atmosphere. The treatment given below will resolve this 
difficulty. 

Our assumption can be stated in the form 

s . T k 


where 2. it the value of £ at the radius r ■ r and V ia a funotion 

o 

of the radius. Inuation (29) holds along the path of a characteristic, 

whieh goes with velocity c ♦ u. The excess velooity of the characteristic 
above the looal sound velocity is given byes if a t time t ■ t°, 

the characteristic is at r « r 0 , the equation of the characteristic is 


, ... r ax. _ f 
' J C vol Yfc " J 




(30) 


'J 

O 


c.H 

C 1 *, 


JL A. 


(31) 
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Jfere r o is assumed s fixed redius end t is then evidently a function 
of 2. °* K appears only in a small term and instead of integrating over the 
path of a characteristic, we can integrate over the path of the shook front. 
Then the integral is independent of 2, °* 


At the shock front we have 

R 


t - t 


- f ¥ - 


/V 


^5 K C R ) 

c * 

6 


(3?) 


Differentiation yields 

dt , ± J. 
dR ' c v 


o CK 

r 5- 


- j 


ikr 

d7// dR C* clR 


(33) 


The shock velocity is given by 

dR .. f = C + ** 

d* 


c + 


2.R 


2. R 


- c+ <a< c. 
2 C" 


dj< «C 
d R Z 'S 


(34) 


Comparison yields 


(K- 


i£ 

dZ‘ 


) 


d 2 - i . d S' d K 

d g + * * R 


(36) 


Let us introduce instead of R the variable Z defined bv 


dZ 

2 


cF 


d a, 

a 


Z - )u 

c o 


(36) 


then 


K - -W ( *-/•%) 


(3 7) 


and 


Cu z/z. 


1 v 

C au C y» » 

q d i \ ci 2L s i -L 

°< TV ' aZ 1 

V : ■ 


£s__ = 0 

Z 


(3 8) 


In this form, the variation of the atmosphere has been eliminated. Hence 

is the same function of Zj E. for a uniform atmosphere is of R» 
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i* c • 


2 (Z) 

OS 


. ( 39 ) 


where ^ (R) Is the for a uniform atmosphere* In order to calculate 

os S 

k we use again the energy dissipation* The work done by the pulse per unit 


solid angle is 
to 


/"pud* = Ir’joc fcrfit 


However, for fixed r we have from (30) 

d* = - C°<£ -4|* ) ±1" 

C o «** Z 


(4c) 


(41) 


Hence 


AjO - 


k>« /‘cr^ ( ^ .i|*.)dx* 


(42) 


We assume that the back shock has not yet started; then the integral 

_ « 

extends to the minimum value of 2~ in the negative phase and then back 
to £ - O . Differentiation of (42> yields 

z" 

is = + cjjd£c_ (« r ( z-/d2, 

d r d r Ki > c V | <* ft ; 

, /t- 6 c e 


(43) 


In the first integral the negative phase does not contribute and the pos- 
itive phase gives • Utilising (35) we find 

dto , d^j°c ) y) . CZ ° f 

duR ,iR k*Pc ft 4 * TR 


(44) 


The energy dissipation occurs at the shock front -nd can be expressed in 
the form (see Chapter 8, Foust icn ( "3 ) ) 
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d-UJ 

d ft. 




tin 3 

eR 


(46) 


Using the definition (3l) of K it is found that this is identical with 
the last term in (44). Hence the first term in (45) vanishes. This requires 


k 


and 


K -'UCT.jZU f 


<Lk 

/c 


o 


The shock pressure is piven by 

With (40), (47) end (39) we find 

cz) 

rs R os 

and the altitude correction factor is 


Ah , nr 2 - 

J — / » gn — — 


CZ) 


CA|> S ), Vy2c # Z #s .O0 

The equation (49) can clearly also be written in the form 

^ "• a 


(46) 


(47) 


(48) 


(49) 


( 60 ) 


(40a) 


where Ap so (r) is the shock pressure in a uniform atmosphere. 

4 

From (3C>) follows that the arrival^f the signal 1° at the radius 
r is delayed with respect to the shook by the time interval 

At = +^T K (») [I* - 2*] (51) 

' = z? JU Cz/*.)p'-2'] ■ 


APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 

II - 14 


For a 

uniform atmosphere, we can replace X° by £ 

and 2 by R* 

Comparison 

With 

Bquaticn (102) of Chapter 8 shows that 




A. t*cz) + I u a. = (*-) 

0 1 0 ° 


(63) 

where 

g(l) is the function Introduced in Chapter 8* 

Hence 



A* = fi [y C2*) - % (Z ) 4 cz s -z.. ) u zj 


(54) 

In view of (39) this can also be written in the form 




Ai = A t C2-° J ) 

0 y ^ 


(54*) 

V 

where 

At 0 is the corresponding function for a uniform atmosphere* 

The 


duration 6 is obtained by putting £° * 0. Then 

e*£ [»<**>♦ z:j»z]*«.«;> 


The eouations we have derived no longer contain the radius r Q except in 
2* It is clear that the equations should be independent of r 0 provided 
we cbm find a value of r 0 sufficiently large that the semi-acoustic theory 
apnlies and at the sane time sufficiently small, so that up to the radius 
r Q the atmosphere can be considered uniform* With that assumption we pro- 
ceed to get rid of r Q in the expression for Z* We write (45) in the form 


k 


Xtv Z 7 Jun n. + 


f 

J 


\ 





«■/ 

fu 

0 


d v 

A. 


( 66 ) 


Since the atmosphere is assumed uniform for r<r Q , the integrand vanishes 
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In this region end we can extend the integral to ft * o 

R * 2 


At -Ar, j {b ~ 


( 67 ) 


It reoeins to show that our equations coincide with those derived in Section 
9.3 for a linear pulse. For this oasa (see Chapter 8, Equation (94)) 


~ viag(R i^y 


(58) 


It should be noted that the R* here hss nothing to do with the r 0 which 
we had before, but Is determined by the properties of, the pulse ns shown 
in Chapter 8. 

From. (57) follows 
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O' - J*t 

If the atmosphere oan be considered uniform up to the radius R* this can 
be written in the form 

R __ 

djL 


X~(Z/R') 

f * 


(so) 


Substituting (58) and ( 60 ) into (5o) , we obtain the formula (27), (28), with 
the correct radius R* • It is now clear that the necessity which arose 
previously of identifying r c with R* was artificial, introduced by the 
limitations of a theory which assumes the pul3e to be linear throughout. 

The results which we have obtained can be expressed in terns of simple scale 
changes* Define 
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tinoe R is proportional to the altitude , it can also be written as 


*->/# m-'Yi 


( 62 ) 
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and depends clearly only on the altitude, but not on the properties of the 
pulse* Then the altitude correction, is taken care of by changing the radius 
scale by X , the pressure soale by and the time scale remains 

unchanged. The factor X appears in the pressure scale because J>c 1 * RAp 
and the 21-scale is changed by the factor k Hence the scale 

factor for RAp is the acoustic correction factor V^c /^ c o , but the scale 
factor for Ap itself contains the length soale factor X (see also Equation 
(49a) ) • 

9.5 EVALUATION OF ALTITUDE CORRECTION FACTORS 

The temperature T is a linear function of the altitude up to about 
10,000 - 12,000 meters. It is given by 


; 1 - * h 


(63) 


a x 2*26 * 10“ 5 / meters 


(64) 


It follows that the density ^ is given by (See for example Durand, 
Aerodynamic Theory, Volume 1, pp 219-223)* 


£1 a. (1 - a h ) 
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( 66 ) 
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k : 2.128 


( 66 ) 


The ‘sound velocity is proportional to the square root of the temperature 
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Substitution into (62) yields 

* ah 

XtvX = / ( 1 I & J. 
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( 68 ) 


Usually ah will be sufficiently small so that an expansion of the integral 
can be used 

U\._ (nX) ^ ( .h ) 2 * (»n) 3< ~ 

f 2 1 ‘ 3 1 ” {691 


Also 


k = ^25 - ( I-** 

" /c 


(70) 


9.6 APPLICATION TO ITIROSHIMA AND NAGASAKI 

* - — n - - rr - — - 

• D 

The experimental data for Kiroshima and Nagasaki are given in Chapter 18. 

Using the semi-acoustic continuation of the l3K~run and the formulae for the 

* 

altitude effect, we can determine the nuclear energy release required to give 
the observed peak pressure or the observed duration of the pulse. The energies 
obtained in this way nave been recorded and discussed in the introduction- 
The altitude correction is quit# appreciable. The soale factor X for 
the length is 2.527 and that for the pressures is 1*442; i.e., the* pulse 
at that altitude was the same as that expected in a uniform atmosphere at 
2.527 times the distance, except that the pressures have to be increased 
by the factor 1.442* 

The predicted pulse Shape is shown in Chapter 19 , Figures 5 and 6 , 
assuming 11,000 tons energy release for Hiroshima and 30,000 tons for 
Nagasaki* The Hiroshima record agrees quite well with the theory. The 
experimental Nagasaki record, on the other hand, looka quite different from 
the expected curve* This, of course, is not surprising, in view of the 
large discrepancy in the energy release estimated either from the peak pressure 
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or from the duration* 

A strange feature of the Nagasaki record is the fact that the shape of 
the negative phase agrees very well with the prediction, as is shown by the 
dotted curve* which is identical with the predicted ourve except that it is 
displaced by 0*4 seconds* 

At present we have no explanation for the shape of the Nagasaki pulse* 
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CHAPTER 10 

THE MACH EFFECT AND THE HEIGHT OF BURST 
J. von Neumann and P. Heines 

10 >1 GENERAL CONSIDERATIONS ON THE PRODUCTION OF BUST DAMAGE 

I 

Bomb damage to structures is largely caused by reflection on the structures 
of the shock wave generated by the bomb, and in the case of long blasts such as 
caused by an atomic bomb, by the ensuing blast wind. In this article we will 
consider the pressure criteria which determine the height at which the bomb 
should be burst so a 3 to maximize the area of blast damage. Hie problem of 
maximizing the incendiary effect of the atomic bomb is treated elsewhere in 
this volume, * • 

the first conclusion one reaches in studying the problem of blast damage 
is that the problem is extremely complex and can only be solved in a statistical 
or average manner. This, is so for two reasons; first, the detailed description 
of a military target can never be completely given, and second, the complete 
analytical solution of even such a relatively simple problem as the behavior of 
a blast wave incident on a wall at an oblique angle has never been obtained for 
all angles. As we shall see later, a solution of the basic problem of shock 
reflection from a rigid wall can be obtained by a combination of theory and 
experiment. This solution is, however, not readily adapted to yielding the 
effect of blast in better than an average sense in a more complicated situation. 
As to the detailed description of the target, not only are the structures of ^ 
odd shape, but they have the additional complicating property of not being rigid 
This means that they do not merely deflect the blast wave without absorbing ' .1 

X *? 

A 
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energy from it, bat take a tariff on the blast at each reflection. 

In addition to being weakened by destroying structures, the blast may, 
of course, also be weakened by imparting kinetic energy to the debris. Ihe 
removal of energy from the blast as it does its job decreases the blast pressure 
at any given distance from the point of detonation to a value somewhat below 
that which it would have in the absence of 'dissipative objects such as buildings. 
Ihe presence of such dissipation makes it necessary to consider somewhat higher 
values of the pressure than would be required if there were only one structure 
set by itself on a rigid plane. The ideal would be to determine for a given 
explosion and given target configuration the required pressures by treating 
these dif fractional and other losses theoretically and from first principles. 
Even this, howeverj would probably be too difficult in any well defined but 
realistic special case, and, furthermore, the necessary statistical theory to 
derive valid ffrerage results for the actual irregular and variable target con- 
figuration is not known. The next best procedure would be to derive from theory 
or from experiment the lethal pressures for a given isolated structure, to 
average this properly over various actual structures, and then to compare this 
with the average pressure level at which the damage under consideration has 
been empirically found to occur. The ratio of these two pressure levels would 
then express the losses in question. This procedure is not practical either, 
mainly because the first mentioned pressure criterion is not sufficiently well 

P' 

known or sufficiently reprodapibly def tried. These influences must therefore 

*■ ■ 

be accounted for in qualitative ways, not on an absolute basis, but rather in 
the sense of comparing them, in two situations - one, the situation of actual 
interest, and the other one, a standard where the pressure level corresponding 
to the observed damage radius is empirically known. 

v- 

Ihe pressures which are actually inflicted upon a given structure will 
have been amplified by reflections from its own or nearby surfaces and decreased 
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by diffractions around openings and corners. The latter forms part of the 
losses referred to above. The former may cause local increases in pressure 
which will, under suitable conditions, be quite considerable. For over- 
pressures of less than, say, 10 pounds per square inch, the acoustic theory 
applies in the main satisfactorily, except for rather glancing reflections. 
According to it any reflection, head-on or oblique, doubles the overpressure 
at the surface. (Concerning the limitations of this assertion - cf. Section 10.3.) 
It is not difficult to find geometrical arrangements where as many as three 
such reflections superpose their amplifications without any losses, e.g. the 
phenomenon of a shock running into a 90° corner.' (Figures la, lb, and lc.) 

Thus, in such a case there is a local increase in pressure of four times the 
initial overpressure in the shock, even in acoustic theory. (The exact shock 
theory gives even higher increases, cf. below.) The case cited is a special 
example of the more general problem of the pressure increase which may be 
obtained when an acoustic shock runs into a corner having an angular opdfting 0 
at an angle of incidence oC. In this general case it can be shown that the 
pressure increase to be expected is np where p is the overpressure in the shocks- 
wave, n is the number of shocks that have traversed the region; n is determined 
by the angle of the wedge 0 and the angle of incidence oC . Reasoning from 
the fact that infinite pressures are expected in a convergent cylindrical shock 
as: 

q — — > 0 for o< se 0/2, n > 00 

The actual pressure amplifications that can be attained are limbed by the 
rarefaction wave from the edges of the necessarily finite wedge. 

These local increases are real and well-known: indeed the blast is noto- 
riously erratic, and frequently affects some parts of a structure much more 
than other closely adjacent ones. In spite of its reality, it is clearly 
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hopeless to attempt to follow this effect into all its ramifications. On the 
other hand this is not absolutely necessary, since what is really needed, is an 
average statistical theory. Now these local reflection amplifications and the 
opposite diffraction shielding effects which necessarily accompany them in 
neighboring areas, are clearly all perturbations of the main blast phenomenon, 
caused by the irregularities represented by the target and other structures. 
Consequently what must be estimated is this: Which formations can be dis- 
regarded as local irregularities not materially affecting the main, average 
evolution of the blast wave? 

We do not possess an exact theory of this phenomenon. There exist, however, 

good guiding analogies with acoustics and optics - that is, with linear wave 

theory. It is well-known (cf., e.g. Wood, Physical Optics) that irregularities 

with dimensions of less than 1/16 of a wave length leave reflection "perfect” 

from the "optical" point of view. Probably a good deal less than this, as 

much as 1/4 ease lengthy will not affect average intensities significantly, i.e, 

they are negligible in the sense outlined above. The linear dimensions of 

typical houses are of the order of 30 feet to 50 feet, they may conceivably go 

up to the order of 100 feet. Hence they are negligible for wave lengths of 

120 feet to 200 feet or over, and even in extreme cases for wave lengths of 

400 feet or over. In the case of a blast wave, a certain difficulty is caused 

by the absence of a well-defined wave length in combination with departures from 

linearity for higher shock strengths. It is clear, however, that Fourier 

analysis, as well as any other possibly preferable decomposition procedure, will 

have to assign to the length of the main blast the essential role of a wave 

length. This main part of the blast wave is the so-called positive phase of 

the blast. Hence, it is felt to be justified in treating houses as statis ti- 
lt* 

cally irrelevant perturbations, if the positive phaset has a spatial tension of 
at least 120 feet to 200 feet, or perhaps at least 400 feet. It is clear that 
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these arguments are loose and rather qualitative, but they probably do Justice 
to the main features of the situation. 

The length of the positive phase of the blast due to an explosion caused 
by W tons of TNT - or any other explosive of equivalent blast energy - varies 
slowly with the shock pressure level at which it is taken, and is proportional 
to wV3. i n the significant region, which will turn out to be at 5 to 10 
pounds per square inch, the duration of the positive phase is about 0.025 seconds 
for W * 1 ton. 

The shock velocity in this range averages about 20 per cent higher than 
sound velocity, i.e. it is about 1.2 x 1100 feet per second - 1320 feet per 
second. Hence the length of the positive phase is 0.025 x 1320 feet = 33 feet. 
For a nuclear explosion, say 20 x leP tons TNT blast equivalent this becomes 
(20 x 103)V3 x 33 feet = 900 feet. 

To sum up: For the nuclear explosions under consideration (and actually 
even for much smaller explosions) houses and other obstructions of comparable 
dimensions may be treated as small perturbations which do not appreciably 
affect the main evolution of the blast. For a later application it is useful 
to point outjthat houses, quite apart from their established "small" sise, are 
also a feeble overall influence because they cover only a small part of the 

! 

ground. Even in "built-in areas" they hardly ever cover more than 25 per cent 
of the ground, and therefore^the progress of a blast wave along the ground 

takes place in the main over a smooth surface. *' 

>. 

y 

This being understood, the one obstruction which cannot be disregarded on 
these counts will now be considered since its dimensions are more properly 
described as infinite. This is the ground itself. In other words: In follow- 
ing the main evolution of the blast, the reflections and diffraction of a blast 
from everything else but the ground can be neglected, but its reflection from 
the ground must be taken into account. According to these principles, 
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the ground can be treated as a plane reflecting surface as long as its formation 

is plane from an average topographical point of view. (Focussing, shielding, 

% 

and other gross effects of hills, valleys, etc. are well-known, but we do not 
propose to consider them here.) Also, since its density is 1,000 times that 
of the air, the transfer of energy through the air into the ground is negligible 
and it may therefore be treated as a rigid reflector. 

There is a great deal of direct experimental evidence confirming these 
views (cf., e.g. data obtained at Woods Hole and Princeton on high burst). 

Thus we have reached a standpoint where the reflection of the blast frfim 
the ground, idealized as a rigid, reflecting plane, is taken fully into account, 
and all target structures are then viewed as immersed into the average pressure 
field thus produced. Two situations (in two different explosions, st two 
different positions) far two targets will be viewed as equivalent, if the 
surrounding pressure fields obtained in this manner appear to be. equivalent. 

Thus far nothing has been said about the characteristics of a pressure 
field which determine damage. A static pressure, i.e. one which lasts forever^ 
will damage a given structure in a given way, if it exceeds a certain minimum 
value which is easily determined experimentally. A pressure which lasts 
only a shorter time, say t, will clearly have to exceed a higher minimum value 
p t ,in order to cause comparable effects. Clearly p^ will not differ signi- 
ficantly from p ao , if t exceeds some time X characteristic of the structure 
in question, being essentially the elastic half period of an elastic structure, 
or the time in which irreversible basic deformation occurs in the case of an 
inelastic structure. For ordinary houses, etc. this X is of the order of 
15 mil li seconds.** For the 900-foot long nuclear explosion blast (cf. above), 
which therefore lasts about 3/U second, this limit is clearly far exceeded. 

In actual blasts p varies with time: p = p(t). When the duration 
experience shows that the significant damage criterion is the impulse f pdt. 
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For a duration which is > x , the beginning of the peak pressure curve up to 
X is in the main p(o), and it is known that the significant quantity from 
the point of view of damage is the peak pressure. Thus for nuclear explosions 
(cf. above), the damage criterion is the peak pressure viewed statically and 
the situation boils down to something fairly simple: reflection from the ground 
and then only the peak pressure to consider. 

10.2 THE HEIGHT OF DETONATION AND A QUALITATIVE DISCUSSION 

OF IKE MACH EFFECT 

A bomb detonated on the ground is certainly closer to the target than an 
air burst. For an air burst bomb, the height of burst has no profound effect 
on the blast received at a point which is several times the height of burst 
from the bomb. At distances which are small or of the same order as the height 
of burst, the fact that the bokib is air Imrst has a profound effect on the blast 
characteristic. In the immediate neighborhood of a ground burst a target suffers 
extremely high pressures which it could not receive if the bomb were air burst. 
For big charges this is, however, not an advantage because it means that the 
immediate neighborhood would be destroyed more radically than is necessary, and 
the energy so wasted would not be available elsewhere. In other words, the 
immediate neighborhood would be overpulverized. 

t 

It is actually practice when using smaller bomba 'to try to make a full hit 
and detonate it exactly on the target. When using larger bombs, however, one 
does not try to make a full hit, but rather to wreck an entire area. In this 
case there is some point in trying not to detonate it on the ground, but at a 
certain altitude. In this way the area nearest the bomb is not overdestroyed, 
more energy is not used on it than is required, and the air burst does not 
permit nearby structures to shield those which are remote. 

For nuclear bombs there is a further reason for detonation at an altitude. 
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It is not desirable to let the enormous temperatures (about 10^ to 10^ degrees) 
immediately around the explosion get in contact with matter in bulk, in 
particular with the ground. Such contact would allow the wastage of much 
energy in evaporating the earth. 

An example of the loss of energy from the explosion by evaporation of 
matter in bulk can be found in the disappearance of the tower which carried the 
atomic boob at Trinity. Figure 2 indicates the asymmetric shock-ball, as 
actually shown in several Fastax photographs of the explosion. It shows that 
the progress of the shock wave has been visibly retarded where it intersects 
the tower, presumably because of the loss of energy incurred in evaporating 
the steelwork. 

There is a further advantage in air burst, a discussion of which constitutes 
the bulk of this report. An air burst is accompanied by certain forms of blast 
reflection which would not occur in the case of a surface burst. 

As an orientation let us first consider the case of a bomb detonated on the 
ground. It might at first be thought that because the shock has a hemispherical 
shape and always touches the ground at right angles there would be no reflection 
from the ground. Actually, this is not true. If the ground were an absolutely 
rigid reflecting surface, then the energy normally transmitted to the lower 
hemisphere will not disappear but will be sent into the upper hemisphere in 
coincidence with the energy normally sent there in the absence of the ground. 

In other words, a bomb detonated on the ground is equivalent to two bombs insofar 
as the blast in the region above the ground is concerned. Now twice the charge 
weight means, if one considers peak pressure, that all pressures are the same 
for the ground burst charge as for the same charge burst in free air if all 
distances from the charges are in the ratio of the similarity^ factor 2"^. 


The similarity factor comes from the fundamental hydrodynamics 1 equations 
in which the distance, R, always occurs in the combination R/wV 3 with the 
blast energy W; i.e. the unit of length is determined by the blast tonnage, 
The time scale i3 similarly extended. 
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Figure 2 
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Comparing pres suras at the same point, in the region in which pressure decreases 
as 1/E or 1/R-^ 2 (R - distance from the bo mb) peak pressures are obtained in the 
presence of the ground which are 2^/3 to 2^/^ higher because of reflection. 

Ihe 1/R law holds roughly at large distances, i.e. low pressures, but the variation 
is like 1/R-*/ 2 in the region of interest, at 5 to 10 pounds per square inch. 

It is clear that one would not get the indicated pressure increase if there 
were cratering because cratering means that some energy does go into the lower 
hemisphere. In other words t If the bon* is burst on the ground the peak 
pressure is increased to anything from 1.26, (2^^) to 1.41, (2^ 2 ) times its 
original value provided that no cratering occurs, and less than this if there 
is cratering. 

Let us now consider the effects of air burst, at distances which are great 
compared to the height of burst. It is common sense to expect that at such 
distances this height itself, and all phenomena caused by it, ‘ought to have a 
small influence, and as the distance increases become entirely negligible. 

Hence one should expect, that for air burst at distances which are very large 
compared to the height of burst (as for ground burst at all distances) the 
pressure is amplified by 1.41 at most if there is no cratering and by appro- 
priately less if there is. 

Now one must observe that this picture is in disagreement with the acoustic 
theory. This is of importance since the above picture is plausible and will 
turn out to be the correct one. 

10.2-1 Acoustic Picture of Air Burst 

We will now describe the action of an air burst charge in the acoustic 
approximation. 

Independently of this approximation, the shock which emerges from the bomb 
bursting above the ground, will be spherical and remain so until it makes contact 
with the ground. Indeed, up to that time it is effectively in free space and 
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behaves accordingly. 

Whan the shock sphere hits the ground, it produces a reflected shock. In 
the acoustic approximation this reflected shock is a part of another, which is 
congruent to the first one, and behaves as if it were the blast ware cooing 
from a "virtual bomb", which is situated at the image point of the real boefe, 
reflected with respect to the ground (Figure 3). When the shock first hits 
the ground, it does so at normal incidence, and hence, by acoustic theory 
doubles the overpressure, i.e. the pressure increases over atmosphere in this 
region. Even later, when the shock sphere intersects the ground at oblique 
angles, the Ifcws of acoustics also call for a doubling of overpressure or 
pressure increase above atsnsphere in the twice shocked region at the refloating 
surface . (Figures 4 and 5). In other words: At all angles 0°, < 90®, 

the reflected overpressure is independently of o( equal to 2p (p « the incident 
overpressure). For c* s 90° the incidence is glancing and no reflection in the 
sense of acoustics dccurs (of. however, the discussion of ground burst above; 
cf. also the non-acoustic discussion of nearly glancing incidence in the 
succeeding pages). 

Let us now try to give a complete chronological account of the sequence of 
events as the blast wave expands, gets reflected, etc. Figure 6 shows the 
appearance, in the acoustic approximation, of the blast pattern as it travels 
outward from the air burst bomb. 

The essential features of the reflection phenomenon are: 

(1) Incident and reflected waves make equal angles with the ground. 

(2) the pressure* increases in the incident and reflected waves are 
equal . 

(3) The total pressure increases exerted by the blast at any point 
on the surface is twice the pressure increase to be expected at 
the same distance in the absence of the ground. 
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Figure 3 
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Spherical continuation of real blast eaves 
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Figure 4 


APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 


INCIDENT SHOCK 

PRESSURE 

p + p 

0 

REFLECTED SHOCK 

A, 

\/ 

PRESSURE ^ 

\ 

/ PRESSURE 

P 


P4-2P 

0 

\ / 

0 


CL Q 



REFLECTING WALL 


APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 



Figure 5 
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Figure 6 

Reflection of a blast warns in the acoustic limit 
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(4) The incident and reflected waves have a constant separation, 

equal to 2h; i.e. ^wice the height of burst at the zenith. At 
an angle from the zenith the separation is smaller; it tends 

to 2h cos 9 as the wave expands out from the bomb (real and 
virtual). For 9 s 90°; i.e. on the ground, the separation is 

of course 0, since the incident and reflected waves are necessarily 
In contact there. Correspondingly, the separation tends to 0 
a» 90°; i.e. as the ground is being approached. 

(5) As the wave expands, 9, the angle of incidence (and reflection) 
starts at 0° and approaches 90° in the limit of large distances* 

Figures 7 to 10 indicate the p, t dependence (overpressure versus time- 
dependence) at selected positions in the blast pattern: Figure 7 shews the 
single rise to 2p everywhere on the ground; Figures 8a, b and c, the double rise 
to p and to p + p' at a fixed height above the ground for increasing projected 
distances from the bomb; Figures 9a, b and c, the corresponding double rise at 
a fixed zenith angle and increasing distances from the bomb; Figure 10, the 
double rise anywhere at the zenith, p is the free air overpressure at the point 
under consideration; p' is the additional reflected blast overpressure which is 
describable as originating from the "virtual bomb* - since any point above the 
ground is closer to the real than the "virtual bomb" p>p', but this difference 

a f 

tends to zero (even relatively to p) as the distance from the bomb increases. 

In Figures 8a, b and c, the time interval A t between the two shock 
(pressure rises) tends to zero. In Figures 9a, b and c. At tends to 
where C is the velocity of sound. In Figure 10 A t is . 

All these sketches are strictly valid for shocks with overpressures of 
infinite duration only (Figure 11a). Actually a decay of the overpressure occurs 

behind the shock: because of the finite duration of the explosion (Figure 11b).- 

* 

Consequently Figure 7 is changed in the manner indicated in Figure 12, and 
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Figure 8 
a,b,o 
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Figure 10 
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Figure 11 
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Figure 12 
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Figures ® to 10 are affected as shown in Figure 13. 

the situation depicted in Figure 13 is of sufficient Importance ^as a 
consequence of the finite duration of the shock overpressure, to deserve one 
more comment. Depending on the height of burst it is quite possible that no 
amplification of the original (free space) pressure by the reflected shock occurs 
everywhere. 

The drop shown in Figure 13 between the pressure rises p and p> may well 
exceed the second pressure rise p'. However, at car near the reflecting surface, 
because of the (exactly or approximately) simultaneous arrival of the two shocks 
in that region, ajapl if i cation must certainly occur. 

jr. 

10,2-2 Criticism of Acoustic Picture - Formation of a Mach Stem 

Certain parts of the above acoustic description of the blast pattern from 
a charge burst above •►rigid ground are in disagreement, first of all with our 
intuition and, secondly, with the facts. 

As to the first point, we have already spoken of replacing the source, 
wall system by a real and virtual source - an explosive dipole. Now one would 
expect that from Large distances, such an explosive dipole having the combined 
mass would look like a single charge. Ibis conclusion is analogous to the result 
in electrostatics, according to which the field produced by two equal electric 
charges of the same sign is essentially the field of the total charge when the 
distance to the point of observation is great compared with the separation 
between the charges. As to the second point, in reality for shocks of finite 
strength, the situation is like the one expected intuitively in the preceding 
section, for the following reason. The permanent separation of the original 
and the reflected shocks in acoustic theory (cf. above) is clearly due to their 
having the same velocity i. : e. sound velocity. Actually shocks .of finite 
strength are faster than sound. Furthermore, the reflected shock is faster than 
the original one because it: travels through air heath# by the former, and hence. 
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has its speed increased relative to it. In fret, we know from hydrodynamics 
(cf. Section 10.3) that if one shock follows another, and is in the region of 
positive overpressures behind that shock, i.e. in the positive phase, then it 
travels faster than the first shock. Consequently, while acoustically the 
reflected shock is not able to catch up with the incident shock* if the shock 
has a finite site it will in reality catch up in regions where the positive 
phases overlap. Since the two shocks are close together near the ground and 
in contact at the ground, the positive phases certainly overlap in this region. 
The reflected shock is therefore faster than the direct shock and since they 
are getting mere and more parallel as time goes on, a merger should sooner or 

■V 

later take place here. 

We 'know in fact from the theory of oblique reflection and from experiment 
that the above is the case and that the overpressure at the fusion shook is 
about twice that at either of the two original shocks. AS the spherical shock 
expands conditions become suitable for fusion further and further from the 
ground. Consequently, the fused portion gradually rises and covers more and 
more of the shock sphere and it is possible to show that eventually the merger 
is complete and the two shocks are finally everywhere fused and form a single 
shock front. In other words, the shock due to the virtual bomb will have 
everywhere overtaken the shock of the true bomb; the two shocks will have 
merged over the sphere^ which corresponds to the double charge. Figures 14 

to — * 1 — ; 

If the charge is burst at such a height that the merger occurs at great 
altitude, say 30,000 feet, because of the variation in density of ths 
pre-shocked air neither the original nor the fused shocked front will 
be spherical. (3ee Chapter on Altitude Effect.) Since we are here 
primarily concerned with the merger near the ground we will not discuss 
this point. 

and 15 shew the appearance for a wave of finite amplitude, of the blast pattern 
as it travels outward from the air burst boob. The essential features of the 
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Figure 14 

Reflection of a spherical blast save of finite amplitude 
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Figure IS 

Growth of Ifech stem' 
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reflection phenomenon are : 

M >. 

(1) Incident and reflected waves do not intersect on the ground 
for all angles of incidence greater than a critical angle . 
Wien incident and reflected waves intersect on the ground they 
do not make equal angles with the ground. 

(2) The overpressure in the incident and reflected waves are unequal. 

(3) The total overpressure exerted by the blast at any point on the 
surface varies with the height of buret, the projected distance 
from the bodb, a nfy the blast energy released by the bomb. It 
is not obtained as in the acoustic approximation, by multiply- 
ing the fjpee air pressure at the point in question by a factor 
of two. 

(4) The incident and reflected waves have a separation at the zenith 
which^as the waves expand, at first varies little from the value 
2h and then decreases to zero as the fusion process proceeds. At 
an angle from the zenith the separation in the early stages of 
the expansion is sn&ller and becomes zero as the two waves fuse, 
forming the Mach stem. 

At distances which are large compared to the height of burst, the 
direct and reflected waves from an air burst bomb have fused and 
proceed outward as a single shock. From complete fusion on, the 
shock wave appears to have come from double the charge detonated 
on the ground. 

(5) As the wave expands, 9 , the angle of incidence starts at 0° 
and the stem becomes perpendicular to the ground in the limit of 
large distances. 

Figures 16 and IT indicate the p, t, dependence (overpressure versus time 
dependence) at selected positions in the blast pattern. Figure 16a shows the 

A <#* 

APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 


x -30 






Figure 16, a,b,c. 
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Figure 1? a,b,C 
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single rise to e< p everywhere on the ground. Figures 16b, e, show the 
dependence of o( on the angle of incidence for strong end weak shocks respec- 
tively. Figures 17a, b, and c show the double rise to p and then to p ♦ p‘ at 
a fixed height above the ground for increasing projected distances frapt. the bead). >, 
Figure 17h, b, and c also apply to the corresponding double rise at a fixed 
zenith angle and Increasing distances fro* the bomb. All these sketches are 

K i' r 

valid only for shocks of infinite duration. If the pressure has a finite dura- 
tion, Figures 16 asdl? should be modified in a manner similar to the modifications 
(Figures 11 to 13) of Figures 7 to 10. 

10.2- 3 Oblique Reflection - Deviation frOt Acoustic Behavior 

A detailed theoretical discussion of the shock wave pattern produced by the 

■* 

reflection of an expanding spherical shock wave from a rigid plane surface must 
start with a description of what happens exactly on the ground where the incident 
and reflected shocks always intersect each other, and where the reflected shock 
was originally produced. It will be recalled that according to the acoustic theory 
a reflection always produces the same increase in overpressure as the incident shock:. 

This result is independent of the angle of incidence up to, but not includ- 
ing 90°, at which it ceases to be valid. If the shock is of finite strength, 
there will be deviations from the acoustic result. Since the acoustic behavion is 
discontinuous at 90°, one would expect deviations from acoustioity t&, set in \ 
earlier than 90° for the true shock must exhibit some kind of continuous behavior 
in the neighborhood of 90°. There are, then, two factors which perturb the 
acoustic behavior: the finite strength of the shock and the obliquity, especially 
in the neighborhood of 90°. 

10.2- 1 Read-on Collision 

Let us first consider the true behavior of a finite shock for a head-on 
reflection. (See Section 10.3 for detailed calculations.) The result is usually 
stated as follows: A shock which hits on absolutely rigid wall will be 
strengthened more than the acoustic theory predicts. The acoustic theory 
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predicts a doubling of overpressure at the rigid wall. This is an incomplete 

description because there are two different ways to measure the strength of a 

shock. One measure is the absolute overpressure in the shock; i.e., the 

difference in pressure immediately before and immediately behind the shock 

front. The other measure of shock strength is given by the relative overpressure; 

i.e., the ratios of the pressures immediately before and behind the shock front 

minus 1. Now, in the limit of small overpressures, it doee not matter which 

criterion is adopted because the reflected shock is equal in strength to the 

incident shock if measured either by the pressure difference or by the pressure 
v 

ratio minus 1. For finite shocks, however, the reflected shock is stronger 
than the incident shock if measured by the pressure difference and weaker if 
measured by the pressure ratio minus 1. 

The clearest example of this is the case of a very strong shock. Here it 
is well-known (cf. Section 10.3) that for air the reflected shock has an over- 
pressure eight times greater than the original one. On a relative scale 
although the reflected pressure is eight times as great as the incident pressure, 
the incident pressure is very great compared with the initial pressure of the 

i- 

unshocked air and therefore the strength of the reflected shock as determined 
by the second criterion is very much less than the strength of the incident 
shock. To summarize, the absolute overpressure is greater and the relatlum 
overpressure is smaller. We will measure the shock strength by the abeo&ttte 
overpressure and on this basis the reflected shock is stronger than the incident 
shock. 

10.2-5 Simple Oblique Reflection , 

In the acoustic limit the incident and reflected shocks form a stable 
configuration regardless of the angle of incidence regardless of the shock 
strength. Since the two shocks have the same velocity component parallel to 
the ground as their point of intersection, 0, the angles of incidence and 
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reflection will also be equal. lhi* ia, of coarse » jaat the application of 
Snell's principle in a very simple case (Figure 18). 

let us non consider the case of a finite shock, where the reflected and 
incident shocks aay be of unequal strengths and hence have unequal velocities 
relative to the air and ground. If the intersection of the tlwo shocks is to 
remain on the ground this difference of velocities requires that the reflected 
and incident shocks stake unequal angles with the ground, if we set up the 
equations of notion (cf. Section 10.3) we find that there are as «»ny equations 
as there are variables and, in general, there are solutions. 

fe mst qualify this last statement. In the first place, when there exists 
a solution, there exists, in general, another one as well. Thus there are two 
solutions; i.e., given an incident shock the reflected shock can be in two 
positions, one less and one more inclined to the surface (figure 19). The 
reflected shock which is steeper, turns out to be faster. This is easy to see 
both qualitatively and mathematically. One, therefore, has to ask which of 
these two shooks exists in reality. It is relatively easy to see that under 
these conditions it is the less strong shock which exists because if we continue 
to decrease the strength of the incident shock and go over to the acoustic ease 
we find that the less steep shock goes over asymptotically to the same strength 
as the incident shock at the same angle, whereas the steeper shock goes over 
into a finite non-acoustic shock and becomes vertical; i.e.. as p incident 
— * 0 , ©3 ► ©]_, ©2 — * ^ • No* sines this latter does not happen, and is 

T“ 

energetically impossible without an external source of energy, we assume at 

«/* ~ 

least in the case of a weak shock that it is the less steep, weak solution 

which exists. One might assume by continuity, although this argument is not 

quite safe, that the strong reflected shock is forbidden for all incident shock 

a. be 

strengths. The experimental observations of all who have worked on this subject, 
in particular the very detailed observations of L. G. Smith, support this 
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Figure 18 

s 6 ^ acoustic case 
q / 0 ^ non-* c oust ic case 
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Figure 19 
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conclusion. If we consider this oblique shock theoretically and let the angle 
at which the incident shock strikes the wall increase from 0° towards 90®, 
the glancing value, we find that as this happens, the less steep of the two 
reflected possibilities gets steeper and steeper, the initially steeper 
solution becomes less steep. In other words., the two solutions for the re- 
flected shock move toward each other. We are satisfied that it is the less 
steep solution which is real. However, when the Incident shock has reached a 
sufficient obliquity the two solutions for the reflected shocks merge; i.e., 
become identical, and beyond this there is no solution;^) so there is an 


For shocks of any reasonable strength this obliquity is far from 90° 

(cf. Chapter 3 ; e.g., for a shock strength of about an atmosphere this 
angle of incidence is about 50° and even for a shock strength of 0,1 
atmosphere this extreme angle is around 80°. 

extreme angle below which there exist two solutions for the reflected shock, 
an angle at which there is just one solution and beyond which there is no solution 
10.2-6 The Critical Angle - Irregular Reflection 

there is then, an extreme angle below which there are two solutions. We 
choose the lower solution for tolerably good theoretical reasons which are very 
well confirmed by experiment, but beyond this extreme angle there simply is no 
reflection of this type, there must, however, be a reflection of some kind. 

We call the reflection of the first kind, the reflection which really is the 
extrapolation of the acoustic case, regular reflection, and all kinds of other 
reflection, irregular. From the above, one can see that the region of regular 
reflection is limited to angles of incidence below a certain critical angle. 

Now the variations with shock pressure of the limiting angle for regular 
reflection are of some interest (Figure 20). It is quite clear that as the 
shock strength is decreased, the acoustic limit is approached and, therefore, 
the limiting angle must be nearer and nearer to 90°, For reasons previously 
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Figure EO 

Critical angle at which regular reflection 
ceases to be possible, as a function of the 
overpressure in the Incident shbck 

Ideal gas ;( m 1.40, 

Nonas! pressure * 14.? psi 
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mentioned, as the angle of incidence becomes close to 90°, deviations from 
acoustic behavior should be expected no matter how weak the incident shock is. 
Theory as well as experiment show that such deviations do occur. As will be 

l 1 

shown later (cf. Section 10.3) the critical angle converges to 90° rather slowly, 
as the square rcot of a shock strength. More interesting is the fact that as 
the shock strength increases, the critical angle decreases from 90°l For a few 
atmospheres overpressure it gets into the neighborhood of 4 0°. For one and 
one-half atmospheres overpressure it reaches 40°. After this it does a peculiar 
thing. It drops a little below 40° to something like 39° which it reaches at 
6 atmospheres overpressure and then it rises again to 40°, which value it retains 
for infinitely strong shocks. This means that for shocks of 1.5 atmospheres or 
greater, the critical angle has already practically reached its limiting value 
of 40°. For half an atmosphere overpressure, it is 50° or 60°. Since moat 
blast damage by large bombs is baaed on controlled pressure criteria and likely 
to occur between 3 and 6 pounds per square inch overpressure, that i», between 
a quarter and one-half atmosphere, we must expect regular reflection to become 
impossible in the neighborhood of 50° to 60°. 

Another interesting characteristic of regular reflection is the variation 
of the absolute overpressure on the surface as a function of the angle of in- 
cidence and the strength of the incident shock. Figure 21 gives the relative 

overpressure as a function of the angle of incidence for acoustic, weak and 

> ■ 

strong shocks. 

10.2- 7 Mach Reflection 

Now we must ask what will happen* if we go past the extreme angle of 
regular reflection. There is a very simple argument which uses the analogy of 
this phenomenon with the collision of a blast with a wedge. 

10.2- 6 Collision of a Super sonic Fleet with a Wedge 
A plane supersonic flow is incident on a wedge of semi- angle 9 W . 
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Flyura 21a 

Relative overpressure on (surface due to reflection 
versus angle of incidence 
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Figure 21b 

Relative overpressure on surface due to reflection 
versus angle of incidence 


Ideal gas Y* 1.4-0 

i'-.iL , 4 - 

fi ' * 

where p * pressure in front of incident shock; 
p 0 * pressure behind incident shock 
p* * pressure in front of reflected shock 


I 


For acoustic oas« £ £ = 
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A — overpressure in Incident shock 
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flout fcolllde* with: the wedge it ia deflected through an aggle^t^ + ) by 

a shock SOS becoming parallel to the plane aide of the wedge. The conditions 
are constant throughout the regions A and B which are separated by the shock 
discontinuity SOS. It is found experimentally, and could be shown by theoretical 
consideration of an oblique shock, that depending on the wedge angle 0^, 

T£5 — — 

Thy lor and MaeColl, "The Mechanics of Compressible Fluids" in W. F. Durand 

"Aerodynamic lheory", Guggenheim Fund, 1934, Vol. Ill, Section H. 

there is a critical value of the supersonic material velocity for which no shock 
exists which is capable of deflecting the stream so that it becomes parallel to 
the wedge wall. Above this critical value of material velocity there are in 
general two theoretical solutions ^ for each wedge angle. Only the solution 
with the smaller value of 0^ actually occurs. At the critical value there is 
just one ' solution and below it there are no solutions of the type pictured in 
Figure 22. Instead, experiment shows a detachment of the shock wave from the 
top of the wedge such as pictured in Figure 23. This is the so-called "detached 
headwave". This means that there is no solution of the type pictured in Figure 22, 
that considering a given value of supersonicity, there exists an upper limit on 
the ftftgle ef incidence (Qft - © w ) for which the flow can be rendered parallel 
to the wedge face immediately behind the shock. Beyond this angle a new type 
of phenomenon is in evidence. This is the detachment of the shock. It is of 
particular interest because it means that signals are sent back from the tip 
at 0 into the region bounded by a now curved shock* there is a propagation of 
signals back against the impinging stream. 

The connection between this phenomenon and the reflection of a blast wave 
from a plane surface can be shown in a simple way. In Figure 24, the incident 
shock impinges on the wall giving rise to a reflected shock. If we adopt the 
frame of reference which reduces the motion of the point P and the two shocks 
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Figure 22 
Flow pest wedge 
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Figure 22* 

3»personic fleer pest sedge 
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Figure 23 

Detached headway* 
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to rest then the amterial will appear t© flow through the shocks from left to 

right as shown. The incident shock causes the material originally flowing 

% 

parallel to the wall to be deflected toward the wall in the direction of X. 

The reflected shock renders the flow again parallel to the wall. If we identify 
the flow in the region between the two shocks with the supersonic stream of 
Figure 22, the reflected shock with SO and the wall to the right of P with the 
wedge wall in the same figure, the analogy is complete (Figure 22»). The un- 
applicability of Figure 22 under the certain conditions discussed above is 
reflected in the failure of Figure 24 under corresponding conditions; i.e*, 
there exists a critical angle of reflection (end hence of incidence) for a given 
strength shock. Beyond these critical conditions we expect from the analogy 
that a signal propagates back from the reflected shock into the region between 
the shocks and causes a fusion of the incident and reflected shocks starting 
in the neighborhood of P. 

> 

This analog; therefore, shows that the reflected shock must be expected 

«■ 

to overtake the Incident shock whenever 0 > © c . This process of overtaking 
originates at the wall and gradually spreads into the volume of the gas. As 
it spreads the two shocks merge and form a single shock for a certain distance 

'%r-' 

from the wall, beyond which they are separate. In the case of irregular 
reflection, then the two shocks will no longer look like a V btrt like a T 
standing on the wall. (See Figures 25a and b.) 

In Mach reflection it ia as if reflection were no longer caused directly 
by the wall but rather by a cushion of air resting on the wall. 1km that this 
species of irregular reflection really occurs, it is a seventy year old experi- 
mental observation of EmstHtach^ 5 ^ , after whom it was named. It has been 

TTJ — ' 

Paper of E. Mech with various collaborators appeared in the Tienna Academy 
"Sitzungsberichte" fol. 72 to 92 (1875-1889), cf. in particular 

Vol. 78 (1878) page 819. 
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Figure 26,a»b> 

25a Regular reflection 
25b Mach reflection 
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explored in great detail by hydrodynamics! research done daring this war. 


TS1 

J. van Neumann, “Oblique Reflection of Shocks", Bureau of Ordnance E.R.R. 
No. 12 (1943). 


H. Polachek and R. J. Seeger, “Interaction of Shock Waves in water like 
Substances", Bureau of Ordnance BJt.B. No. 14 (1944). 


P. C. Keenan and R. J. Seeger, "Analysis of Data an Shook Intersections, 
Progress Report No. 1”, Bureau of Ordnance E.R.R. No. 15 (1944). 


H. Polachek and R. J. Seeger, NavCrd No. 88-46, “Analysis of Data in c -* s. 
Shock Intersections, Progress Report No. 2" (1946). 


KavOrd Report No. 74-46, "fliclwiical Conference on Optical Phenomena is 
Supersonic Plow" (1945) C^enf identigl) . 



Lincoln Q. Sadth. "Photographic Investigation ’of the Reflection of Plain 
Shocks in Air" ' (1945) (Final Report] v OSRD Report No. 6271. 

— 7 —— — — — — — — ■ * — — " ■- 

The qualitative picture of Mach reflection is quite single for the case 

i . -"If 

Of a plane shock incident on a wall. Here the incident shock wakes a constant 
angle with the wall. The situation is more complicated in the case of blast 
produced by a bomb. First of all, the angle of incidence changes as the blast 
wave proceeds outward from the bomb. In addition, the sphericity of the blast 
wave makes a quantitative difference where Mach reflection occurs, while it 
introduces no additional features in the case of regular reflection. This is 

j 

so because regular reflection takes place entirely in the neighborhood of a 
single point of the wall and, therefore, only local conditions at that point 
enter. The definable si*e Y type reflection, on the other hand, extend* over 
a finite area and grows up on the shock. Therefore, the properties of the shock 
in the large area now become relevant. As indicated earlier, the reflected 
shock merges with a continuously increasing fraction of the incident shock, and 
eventually the incident and reflected shock may even coincide completely . Most 
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of the theoretical and experimental work on the phenomenon, but not all of it, 
was done in the case which is somewhat simpler; i.e., a plane wave incident 
on a plane wall. It is found that, in general, the dimensions of the Y in 
this case are not constant, but that they grow with time. In other words, while 
the regular reflection, which produces a V, is stationary, that ie, the V never 
changes, this new kind of reflection is not stationary, i.e., the Y grows as 
time gees an. Of coarse one must admit that the I contains a length, and hence 
a size can be attributed to it, in a manner which is not possible for the V. 
the length associated with the I is tiie length of tbs stem; the V has no stem 
and definable size. Since all experimentation, as well as theoretical consider- 
ations, show that in the plane case the stem of the I grows proportionally in 
time, the Mach effect must have a well-defined beginning. The length of the Y 
stem at any moment defines the duration which the phenomenon must have had from 
the time of its inception to the time of observation. 

10.2-9 Difficulties with Irregular Reflection 

As was shown at a certain obliquity 9 C , regular reflection ceases. All 
forms of reflection occurring after this, i.e. for 0> C^, are by definition, 
irregular. It is believed that irregular reflection at its very inception, i.e. 
for the 0 immediately following , belongs to the type described above as Mach 
reflsctJbon, although if one goes into the minutiae, there is an interval of a 
few degrees where one might conceivably have doubts, these doubts as to the 
nature of early irregular reflection are based on experiments by £. G. Smith. ^ 

E 

vn 

Lincoln G. Smith, "Photographic Investigation of the Reflection of Plain 
Shocks in Air" (1945) (Pinal Report) OSRD Report Ho. 6271. 

We would expect that the small angle between the trajectory of the triple point 
and the wall should be zero exactly where regular reflection ceases. According 
to Smith's experiments the triple point is hot visible for about 2° after 
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regular reflectien ceases. If the angle between the triple point trajectory 
and the wall la measured In the region of 2° to 15° past the critical angle of 
incidence and then extrapolate back to zero, it hits zero at about 1° ft* on the 
extreme angle. Ifcether this effect ie real we cannot be certain as yet. It 
■ay be that a higher resolution will settle this point. 

A further co«plication is caused by the fact that the reflection phenomena 
must necessarily have a beginning at a definite point, where the shock first 
■sets the oblique wall. This point is clearly a corner which can belong to 

s 

any one of several types. Figures 26a, b, and c are some examples of such 
corners. (It will be noted that 26b and c represent alternatives, which are 
equivalent since viscosity and wall friction are disregarded . ) Both theory and 
experiment show that this comer must be the source of a disturbance behind the 

• X . 

reflected shock R. Smith's schlieren photographs indicated that this disturbance 
always a rarefaction wave. Ihls rarefaction's edge shows up as a curved wave 
ft* terminating (an the back aide) region between R and the 

wail. It makes contact with S at S and beyond S it has "eaten into" R, and 

thereby replaced the straight shock R by a weaker curved shock R*. It will be 

* , 

noted that ft' does not catch up with R along the wall. This is so because the 
above mentioned homogeneous flow immediately behind R is supersonic, and R' is 
4 rarefaction, and hence precisely sonic. See Figure 27. 

• As the obliquity of I, 0 , increases, the flow behind R becomes increasingly 
nearly sonic, and correspondingly R» moves close to R, and S' close to the wall. 

Consequently the straight piece R shortens. At a certain angle 0 Cvt the flow 

* 

becomes exactly sonic, S reaches the wall and the entire reflected shock becomes 
curved, i.e. R* replaces R in its entirety. Thus Qt u ha3 been computed, it 
lies very cloee to Qfe but it is definitely less. (9 C - © cu varies with the 
shock strength, in air it is never more than about .6°.) From the point of 
view of the empirical evidence one cannot even be quite certain, whether the 


APPROVED FOR PUBLIC RELEASE 



APPROVED FOR PUBLIC RELEASE 


X- 52 


i 


Figure 26,a ( b,c 

26b Obtuse angle obstacle 
26c Acute angle obstacle 
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Figur# 27 


i 


APPROVED FOR PUBLIC RELEASE 



APPROVED FOR PUBLIC RELEASE 





APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 


X -54 

regular reflection still exists beyond , i.e. whether it ceases at 0 C or 
at • 

In the plane cases studied, it was found both theoretically and experi- 
mentally that the shock configuration remains similar to itself in time, i.e. 
it can be described in terms of two variables, r/t, and 8 , where r and 8 are 
polar coordinates with the origin at the corner, and t is the time of travel of 
the incident shock from the corner to the position given by r , 8 . We repeat, 
the configuration remains similar to itself as it proceeds, the triple point 
travelling a linear path (Figure 28). 

In the region of regular reflection it is clear that the angle of the 
triple point trajectory is 0®. We know from Smith's results that with the onset 
of irregular reflection the angle y? made by the trajectoiy with the plane 
surface (Figure 29) is small but becomes greater as the critical angle is 
exceeded by larger amounts. Actually, it is very small for a considerable 
angular interval. For overpressures of about 1/4 atmosphere, the range in which 
the most detailed observations are available, the Mach effect begins when the 
angle between shock and wall is 56°. Even at 66° incidence the angle of the 
triple point trajectory with the wall is only 1 or 2°. 

10.2-10 Spherical Mach Effect * 

* 

In the case of a spherical wave, there is no question of edge effects but 
additional eompli cations arise from the fact that the angle of incidence changes 
and the shock weakens as it proceeds outward from the center. The increase in 
the angle of incidence of the expanding spherical shock is, therefore, expected 
to be accompanied by an increase In the rate of rise of the triple point, as 
already indicated in Pigure 14. The rise of the stem of the Y takes place in 
a manner for ilhich there is a good qualitative description. Halverson and Tbub 
have made extensive studies of the experimental data on the Mach effect for 
apherical waves. Despite the lack of a satisfactory theory, the use of available 

J ' 
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Figure 28 

Linear travel of triple point 
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Figure 29 

Angle made by triple point trajectory with pinna surface 
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experimental data togetbar. with tha »***' scaling law enables us to predict the 
sain features of the behaviour of a reflected spherical blast wave. 

Returning again to the case of a plane wave, we note that certain diffi- 
culties appear in the study of irregular reflection as the angle of incidence 
approaches 90°. As the angle of incidence gets close to 90° the appearance of 
the reflection changes considerably. Specifically, the contribution of the 
reflected shock begins to becoae less and less, and in the neighborhood of the 
triple point, tits back branch of the I, i.e. the reflected shock, gets weaker 
and ‘weaker. In other words, a situation develops where it gets progressively 
mord difficult to tell the forward branch of the upper part of the T from its 

f 

stem and to observe the backward branch of the upper part of the T. 

If one is completely phenomenological, if one talks only about what one 
sees and not what one expects theoretically, then one must admit that in the 
Mach effect, after 80° the back side of a I has not been observed. It looks 
as though the incident shock simply makes a turn and gets deflected away from 
the wall without the benefit of the reflected shock. Further, from the triple 
point the reflected shock is observable but it does not reach the other shocjpn , 

,io fora a T. This may be so, not because it is absent in this region, but | 

'£ 

because it gets too weak to be observed with the experimental technique employed. 
However, we really do not know whether late reflection is of the Mach type or 
has soma other fora. 

In any case, this late form of irregular reflection is the one which 
corresponds for the spherical wave produced by the bomb to the double point 
charge, and which, in the end, yields a pressure increase by a factor of 2^3 
or 2 ^/^, depending on the pressure distance curve in the region of interest. 

This concludes our qualitative description of the reflection of plane and 
spherical shock waves from a rigid plane surface. 
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We no w proceed to a more detailed discussion of the present status of 
the theory of regular and Mach reflection. We Hill refer frequently to 
experimental results. 

10.3 THEORY OF REFLECTION 

In this section me treat in some analytical detail certain relevant topics 
in the theory of reflection, some of which have been mentioned in the general 
discussion of Section 10.2. We assume step shocks in non-viscous gases 
throughout, although such stringent restrictions are not necessary to all ensuing 
discussion. First, we discuss the one dimensional case: a step shock nornmlly 
incident on a wall, and the phenomenon of "eatchup" which occurs when the 
positive pressure region of two shocks overlap. Next, we consider questions 
related to the regular reflection of a plane shock such%s the pressure multi- 
plication as a function of the angle of incidence and shock strength, and the 
critical angle beyond which a two-shock solution no longer exists. Finally, 
there is a discussion of the irregular reflection of a plane wave and attendant 
difficulties. 

10 , 3-1 Shocks In One Dimension 

We will first considor a single plane shock, starting with the classical 
Bankine-Hugoniot formulae. Consider such a shock S (Figure 30). It separates 

T® 1 '* 

For a derivation of these formulae, see Vol. XIII, Chapter 3 on Shock 

Wives, by K. Fuchs. 

two domains M 0 and M, in which the physical characteristics of the substance 
pressure and specific volume have the values p 0 , v 0 , and p, v, respectively. 

The corresponding sound velocities are c 0 and c, respectively. It is helpful 
to tie the frame of reference to the shock front S. Then the substance of U Q 
flows inio the shock with the shock velocity V OJ while the substance M flows 
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Figure 30 
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out of it with the velocity V. the shock S thus ceases e decrease in velocity* 


U a V 6 - V 


and its compression ratio is J 5 


the formulae of Rankine and Hugoniot express the c erase rvation of mss* 
momentum and energy . the first part of this system correlates the inner 
properties of the substance on both sides of the shock, i.e. p Q , v Q and p, v: 

E' E . = i(P + P.)U-^ (i) 

Where E is the energy density of the substance in the state M, and E 0 is the 
energy density in the state U 0 ; ^ 


E--e:(f>,v) , E,= E.fo,%) 


( 2 ) 


The second part of the system expresses the velocities V» -V, U in 

* T 

terns of p 0> v Q and p, v» 

^/v (£ is the sign of V* - V ) O) 

and hence the difference in material velocity, relative to the shock front is 


U »V 0 ~V - + 

Equations (1), (3) and (4) imply that 


(4) 


f>l;i.e. implies lr< ir„ andvV»>V>oj (5) 

In Figure 25 this implies a shock "facing* right / 

|<jji.e. }><)». implies \r)vi, and VX\£<0\ (6) 

In Figure 25 this implies a shock "facing" left J 
For p— ► p 0 the shock becomes a sound wave. (1) shows that in this limit 
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ci £ > £ 

Jo ~ . to the other hand since (identically) at constant entropy, 

an infinitesimal shock leaves the entropy asymptotically constant. Now (3) 
gives V e ~V~ : »V|F * *^ 1 constant entropy, i.e. adiabatically the 

shock velocity becomes sound velocity. Therefore, with respect to an infinite- 
simally weak shock the flow on both sides becomes asymptotically sonic. 

For a shock of .finite strength the following rule is generally valid: 

(A) The shock velocity is supersonic with respect to the lower pressure side 
and subsonic with respect to the high pressure side. 

For (5) these two sides ard Uq, M, respectively. For (6) they are M, M 0 , 
respectively. So we have: 


5>l 

implies | V* j > C* , J J 4. C 

(7) 

S< 1 

implies | V 0 | < ^ |V.| ^ C- 

(a) 




lhis rule can be demonstrated quite readily for an ideal gas but it is 
also generally valid for a qy substance with a p, v characteristic that is concave 
upward. ^ 

to 

H. A. Be the , "The Theory of Shock Waves for an Arbitrary Equation of 
State," OSRD Report No. 545. 


Consider a shock moving to the right. Restating (3) and (5) 




p- P» 
v-.-v 


But 


]] ZJL) 


c --\ 

1 d / W (J0na tant 

y 

v-o- ^ 


entropy 


I constant 
entropy 


Far an ideal gas the adiabatic law and tile equation of state, and the 
internal energy density are respectively 


(a) 


jjv* ~ constant , (b) ,pv * RT , (c)£s (9) 
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So C * 


i*p * c* - ,/ ^ v 


that ore with to prove is that 

c>V>c. 

Substituting the expressions forA4 1 C,C 0 » this Inequality becomes 

vtr > v.yEE- > 


( 10 ) 


or 




where 


i.- 


ui) 


An additional relation between ^ir 0 > ^ and }f is required so that the 
correctness of this inequality can be demonstrated. Substituting (9c) into (1) 
we obtain 


v _ 

V. 


(ir-0 5 

(tf-H)f 4 C if- 0 


(12) 


Combining (11) and (12) and simplifying we get in each case the same inequality, 

aw. 

namely 


• t > l 


(13) 


Dlls condition I was already required by the energy density condition (9c). 

This proves the statement (A) for an ideal gas. 

10.3-2 Shock Catchup 

Having demonstrated (10) we are now in a position to show under what 
conditions two successive shocks will merge. Consider two step shocks (Figure 31). 
Ey (10) shock I is supersonic with respect to undisturbed medium ® , and is 


subsonic with respect to the once-shocked medium (1), i.e. C t > V, > C„ 
Similarly, shock II is supersonic with respect to once-shocked medium (1) and 
subsonic with respect to twice-shocked medium (2) i.e. C t >V„ >C ( 

Combining these two statements W 0 > C, >\^ or v. > v. 
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Figure 31 
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anti Shock II will overtake Shock I. This is true as long as | and hence 
catchup' will certainly occur whenever the positive phases of two shock waves 
overlap .(Figure 32). 

A more detailed investigation would show that catchup would occur even if 
the second shock were in the early part of the rarefaction region of the first 
shock,' but we will not go into this question here. 

10.3-3 Normal Reflection by a Rigid Wall (Ideal gas) 

Ihe Figures 33 have three domains A, B, and C which are defined in terms 
of the number of times that a shock has passed through them. The pressure, 

a . 

specific volume and sound velocity in these regions are p, v, cj p Q , v e , c 0 j 
p'Vv', c', respectively. 

. v. 

Since reflection increases the overpressure 

P < P 0 < P' 

or § < I < 5 ' 04). 

where * §'* 

In the present setup there occur no movements parallel to the rigid wall. 
Hence the substance of the domains A and C which cannot move normally to the 
wall either is necessarily at rest. The substance of the domain B has therefore 
the velocity U from Figure 33a and the velocity U' from Figure 33b. The 
condition relating the strengths of Shocks I and R is therefore 


U 3 U» 

Now for an ideal gas characterized by (9&, b, c) we recall that 

. - S ♦ (W) 

/v * ' (#+0 $ +( a -0 


Solving for 0, U' from (2) and (4) we find that 

u _ , 

'c. * + -y 2-r [' 


(14') 


( 12 ) 


(15a) 
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Figure 32 
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Figure 38 e,b 
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' c * " Vi^oTTTT^T <15b 

Ihe minus sign in (15b) appears because of the opposite directions of travel 
of the incident and reflected shocks. Equating (15a) and (15b) as required 
by (14) and solving for £ 1 in terms of £ 


6*-0 ~ ( | 
(**0 5 * ( v -0 



consequently 0 < J-J < l corresponds to ) < $* < 


5 = %. 

it - 1 
v + i 


(16) * 


For air (at moderate temperatures and densities) JT * 1.4 and *8. these 

T* I 

results deserve a brief discussion because they contain the first quantitative 


indications about the “unaeouetic" effects in shock reflection: that is, the 
deviations which shocks of finite strength present from "acoustic* 1 laws, whi ch 
hold asymptotically for shocks of infinitesimal strength. 

From (16) we infer easily that 


< - s < s’- 1 , 


(17) 


This means (ss was stated without proof in a previous part of this chapter) that 

< K-N * % > (18) 

The reflected shock is stronger than the incident shock if strength is measured 
by the absolute compression (pressure difference), but it is weaker if the 

strength is measured by the relative compression (pressure ratio). 

> 

In the limit of infinitesimal shocks, the acoustic limit (16) can be 
manipulated to show thatji-)>~)>-j» # ,^j^', i.e. that both absolute and relative 
criteria give the same result in the. acoustic limit. 

Thus the "unacoustlc" theory is so far just a plausible extension of the 
acoustic theory with little individuality of its own. We shall see how radically 
thi? changes when oblique reflection is considered. 
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10.3-4 Regular Oblique Inflection 

We will now consider the case of the oblique reflection of a shock of 
finite strength (Figure 34). 

A single plane shock wave, I, 10 incident on the wall, the included angle 
being U . It produces a secondary shock wave, the reflected wave, R, which 

I 

includes the angle oi with the wall. We will derive relationships between 
the strength of I (measured by its compression ratio), the value of oC , the 
strength of R and the value of ot . 

The shock waves I and R and the wall define three domains in the substance 
the unshocked material ahead of I; the once shocked material between I and R; 
the twice shocked material behind R. We denote these domains by A, B, C, 
respectively. 

The ordinary frame of reference has the unshocked material A at rest, 
while the waves I and R move with shock velocities u and U* respectively normal 
to themselves. A more advantageous frame of reference is one which is tied to 
the reflection point, P, where I, R and the wall meet. Here the substance A 
comes in with a velocity £? , parallel to the wall, while the substance B , 
leaves with a velocity g' , parallel to the wall. From this point of view, the 
first shock causes flow which is initially parallel to the wall to be deflected 

toward the wall and the second shock renders this oblique flow again parallel. 

m 

In this frame of reference, the shocks I and R and the reflection point P are 
at rest. 

Clearly the pressures in the three regions increase with the number of 
times the regions received a compressive shock, so that j^< )a 6 < p where 

I ' 

|?,)3, ,jp are the pressures in regions A, B, aad C respectively. Or: 

— 1 

5<l< 5 
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Figure 54 
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The velocity components in A * C/yunits of the eound velocity in B, 
normal to I and B: are |jg| sin ct * T > |jg'| sin oi' * ■£* i.e. 

r 


1*1 ■ 


Ua\°c 


» I H I 


W ot 


( 20 ) 


In our present units I, R modify the velocities of the substance which 
crosses the^by the amounts 


Ur s 


lul 

Oo 


„'.M 

Co 


respectively. 

If we denote the components of the velocity X in B by i and y and then 
calculate these compsnents by passing into I through I and then through R, 
we obtain: 

Through I: 


* * nsnoc , 


Through R: 


X a 


V 

JA* 

T 1 


•Cjr COS at 


(O' at,* » y » + ur' C^O-oL* 


Combining theee relations 


T + <*T sin 2 pC 


Sim. oc 


» i . a. .t 

r + t*s sm eC 


( 21 ) 


— Ur can oC 


aj ’ con oc 


( 22 ) 


where by (15) far an ideal gas 

U z(i-t) i V 

yiVfff+l) J Vfr- »>T c 

and from (3) end (12) for an ideal gas 

t . v _ kzhx± ii±il 


2 I ) 




c « ‘V^vrcx+os+^'OT 


c e ^2t[(t+o $’+ (Y-0]‘ 


(23) 


(24) 


A more convenient form of (21) is (Attained if we introduce a new variable (r 
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where g- mcJ+T: i.®. by the definitions and by (3) O' T 

<T - cxm.^oC _ <r' - <jJ ' ur*, X oC 


SA/W OC 


Swi oC 


( 21 ’) 


where for an ideal gas 


1 fiWQs*(y-0 > 

^ y z * V xr 


(25) 


Equations (21') and (22) in conjunction with (23) and (24) relate the 
shock strengths 1* , ^ , to the angles of incidence and reflection «L , oc 
so that given j£ , ec , the quantities j ' , <*.» can be obtained. 

* k 

Eliminating ' from (21') and (22)^ using (23) we find an expression for 
>'8r' in terms of <r and oC . 

(y+i)«r 4 » cr' 4 - ^ z ("e+if <r 4 + 4 * <r*(k-i)<W*oi. + r* 

+ { (itf4l)V 4 -4(«r-l)WV + [(^l)<r%l]^UV} -O (26) 


The condition (22) in conjunction with (25) yields the relation for oc' in 
terms of <y • , c T and o L . 


Coil, eC 


/r -«r 


Cot* oo 


(27) 


There are, in general, two solutions for given jr , oc : one has 

greater values for ^ ^ , the other lower values for jr • , 
value of the shock strength £ , however, there is a maximum value of (i.e. 
*-ext. ) at which these two solutions coalesce and beyond which there is no two 
shock solutions of the type predicated. For reasons cited in Section 10.2, 
the less steep solution is the one considered physically admissible. 

We will now discuss briefly the question of the existence of two solutions 
and the coalescence of these two solutions as the critical angle of incidence 
o6*xt. is approached. Rewriting (26) 

X 1 “ Bx ■* C = O (28) 


For eao£i 
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.ter.: 

6=1 + eteSt + -k rfcJlfiiaf cfttV 

fr+0 <r*- 9* a L fr+O j 

x , , \ r(» u 0rSa.1 l .te tj 

°» <+ ifer-J ^ * 

> >>'>VW 


<r > 


C = 1-4 


<r x - 


(tr+t)* 4 

(* + i) 5 + (*-») 


i * 


X.<r'\ f } (29) 

New if ^ ^ | ; < two solutions exist ( corresponds to regular reflection). 

r one solution exists (corresponds to critical angle). 

I 

> no real solution exists (30) 

These three conditions can be met. To see this, since C and B are continuous 
functions, we need only show that the conditions are net at least at two 
points. r 

At oC - 90°, B r 1, C - 1, so that > 1 and for 

<90°, Ba^[-(7Vo -J ** * > c ‘[- fr+ft "J * 

4 


oC <■ 


so that 


4C ~ 


B i. ~ cCrw»-«C 

77+7) ' 


< 


and the three criteria (30) can be met. 

Using the values of cr 1 and the assumed value of or and «. we can find 
the corresponding values of *c' by means of (27). 
cC ext obtained from * in (30) is given by the following equation 


i6*V<rM) 4 ^ 3 - fiefr +!)*■(* ^ 

-[(«.,) «r‘ + *] 4 =tt 
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where y * sin^ot 

ext 

These equations have been investigated in numerical detail by H. Polachek 
and R. J. Seeger.^^ Some of their results are embodied in Section 10.2, 

m 

Regular Reflection of Shocks in Ideal Gases (AM-524) Feb. 12, 1944. 

Figures 20 and 21. 

It is. v of interest to see how the critical angle approaches 90° as the 
shock becomes sonic. 

f * * i -a , a << i 

and by (25) <r l -l = - & . Inserting in (31) and neglecting 

terms in 

y 9 AA-n, cc ext - I - — of 
90-oc ext ' 

a 

The critical angle converges to 90° as the shock overpressure approaches 
zero, the approach to 90° varying as the square root of the shock overpressure. 
10.3-5 Mach Reflection 

Let us now turn to the case in which regular reflection ceases to be 
possible . 

When od increases beyond od,^. We wish to investigate the case in which 
regular reflection continues as long as it is possible, i.e., up to od = oC e>ti , , 
and we shall try to determine what happens beyond this point. 

That regular reflection is impossible, means that the oblique flow in 
region B (Figure 34) cannot be deflected by a standing oblique shock R so that 
it becomes parallel to the wall. This corresponds to the propagation of signals 
from the shock R, back against the flow X into the region B, a fact well 
substantiated by the experimentally observed fusion of the incident and reflected 
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shocks I and R (cf. 3ection 10.2). Another way of describing the critical angle 
is to say that it is the greatest angle at which the components of U ' and O 
orthogonally to the nail can be made equal, i.e. 

O' cos oc.' * U cos ot (32) 

Beyond this angle, the compensation of the noreal components of the floe 
velocity is no longer possible by such a simple, stationary process. 'Hie 
phenomenon of irregular reflection which is thus produced will not be atationary - 
even in the frame of reference in which P (and I) is at rest. The wave of 
impact of X on the wall will propagate from P in all directions and thus change 
in "size" at all times and never reach a final equilibrium state. However, 
while the shock configuration shows a continuous change in size, its shape will 
be permanent (cf. the comments connected with Figure 28 in Section 10.2). Its 
size is proportional to the time t that has lapsed since its formation: For 
t « 0 it begins concentrated into a single point. 

Since the shocks R and R are advancing into two different media A and B, 
a discontinuous change of direction of the RU front mast be expected at T* (Figure 
That is, T is really the contact point for three different shocks - I, R, U. 

The point p 0 is the so-called "center of similitude" from which the phenomenon 
originates (the "comer" in the sense of- the discussion in Section 10.2). 

Now consider a frame of reference in which I and T are at rest (Figure 36). 
The flow £ * of the substance in A crosses the shock system I, R, U, to 
finally reach the domain C. However, this process operates in two different 
ways even in the immediate neighborhood of T: the substance in A above the 
line t i crosses into C through two shocks I end R; while- the substance in A 
below XI crosses into C through one shock M. Since we assume no shocks beyond 
R, U (in C), both processes must compress the substance to the same pressure. 

But this compression occurs in the "upper" half in two stages (I, R) and in 
the "lower" half in one (M) . The former process is less irreversible than the 
latter - essentially because it is less abrupt. Hence the substance which 
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Figure 35 
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Figure 56 


APPROVED FOR PUBLIC RELEASE 



APPROVED FOR PUBLIC RELEASE 


* - 7T 

' t- * 

crossed "above* T may be expected to have a lover entropy than that which 
crossed *below* T. Thus we have in C near T two flows: both with sane pressure 
bat the "upper* flow has lower entropy. Since the two flows have the same 
pressure, but the "upper* flow has the lower entropy, therefore it joust also 
have the lower temperature, i.e. inner energy, and the higher density. Again, 
since it has the lower irmer energy, it follows from Bernouilli's principle 
that it must have (in this frame of reference) higher kinetic energy, i.e. 
velocity. Therefore these two flows mist be gliding past each other along a 
dividing line D issuing from T. So D is a slipstream, and also a discontinuity- 

line for entropy, temperature anddensity, but not for pressure. We therefore 

* ';£* """ ' 

have a pattern consisting of 

in T: the three shocks 1, H, 

Hi 

That D is a slipstream and nit a shock is indicated by the experimentally 
observed large angles between B and D and the 're^uiremsnt that the flow into a 
shock be supersonic With respect to the medium into which the shock travels 
(cf. rule A). The flow is subsonic behind U and it Is inconceivable that it 
can become rapid enough before reaching 0 that Its component across 0 should be 
supersonic* 

Since the existence of the extra shock, U, and the irregolar reflection 
its presence implied was first recognised by Mach, it seems appropriate to give 
his name to this typo reflection, fe shall therefore designate irregular 
reflection according to the scheme of Figure 37 as Bach reflection. 

10.3-6 Threw Shock Solutions 

On the basis of tbs above it would seen indicated that a theoretical 
investigation of Mach reflection should be undertaken. For regular reflection 
Figure 34 provided the basis for theoretical discussion . For Mach reflection 
Figure 37 ms y seem to provide the basis in the same sense, and therefore we 
hgvs to eee what possibilities of a theoretical treatment It offers. 
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M* ijfid the slipstream D (Figure 37). 
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Figure 37 
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The theoretical treatment of jh^ular reflection was greatly simplified 
by the fact that conditions in each one Of the three domains A, B, C of 
Figure 34 were constant. We needed only to discuss the purely algebraical 
connections between them by applying the Rankine-Hugoai ofe conditions to the 
straight shocks I f R, which separate them. 

In Figure 37, conditions in the domain C u are certainly not constant and 
the shock R is certainly curved. 

Probably the same is true for the domain Cj , the shock U and the slip- 
stream D. Thus we must resort to the differential equations of compressible 
fluid dynamics, aggravated by the appearance of curved shocks, which introduce 
varying amounts of irreversibility* ^ 

It would, however, be pointless to attack these difficult partial differen- 
tial equations, without having first ascertained She nature of the conditions 
at the boundaries of the areas in which they apply. These areas are C u and c^. 
The boundaries are the lines R, M, D, and the point T at which they all meet 
is one of particular interest. 

The Rankine-Hugoniot conditions take care of the situation along S and U. 

For the slipstream D we must require that it be a streamline in C u as well as 

in C £ and that at each point the pressure on its C u side be the same as on its 

Cg aide, but we must also reconcile all these requirements at the point T, 

where R, M, D meet. Therefore a discussion of the conditions at*T is a , 

• 

necessary preliminary of any theoretical treatment of Figure 37. 

The discussion of conditions in the immediate neighborhood of T is best 
done In the frame of reference where T (and I) is at rest. This is the scheme 
of Figure 37 except that we may replace the lines I, R, M, & by their tangents 
at T - i.e. we may assume them to be straight. (Figure 3&). 

The following remarks are now in order; 

(1) The velocity vector Z* of the incoming flow in A is determined 
by the strength of the incident shock I and the orientation* of I and the wall. 
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Figure 1# 
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Indeed, the former determines the velocity of the flow normal to I on the A 
side. We know that Z* has the direction of ££ , and so the statement about the 
normal component amounts to jz*| !});Vcorapleting the determination 

of Z*. 


(2) The four discontinuity lines I, R, M, D divide the field into 
four sectors A, B, C u , Cg . Assuming continuity within each one of these 
sections, we may even treat theofas domains of constant conditions, since we 
are investigating an immediate (infinitesimal) neighborhood of T only. Thus 
we have a situation in which straight shocks delimit domains of constant 
conditions. The resulting problem is therefore again only one of applying the 
Rankine-Hugoniot. conditions (and those of a Slipstream) , involving no 
differential equations. 

Considering the importance of this conclusion, it is essential to re- 
emphasise the aasuoptioa on which it is based | continuity in each section 
A, B, at T. Let us see what the baai^ of this assumption is* 

The main reason is the experimentally established aspect of Figure 37 
which shews the lines I, R, M, D, and no others. Renee there are certainly 
no lines of discontinuity across either sector A, B, C u , Cg . However, this is 
not the only possible type of discontinuity. Thus a supersonic flow around a 
convex corner "turns" in a manner, where the state at each point is a function 
of the direction from the corner to that point. In this way there is a dis- 
continuity at the corner, but nowhere else. This Meyer discontinuity or 
angular discontinuity was discovered by Prandtl^) and Meyer^). 


my 

my 


Phys. Zeitschrift, Bd. fe, p. 23 (1907). 


Forschongsarbelten auf . dam Gebiete des Inge nieurwe sens, V.D.I. Heft 62 

( 1908 ). 
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There are reasons which make the appearance of an angular discontinuity 
under the conditions prevailing at T not at all improbable. Without discussing 
them, we point out this: the angular discontinuity is unlikely in the sector 
(^because Jhe flew there is in many cases subsonic It is unlikely in A, 

TB7 — — 

cf. discussion of Prandtl-Meyer angular discontinuity in article by 

Taylor and MacColl, previously cited (footnote 4). 

because A is ahead of the incident shock I, which is the "signal** of the 
approaching disturbance, anti therefore A ought to be entirely undisturbed. 

* 

It is unlikely in B, be cause B (while behind I) is ahead of the reflected shock 
R, which is the "signal" of I having run into an obstacle, and therefore B 
ought to be unaffected by any reflection phenomena. Hence the angular dis- 
continuity, if there is one at all, should be in C u . 

We shall point out that the assumption of continuity at T is in »ny cases 
untenable, because the conclusion conflicts with experience. We shall point 
out further that an angular rarefaction in C u does not appear to be able to 
resolve the difficulty. It may be, that there is one in B, i.e. that R is not 
the "first signal" of reflection. (In some, but not in all caees even A and 
may be questioned.) Alternatively, existing theory does not allow us to rule 
out entirely the possibility of point-discontinuities at T which are not angular 
discontinuities. The situation is obscure. However, in order to understand 
the situation and its difficulties, we must first follow up the assu^>tion of 
continuity at T, following the scheme pictured in Figure 38. 

The general procedure to be followed would be this: Knowing the strength 
of the incident shock I and the orientations of all five lines I, R, M, D and 
j It , we can follow Figures 37 and 38 and apply the conditions of Rankine-Hugoniot 
and those of a slipstream. The former determine by themselves the situation in 
C u and in , i.e. on both sides of D. The slipstream condition then requires 
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that the substance velocity on both sides of D must be parallel to f) and that 
the pressure must be the same on both sides of D. Thus we obtain three equations 
If the orientation of/j? is unknown, we can eliminate it, and still have 
two equations. If the strength of I is unknown (but the "undisturbed 1 * state in 
A - p,v, but not the velocity Z* - is known), we can eliminate it, too, and 
still have one equation. 

We call the solutions for the immediate neighborhood of T - according to 
Figure 38, and with the assumption of continuity - the three shock solutions. 

Thus measurements made on a shadow photograph of the type of Figure 32 
provide data which can be fitted to a three-shock solution only if they fulfill 
one equation. If the strength of the incident shock I is known, they must even 
fulfill two equations. Thus a determination of all three-shock solutions allows 
for a direct empirical verification. 

The best specific instance for such a determination is that one of air at 

moderate pressures (1 to 5 atmospheres), i.e. of an ideal gas with i = 1.1*. 

•** ' 

Actually the algebra of this case is rather cumbersome, but a numerical approach 
is practicable. A numerical survey of three-shock solutions was obtained by 
S. Chandrasekhar k. Friedrichs^) , and H. Polachek and^t^. Seeger^^). 

TO 

S. Chandrasekhar "On the Conditions for the Existence of Three Shock 
Waves", Ballistic Research Laboratory, Aberdeen, Report No. 36? (1943). 

TO 7 

K. Friedrichs, "Remarks on the Mach Effect", Div. 8 and Applied Math. 

Panel, NDRC (1943). 

TO 

H. Polachek and R. J. Seeger: Reg. Reflect, of Shocks in Ideal Gases 
AM-524 (1944) . 

i - • 

Experiments show disagreements with three-shock solutions. 

Thus a conflict between theory and experiment exists at the point T which 
seem3 to justify our dropping the assumption of continuity. 
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In dropping th« assumption of continuity we must try to introduce point 
discontinuities at T. As pointed out an angular rarefaction in C u would seem 
to be the most natural solution, but preliminary investigations indicate that 
this device produces no solution - not even for a weak shock in an ideal gas* 

Thus far we have been discussing an infinitesimal region of the shock 
configuration due to the reflection of a plane shock from a plane surface in 
the range of irregular reflection. This was a discussion of the boundary 
conditions designed to indicate the direction the solution of the partial 

V 

differential equation of compressible fluid flow would take. The essential 
result is that the nature of the boundary conditions is not yet understood. 

Our real interest, however, is not, except in an exploratory sense in the 
case of an irregularly reflected plane wave but in the irregular reflection of 
a spherical wave. The discussion for a regular reflection is equally valid in 
both cases. The theory of the irregular reflection phenomenon is even less 
understood in this more complicated spherical case. We are therefore forced 
in out present state of knowledge to rely on certain qualitative notions about 
the Mach effect coupled with experimental information. The experimental observa- 
tions on the reflection of spherical blast waves from a plane surface are discussed 
in the next section. 

10.4 EXPERIMENTAL DETERMINATION OF THE HEIGHT OF BURST 

The data on which the following discussion is based was obtained from TNT 
explosions, not nuclear explosions. Therefore, inasmuch as the free air pressure 
distance curve aftep taking scaling into account differs in these two 

classes of explosions, the results discussed below are not applicable to the 
blast produced by nuclear bombs. However, it is found that in the range of 
interest, i.e. 5 to 10 pounds per square inch overpressure a suitable value of 
W can be found to make the ratio of the peak overpressure produced in the nuclear 


APPROVED TOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 


* * oa 


explosion to the overpressure produced in the •equivalent* TNT explosion nearly 

t; *:'■ • 

unity. We will be interested in peek pressure only, since this date mines the 

, * 

danage caused by large bosfca (of 0 Section 10.1). For TNT explosions the free air' 
pressure, p (pounds per square inch), in the range of interest is given by 


P S3 -jOg. 




r.r 


(33) 

for W pounds of TNT 
r m distance froa explosion feet 


Nuclear explosions as calculated on the IBM have a free air pressure 
distance dependence which can be represented in the region 5 to 10 pounds per 
square inch to within 7 per cent toy 


p * 


C 

(r/W ’ /l ) ^ 



( 84 ) 


If we choose o m 206, then K is the •equivalent* tannage for the pressure range 
5 to 10 pounds per square inch. There is no equivalent tonnage in the sense that 
the press ure- di sta n e* curves for TNT and a nuclear explosion oan be made everywhere 
Identical by the choice of an appropriate value of energy released in the fora 
of blast. With this in mind we, will now proceed to construct tables of heights 
of burst by scaling experimental data on. the reflection of blast waves due to 
TNT explosions. The following 41scusslca refers to 1 pound of TNT* lengths end 


-durations are greater for an explosion of W pounds of TNT by a factor 

10 f 4~1 Determination of Triple Point Trajectory^ 17 ^ (See Figtre 89) 

,l f * ‘ Vw 11 ’ " * " ~ r ' ~ "■ *• 

The e zperiasntal data suanarlsed by Halverson Is used here, QSRD Report 4899. 
The rssilir is also referred to A. H. Taub, QSRD Report 6660 • The Effect 
k sf Mifcfc of Burst on the Blast fro* Explosives* Confidential, 1946, and 
> OBI©- Report (943 •Airtmrst froa Blast Baabs* Syaposiua, 1945 0 


In the region of regular reflection D<B„ («C 4oc extreme), the sphericity 

w 

of the shook front introduces no deviations from tha results obtained far 
plane waves because of the local character of the reflection phenomenon,. 


APPROVED FOR PUBLIC RELEASE 


APPROVED FOR PUBLIC RELEASE 


X - 80 


Figure 59 


Spfterieal shook reflected from rigid mall 



, given by theory of regular 
reflection 
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Therefore, we can use regular reflection theory and shock tube experiments 
to determine the dependence of D 0 on the height of burst (Figure 40) . 

For the case of irregular reflection the situation was studied experi- 
mentally. In particular, the triple point was located as follows: gages 
were placed at various heights (H *■ Y) and a fixed horizontal distance, D, 
from the explosion and the differences in the time of arrival between the 
direct (I) and reflected (R) shocks were noted. An extrapolation of this time 

f 

to zero gave the height of the^ triple point for each height of burst. By 
repeating the procedure for various heights of burst and then scaling the re- 
sults down to 1 pound of TNT curves of Y versus h c for various values of D 
were obtained (Figure 41). 

It is more convenient to have a plot of the height of burst he as a 
function of the horizontal distance from the explosion for various selected 
stem heights y (Figure 42). The stem heights were so chosen that they scaled 
up to 30 and 100 feet for various tonnages of TNT. 

10,4-2 Optimum Height of Burst for Peak Pressure 

Knowing the geometry of the liach effect, the next problem is to connect 
it with a choice of the heights of burst which maximize the area over which 
1 the pressure exceeds a chosen set of values. 

Suppose we are interested in the height of burst, be , which makes the 
peak overpressure, p, on the ground, occur at the greatest horizontal projected 
distance, D. The procedure would be to measure p (D) for each value of h c and 
then to plot D(h c ) for selected values of p. Such experiments have been per- 
formed using TNT^^ and the results D(h c ) for selected values of p, are 

m 

Various reports of Division 2 NDRC, by A. H. Tfcub and W. T. Read. 

1 — . . . - - 

incorporated in Figure 42. That h c , for which D a Dmav . is then the height 

of burst which yields the greatest distance to a point on the ground for which 
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Figure 40 

Theoretical limits for regtHar reflection versus charge height 
(» *1 lb TNT) (after Halverson) 
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Figure 41 

Stem height versus har is ont*l distance for a given charge height(h ) 
(I c 1 lb TNT) after Halverson 
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Figure 42 

y 

Haight of burst versus horizontal distance for a given 
stem height and overpressure 

- 1 b 1 lb TNT 

D * horizontal projected distance (feet) 
h c « charge height (feet) 
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the peak pressure has a prescribed value. 


^^^Ihis value will. 


in general, be 


7I9] ' 

It might be stated that a better scheme now exists than that which was 
used for obtaining the pressure (on the ground) versus distance curves 
for an elevated nuclear bomb. The idea is /jo \firstt deduce the reflection 
coefficients for chosen overpressures and angles of incidence from the 
measurements made on TNT and then apply these results to the free air 
blast curve from the nuclear bomb. The main reason this was not dept in 
determining the height of burst is that the IBM runs which give the free 
air overpressure, distance curves for the nuclear bomb were not yet 
carried out when the height of burst tables were made up. 

in the region of Mach reflection. If we choose, instead, to determine the 
height of burst by requiring that the stem of the Mach Y have a prescribed 
height, y, at a given peak pressure, then for this value of h c it is, in general, 
true that D< D mny . The advantage gained by basing h c on y is that the pressure 
is increased, not only on the ground but over a vertical region coinciding 
with the Mach I as well. In this way, the average pressure exerted by the blast 
on a structure is increased, resulting in increased destruction in regions 
where the pressure is marginal. The problem is somewhat complicated by the 
variation of pressure along the stem of the Y. A 15 to 25 per cent decrease 
in pressure occurs in traversing the stem of the Y from the ground to the triple 
point. Because of this variation the mean pressure along a chosen vertical 
strip is not rigorously maximized by making the stem of the Y just tall enough 
to cover it. As a working approximation, however, we will choose the height 
of burst so as to achieve a desired stem hpight at a chosen peak overpressure. 

By using a scale factor, T&bles 10.4-2 (a), 10.4-2 (b) were prepared for ' 

various tonnages of TNT. These tables give the heights of burst necessary 
to obtain stem heights of 30 and 100 feet at various chosen peak overpressures. 
The distances at which these stem heights are obtained are also listed. 
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Tafcle 10.4-2 (*> 
Y - 30 ft. 


Height of Buret end Radius at which Stem of Mach Y * 30 ft* and 
Overpressure Exceeds Given Values for Various Xilotons of TN? 



I 
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Table 10.4-2 (b) 

Y « 100 Ft. 

Height of Burst end Radius at which Stem of Mach Y ■ 100 ft* and 
Overpressure Exoeeda Given Values for Varioua Kilotone of TUT 


1 

1 P (p«i) 

W - 0.6 

(Kilotons) 

W « 1.0 

— ■ * 

^ W * 2.0 

W » 6.0 

- i 

W * 10 

vr - 20 

7/ « 100 

1 18 h c 

(ft) | 

160 

200 

300 

500 ] 

800 ; 

1200 

. 

2300 

1 d 

(ft) . 

600 

800 

1100 

1500 ! 

1 

2000 

2700 

4900 

;h h 0 


200 

300 

460 

700 ; 

1 

lioo : 

1600 

2900 

! d 


800 

1000 

1400 ; 1900 j 

-v* - 4 - 1 

65o ; looo ! 

2600 

3300 

5900 

j 10 h c 


300 

460 

1400 

1900 

5700 

j d 


1100 

1400 

1800 

2600 1 

3300 

4300 

7600 

1 9 h c 


350 

600 

760 

1100 1 

j 

1600 

2100 

3600 

d 

j 

1200 

1600 

2000 

2800 

3600 

4700 

8200 

| 8h c 


400 

660 

800 

1200 

1700 

2300 

4300 

d 


1300 

1700 

2200 

3100 

4000 

6200 

9000 

r f h c 


450 

600 

900 

1360 

•* 

19Q0 

2500 

4600 

d 


1400 

1900 

2500 

3400 

4400 

6 700 

10000 

6 h c 


600 

700 

1000 * 

1600 

2000 

2800 

6200 

d 


1600 

2100 

2800 

3900 j 

5000 

6500 

11600 

| 5 h o 

i 

650 

.800 

1160 

1700 

2300 

3100 

— 

d 


1900 

2600 

3300 

4600 

6000 

7800 

— 

4l h c 

i 

600 

900 

1250 

. 

1900 

— 

— 

— 

d 

i 

1 

i 

j. 

2000 

.. 

2800 

3700 

1 

-J 

6200 

i 

i 

L 

— 

— — 

— — 
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10.S CONCLUSION: THE HEIGHT OF BURST 

In this concluding section we will bring the amterial discussed in the 
preceding four sections to bear on the problem of determining the height of 

burst which results in the greatest area of blast damage. There are two 

* 

arguments which favor an air burst quite apart from the influence of oblique 

% 

reflection. First, a bomb burst close to the ground is accompanied by 
cratering and melting of the ground and hence a loss of energy to the blast. 
Second, an air burst avoids much shielding of one structure by another. An 
undesirable feature of air burst is, of course, the fact that the bomb is 
further removed from the target than it would be if it were burst on the 
ground. A compensating feature is the fact that the high pressure region of 
a bomb burst on or close to the ground would over-destroy the target in the 
near vicinity of the bomb. This local overdestruction represents an unnecessary 
expenditure of energy on nearby parts of the target region which decreases the 
destruction inflicted on more remote structures. The reduction in blast 
pressure due to elevating the bomb is of course more eerious for parts of the 
target which were in immediate contact with the ground burst bomb — they 
become removed by at least the height of burst. For more distant parts of J the 
target the effect of raising the boafe, off the ground is less important, and at 
distances which are two or three times greater than the height of burst the 

change in distamse from bomb to target due to elevating the bomb is completely 

•u- 

unimportant in its effect on the pressure at the target. 

Judging frost the results obtained in the low burst (100 feet) at Trinity, 
it is possible to set reasonable lower limits on the height of burst required 
to minimise some of, the above blast reducing effects due to the proximity of 
the ground . 

If it is desired to avoid fusing earth and structural materials, then 
since the radius of the ardi over which the earth was fused at Trinity was 
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about 1,000 feet, the height of burst hf which will avoid such fusing is 

* 

hf > 1,000 feet'" (35) 

This number ia for an energy release in the form of radiant energy of 3 kilo- 
tons of TUT and, since one may use an inverse square law for such radiative 
effects, thd height of burst which will avoid fusing will be related to the 
tonnage released as radiant energy f r (kilotons BUT) by the inequality 

hf > 1,000 1/3 (36) 

» 

This calculation assumes no attenuation of the bean due to absorption. It 
is not possible to state what proportion of the nuclear energy will be released 
as radiant energy without knowing the design details of the bomb. To date no 
such calculation has been carried out because of the extreme complexity of 
the problem. However, the Trinity figures give a useful indication of the 
proportion of energy that appears as radiation. 

(37) 

-^>3/4 \A £ (38> 

where hf s height of burst in feet to avoid fusing, 

Wb = blast energy in kilotons TUT. 

The available evidence on cratering from air burst bombs is very frag- 
mentary. Indeed,, because of the extremely high pressures and great duration 
* * 

of the blast from a nuclear explosion it is not possible in our present state 
of knowledge to interpret, in any complete way, data on cratering from ordinary 
explosives so 'that it will apply to nuclear explosives. The only data on 
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As a rough ruld then, to avoid fusing 
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cratering by a nuclear explosion is that obtained at Trinity where the boat 
was detonated at a height of 100 feet; A compression crater 10 feet deep in 
the center and f>00 feet in radius was formed in the close packed sand of the 
New Mexico desert. Using this point and the data in the Weapons Manual as a 
guide, it appears extremely unlikely that any cratering at all would have 
occurred had the charge been detonated at a height of 250 feet. 

h n6 crater > 250 feet for Wb « 10 kilo tens 


It is not poesible to match everywhere the blast from a nuclear explosion 
by the blast from a suitable quantity of TNT. At small distances the 
pressures developed in the nuclear explosion greatly exceed any pressures 
developed in a chemical explosion. In addition, the nuclear explosion is 
very much more rapid than a chemical explosion and does not feature after- 
burning of the constitutents so characteristic of the latter. Because of 
this, and the finite size of the mass of TNT as opposed to that of the 
nuclear explosive, the shape of the blast wave anji hence its decay as it 
travels outward is also different in the two cases. Despite this, it is 
possible to find (cf. Section 10.4) a quantity of TNT which is equivalent 
to the nuclear explosive in the sense that the peak overpressure, distance 
characteristic is nearly the same over a small range of overpressure, say 
from 5 to 10 pounds per square inch. The Trinity value Wb s 10 kilotons 
TNT equivalent is for distances where the overpressure is in the range 
5 to 20 pounds per square inch. Since the ground was close, the energy 
effective in producing a crater was greater than 10 kilotons and hence 
the assumption of 10 kilotons amounts to saying that the ground is easier 
to crater. The value for h n>c# is therefore probably too high. 


scaling for an explosion of blast tonnage Wb (kilotons TNT) 

1/3 

h n .c.> 120 W b ^ (39) 

where h n>c# - height of burst in feet to avoid cratering. 

From the above considerations the height of burst h^n to minimize the 

reduction in blast due to the proximity of the ground may be estimated as 

, • \ 

n > h n#c , or hf, whichever is greater. (40) 

The question of the value of h required to minimize overdestruction of the 
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target is sensitively dependent qp £h« jetmile of the target. Given the 
pressure which is considered as the limit beyond which overdestructicn sets 
in, a reasonable value of the minimum height of burst can be obtained from the 
pressure distance curve in free, air (cf. Chapter 7 of this volume)* and the 
multiplication of pressure on reflection from the ground, considering the blast 
wave to be normally incident on a rigid ground. 

Vow let us consider the effect of reflection on the pressure in the 
blast from a bomb burst high (h > h c ) in the air. Directly under the bocfc 
a reflection from the ground partly compensates for the loss in overpressure 
due to the increase in distance from the boipb to the target area which ac- 
companies air burst. the gain in overpressure occasioned by the reflection 
of a normally incident shock is a factor which would be 2 if the shock were 
weak, and between 2 and 8 if the shock is of finite strength (cf. Figure 21, 
Section 10.2). For shock strength in the interesting region, 5 to 10 pounds 
per square inch* thie factor is only a little above 2. This, then, is the 
effect of head-on reflection. As one departs from the point immediately under 
the bomb, the increase in the Overpressure g^fcs even more favorable because 
of the properties of oblique reflection mentioned previously. The highest 
amplification occurs'soon after Mach reflection sets in. After this it drops 
again as incidence becomes more and more glancing. Since the blast decays 
with distance and the free air peak overpressure drops, it is clearly most 


advantageous to get the greatest boosting factor where the blast pressure is 

Just marginal for the desired type of damage. One should , therefore, choose 

•' % 

the height of burst so that the maximum amplification occurs at that point. 
Since the optimum amplification occurs for early Mach reflection the 


height of burst is to be determined by the requirement that Mach reflection 
sets in at about the limit of B damage. At this point the amplification 


If it is desired to maximize A ciamage’/’the height of burst should be 
modified so that Mach ref lectidAn dSfidrtV lornespohdingly earlier. 
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factor for overpressure can bejkslbiijgh a* 3*(tifS Section 10.3), but it seems 
safer to count on a somewhat smaller value and consider higher values of peak 
overpressure, because of variowm dissipative mechanisms (cf. Section 10.1). 

The method by which the optimum Salification is properly positioned requires 

some further discussion. Actually, when the target, for example the wall of 

* 

a house, is struck, it receives two blows if there is regular reflection by 
the ground in its vicinity: on# by the direct and one by the reflected blast 
wave. If these two waves are close together- they both act as one blast. If 
they are far apart, i.e. , the angle of reflection is far from 90°, then these 
two shock waves hit the larger part of the wall with a considerable lag 
between their times of arrival. (As long as the reflection is regular the 
two shocks would arrive simultaneously at the ground but would be separate at 

t 

all points above the ground* The separation between shocks increases with 
distance off the ground.) ,In this case dissipative and other unfavorable 
effects may act between the tjk> shocks. Clearly, the beat situation from this 


point of view is one in which the two Shocks are merged together, which happens 
in the stem of the Y. Hence the house should be hit by this stem. As was 
pointed out, it is desirable to have the Mach effect in its early stages just 
as the distance at which the type of damage under consideration ceases. Now 
if one wishes tc destroy a wall, then the stem of the Y should cover the entire 

t 

wall. Consequently the height of burst should be so chosen that the stem is A 

i A r 

about as high, or higher than the target, about at the distance B damage ceases 
In view of these considerations and the validity of the peak pressure 
criterion (cf. Section 10.1) we determine the proper height of burst as 
follows: 


(1) The peak pressure required to inflict a given type of damage, 

B damage, for example, ia known: from 3 to 9 pounds per square 

• ••• • ••• •• 

. , • • • * Z Z l z 

inch. 
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(2) An amplifies tlon*JLnjtfcg» ov«r£*«4s&re of at least 2 can be 


expected 


( 22 ) 


WY 


This is approximate; the actual amplification is determined experiment- 
ally, ae described in Section 10. A* 


(3) The stem of the Mach Y should be about equal to the target 
height. 

By 1, 2 the bomb oust be buret at a distance where the peak preesure is 
something less than half of the peak pressure required according to the above. 
The height of burst is determined by 1, 2 in conjunction with requirement 3 
according to the procedure described in Section 3.0. 4. 

We can now proceed to a statement of the helghi of burst and the damage 

radius in a few selected cases. On the basis of the Trinity test and the 

■ * %> 

combat drop at Nagasaki we may assume an equivalent blast tonnage of 10 and 
20 kilotons. The peak pressures required for B damage may vary according to 
the structures involved, and even for typical residential property from country 
to country. In Japan 3 pounds per square inch may be critical; in the U.S.A., 
England, and Germany, 6 to 9 pounds per square inch might be required. If one 
is conservative one may use a higher value. (23) 


TOT 


In connection with the planning of this project, values of between 3 and 
6 pounds per square inch were usually talked about. The height of bursty 
was finally determined on the basis of * 6 pounds per square Inch estimate, 
partly in order to be conservative and partly because we were not certairf 
what the tonnage of the blast would be. The uncertainty in tennage was 
due to the probability of predetonation and general uncertainties neces- 
sarily affecting the first trials. 


The following table gives the heights of burst consistent with the pressure 
levels of 10 and 6 pounds per square inch and stem heights of 30 and 100 feet. 
The table is obtained from experiments as discussed in Section 10.4. Table 10.5 
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gives the heights of burst fon 10 Sand 20 IkfiLotbrts TNT which should be used if 

• ••• * im ee 

two heights (30 and 100 feet) are to be attained by the stem of the Uach Y at 
the positions, R, where the peak overpressures are 10 and 6 pounds per square 
inch. A mere complete table can be found in Section 10.4. 

■Bible 10.5 

SjiiSJsiMflaa 


Overpressure 
Pr (psi) 

Y . 30 ft 

Distance 

R (ft) 

. Y « 100 ft. 

Height of burst 

H <ft> 

Y • 30 ft. Y a 100 ft. 

10 

3600 

3300 

2000 

1400 

6 

5300 

5000 

2600 

2000 



w, “ .?P. 



10 

4600 

4300 

2500 

1900 

6 

6600 

6500 

3300 

2300 


On the baeis of such considerations it was found that a nuclear bomb in 
the 10-20 kiloton range should be burst at an altitude of 1500-3000 feet to 
maximize blast damage. Such heights of burst would, in addition, avoid both 
fusing and cratering the earth. 

10.5-1 Accuracy of Heigh t of Burst 

It i# appropriate here to insert a few remarks about the accuracy to 

« 

which it is desirable to fix the height of burst since this requirement is 
immediately reflected in the complexity of the fusing apparatus required. An 
accuracy of ± 50 feet, for ‘example, dictates a radar activated fuse whereas 

an accuracy of i 250 feet might be attained by a relatively simple clock 

, .... ........ 

mechanism which could be pre-set in the airplane; assuming that the bomb is 

. • •• • : 
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carried by airplane. “ 

There are several factors which limit the accuracy to which the bomb W 

h 

should be fused. They are; j ^2 

D X 

(1) Experimental errors in measuring triple point trajectories Q ^ 

os d 

**•4 Q 

and pressures in the blast configuration. (J 3 

(2) Lack of criteria as tc the precise height the stem of the O S 

Mach Y should attain at the limiting pressure which it is Q 

rr a 

. . desired to enhance by irregular reflection. 

i f 5 . 

‘ (3) Zwparfact knowledge of the target configuration and hence of 

•— J . the pressure which is required to produce marginal destruction. 

U 

SS3 (4} i The variable nuclear and hence blast performance of the bomb. 

• It can be stated with reasonable certainty, say with a prob- 
. 4 ability of 0.95, that the nuclear efficiency will he greater 


than 1/2 the rated efficiency. 


| I 

U «•» 


This was calculated by R. F. Christy for a Christy type gadget, but is 
not sensitively dependent on the specific implosion design. 


One feature which makes the choice of the heights of burst less seriously 
dependent on the above facts is the relative insensitivity of the value of 
the area over Which the pressure exceeds a certain prescribed value to the 
height at which the bomb is burst. A variation in height of burst of ± 450 
feet, for secampls, produces no more than a 20 per cent variation in the area 
defined above in the pressure region 5 to 10 poundf^per square inch for blast 
tonnages of 10 to 20 kilotons. In view of this fact, and those cited above, a 
very contwrvativo limit cm the accuracy to which the. bomb should be burst when 


used in area attack is t 250 feet. 
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